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HAUPTAUFSÄTZE 


Die turbulente Strömung an welligen Wänden. 
Von Heinz Motzfeld in Göttingen. 


1. Einleitung. Für die Frage nach der Entstehung von Wasserwellen durch den Wind 
ist es von Bedeutung, die Verteilung der Windgeschwindigkeiten über dem welligen Wasser- 
spiegel zu kennen, ebenso den Widerstand, den eine solche Fläche dem Winde entgegensetzt, 
oder auch die Druckverteilung, die der Wind an der welligen Oberfläche hervorruft. Es 
scheint nieht sicher, die Ergebnisse der turbulenten Strömung an glatten und rauhen Wänden, 
wie sie sich aus theoretischen Überlegungen von Prandtl') und von v. Käarmän?) und aus 
Messungen im glatten und rauhen Kreisrohr von Nikuradse?) ergeben, auf die wellige 
Wand zu übertragen, da nach Hopf’) und Fromm’) zwischen der „Wandwelligkeit* und der 
„Wandrauhigkeit“ scharf zu unterscheiden ist. Im folgenden soll nun über Strömungsversuche 
berichtet werden, die an festen Wellenmodellen durchgeführt wurden. 

Die Anregung zu dieser Arbeit stammt von Herrn Prof. Dr. L. Prandtl, dem ich auch für 
die vielen weiteren Hinweise bei der Durchführung der Versuche zu Dank verpflichtet bin. 

Die Messungen wurden im Kaiser Wilhelm-Institut für Strömungsforschung in Göttingen 
durchgeführt, wo auch die Apparatur hergestellt wurde. 


2. Versuchseinrichtung. Das den Luftstrom erzeugende zweistufige Axialgebläse wurde 
von einem Elektromotor von 1O PS angetrieben, dessen Tourenzahl mittels eines veränderlichen 
Feldwiderstandes gut reguliert werden konnte. Die Luft gelangte aus dem Gebläse durch 
einen Diffusor und einen Gleichrichter in eine Beruhigungsstrecke mit eingebautem Druckregler 
Schrenkscher Bauart®). Hieran schloß sich ein Übergangsstück an, welches von dem 
kreisförmigen (@uerschnitt zu einem rechteckigen Kanal von der Höhe 200 mm, der Breite 
600 mm und der Länge 6000 mm überleitete (Abb. 1). 

Bei der Prüfung der Strömung am Ende des Kanals zeigte sich eine unregelmäßige 
Geschwindigkeitsverteilung. Da sich ergab, daß die Störungen durch ungleiehmäßigen 
urbulenzbeginn in der Gleitschicht beim Eintritt in den Kanal entstanden, wurde, um ein 

1), L,. Prandtl: Ausgebildete Turbulenz. Verhandl. II. Intern. Kongr. für Techn. Mechanik, Zürich 1926. 

?2), Th. v. Kärmän: Mechanische Ähnlichkeit und Turbulenz. Nachr. d. Ges. d. Wiss. Göttingen, 1930, 8. 58. 

3) J. Nikuradse: Gesetzmäßigkeiten der turbulenten Strömung in glatten Rohren. VDI-Forschungsheft Nr. 356, 
serlin 1932. Strömungsgesetze in rauhen Rohren. VDI-Forschungsheft Nr. 361, Berlin 1933. 

4%, J,.. Hopf: Die Messung der hydraulischen Rauhigkeit. ZAMM, Bd. 3, 1923. 


5) K. Fromm: Strömungswiderstand in rauhen Rohren. ZAMM, Bd. 3, 1923. 
‘,0O.Schrenk: Ing.-Arch. I (1030), S. 350. 
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normales Geschwindigkeitsprofil zu bekommen (Abb. 2), am Anfang des Kanals an den beiden 
lanzen Reehteekseiten ein Einlauf angebracht, wie ihn untenstehende Skizze (Abb. 3) im 
(Ouersehnitt veranschaulicht. Dieser Einlauf wurde derart hergestellt, daß zuerst ein Rund- 
stab von 6 mm Durehmesser einzelötet wurde. An diesen wurde ein zweiter von 2 mm 
Durchmesser angelötet und dann der Nachlauf hinter den Stäben mit Plastılin ausgefüllt. 


Durch die an diesem Einlauf erzeuzrten Wirbel wurde ein einheitlicher Turbulenzbeginn erreicht. 
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\bb. 2. Geschwindigkeitsverteilung im rechteckigen Kanal. Abb. 3. Einlauf. 


Die Meßstrecke (Abb. 4) bestand aus einem rechteckigen an der Unterseite offenen 
Holzkasten von der Höhe und Breite des Kanals und einer Länge von 1200 mm. Dieser 
Kasten wurde am Ende des Kanals mit diesem fest und dieht verbunden. Als Diehtung 
dienten eineelegte Filzstreifen. Um das Einstellen der Drucksonden zu erleichtern, waren an 
den Seitenwänden Glasfenster angebracht. An der Unterseite trug der Kasten am vorderen 
und hinteren Ende je eine Versteifungsleiste von 150 mm Breite. In die bleibende Öffnune 
von ISOOX6O0 mm wurden jeweilig die Wellenmodelle eingesetzt und mit dem Kasten ver- 
schraubt. Vorhandene Undiehtiekeiten wurden mit Plastıilin verkittet. 


In 


N U 
| | 
| 
ij 
ws | 
— e— N 
” u z ı ! 
| I AA I 
_ U N 
2: „B UN 
ed | X 2 SITES 
| i , 
bl us - cn ® Min Wim WE mem mm elf mn 2 > - zu _ mn min _ ehem fm mr „4 | ne — ee m | un 
% Fu 
rs < ; | VE DR 
x 4) v 
„u Pr Pe A 
= — _ En 
L | = ar \ 
1 “ 


Abb. 4. Meßstreeke mit Einstellvorriehtung für die Sonden. 


Die obere Wand des Meßkastens war in der Mitte mit einem in Strömungsrichtung 
verlaufenden Spalt versehen, dureh den sich die Sonden einführen ließen. Während der 
Messungen verschlossen Sperrholzbrettehen das jeweilig offene Stück desselben, die mittels 
einer Messingschiene fest angedrückt wurden. Parallel zu dem Spalt verlief in der Mitte 
des Kastens, auf zwei @uerbalken ruhend, eine optische Schiene, an deren Schlitten die 
Kinstellvorriehtung für die Drucksonden angebracht war. Durch diese Art der Anordnung 
wurde ein Schwingen der Sonden während der Messungen vermieden. 
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Die Einstellvorricehtung (Abb. 4) bestand aus einer Führung, in der ein Profilrohr mittels 
einer Zahnstange und einem Handrad auf und ab bewegt werden konnte. Eine Millimeterskala 
mit Nonius ermöglichte eine auf "ji mm genaue Einstellung. Am unteren sich verjüngenden 
Ende trug das Profilrohr einen um eine horizontale und zur Strömungsrichtung senkrechte 
Achse drehbaren 50 mm langen profilierten Stutzen, in dessen äußeres Ende sich die Sonden 
einschrauben ließen. Die ganze Anordnung war so getroffen, daß die Sondenöffnungen in die 
senkreehte in Strömungsrichtung verlaufende Symmetrieebene des Kanals zu liegen kamen. 
Kin Ventilschlauch stellte die Verbindung der Druckleitungen des Stutzens und des Profil- 
rohres her. 

Der Gesamtdruck 9 wurde mit einem Pitotrohr gemessen von 0,5 mm Innendurehmesser. 
welches in der Mitte leicht gewinkelt war (Abb. 5a). Da es bei zwei Modellen nicht eelang, 
mit diesem Pitotrohr in die Wellentäler hineinzukommen, wurde hier ein zweites Rohr ver- 
wandt. Dieses war rechtwinklig gebogen. Der Schenkel mit der Meßöffnung hatte eine 
Länze von 10 mm (Abb. 5b). Beide Pitotrohre hatten den Eichfaktor 1 und waren bis + 18° 
richtungsunempfindlich. In Meßquerschnitten, wo die Tangente an das Wellenmodell mit der 
Horizontalen einen größeren Winkel als 15° bildete, wurden die Drucke gemessen einmal mit 
zur Tangente parallel gestelltem Pitotrohr, ein andermal mit horizontalem Pitotrohr. Eine 
dieser Messungen lag dann immer innerhalb der Grenze der Riehtungsunempfindliehkeit, und 
zwar jeweilig der größere Wert. 


ZEN 














Abb. 5. a), b) Pitotrohre, «e) statische Sonde. 


Wegen der durch die Wellen hervorgerufenen Änderung der Strömungsriehtung wurde 
zur Bestimmung des statischen Druckes p eine linsenförmige Sonde (Abb. 5e) von S mm 
Durchmesser und I mm maximaler Stärke, die in der Mitte der gewölbten Flächen je eine 
Anbohrung von 0,6 mm Durchmesser hatte, verwandt. Die Sonde wurde derart in die Strömung 
eingeführt, daß die Linsenebene mit der senkrechten in Strömungsrichtung verlaufenden 
Symmetrieebene des Meßkastens zusammenfiel. 

Die Eichung der statischen Sonde wurde in dem Freistrahl am Ende des 6 m langen 
Kanals durchgeführt. Hierbei wurde der statische Druck im Freistrahl gleich Null angenommen. 
Die Sonde zeigte einen Prozentsatz des Staudruckes q als Unterdruck an, wobei sich ein 
kleiner Gang mit der Reynoldsschen Zahl ergab, der aber zu vernachlässigen war. 
Der Mittelwert betrug - 0,14 des Staudruckes. Zwar zeigte sich eine gewisse Abhängigkeit 
des Eichfaktors von der Art der Turbulenz, denn im Mischgebiet des Strahles stieg der 
Koeffizient absolut genommen mit Annäherung an die ruhende Luft. Eine zweite Eichung 
in einem Kreisrohr mit turbulenter Strömung, wobei der statische Druck durch Wand- 
anbohrungen bestimmt werden konnte, ergab ebenfalls den Faktor 0,14. Die Prüfung der 
Kichtungsunempfindlichkeit ergab eine Unabhängigkeit von der Anblaserichtung, die mit der 
Linsenebene einen Winkel einschließt, bis + 4° für diesen Winkel. Der mit der Sonde ge- 
messene Wert ist also 


p=p Vd,lHHg, 


das Pitotrohr ergibt 


g=ptrg, 
somit ıst 
9 pP Ä 0,4 g+p' 
4 — 1.14 une p= 1,14 ; 


Die Drucke wurden mit Prandtl-Manometern gegen den Außendruck gemessen. 
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Die Wellenmodelle wurden derart hergestellt, daß zuerst aus Eisenblech die ungefähre 

Form gebogen wurde. Dieses so vorbereitete Blech schweißte man zur Versteifung in einen 
rechteckigen Winkeleisenrahmen von 600 X 00 mm Größe, der, um ein Verziehen zu vermeiden, 
noch mit zwei Rundeisen parallel zur kurzen Rechteckseite und zwei Flacheisen mit aus- 
gefeiltem Wellenprofil parallel zur langen Rechteckseite versehen war. Dann wurde Gips 
aufgetragen und mit einer Schablone die genaue Wellenform hergestellt. Durch Lacküberzug 
und Politur erzielte man eine einwandfreie glatte Oberfläche. Bei dem zweiten, dritten und 
vierten Modell befand sieh zur Bestimmung des Druckverlaufes an der Wand in der Mitte 
senkreeht zu den Wellenkämmen, d. h. in Strömungsrichtung, ein Messingstreifen mit An- 
bohrungen, die nach unten zu Schlauchtüllen führten. 


- /:=300mm - 
- W—— - 
30 
4 
_ 70 ö Fr 
195 
= 750 
x > 20 
Abu7E =. i 
> 720° 
Abh. 6. Profile der untersnehten Wellen. } 


Es wurden Messungen an 4 Modellen durchgeführt (Abb. 6). Das erste bestand aus 
einem Zug von 3 Sinuswellen mit einer Wellenlänge 4=300 mm und einer Höhe 2a = 15mm. 
Das zweite war entsprechend, hatte aber eine Höhe 2a = 30 mm. Aus einem Zug von 6 Wellen 
bestand das dritte Modell. Diese hatten nahezu Trochoidenform. Die Wellenlänge war 
; = 150 mm, die Höhe 2a = 14,5 mm. Ebenfalls bestand das vierte Modell aus 6 Einzelwellen 
der Länze = 150mm und einer Höhe 2a=20 mm. Diese Wellen hatten scharfe Kämme. 
Der einfacheren Herstellung wegen waren es Kreisbögen, die unter einem Winkel von 120° 
aneinanderstießen. Die von Mitchell’, nach der Theorie von Stokes*) berechnete höchst- 
mögliche Wasserwelle weicht etwas von dem Kreisbogen ab. 





3. Versuchsergebnisse. a) Druckverlauf an der Wand und Widerstand der Wellen. 
Da ın der vorliegenden Arbeit das zweidimensionale Problem der turbulenten Strömung an 
wellizer Wand untersucht wird, so sei die Koordinatenebene die in Strömungsriehtung ver- 
laufende lotreehte Symmetrieebene des Kanals. Die positive »-Achse zeige in die Haupt- 
strömungsrichtung, die positive y-Achse senkrecht nach oben. Der Ursprung liege auf der 
Verlängerung des unteren Kanalbodens senkrecht unter dem vorletzten Wellenkamm. Die 
Auslenkung der Wellen aus der #-Achse als Ruhelage sei mit 7 (#) bezeichnet. 
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Abb. 7. Druekverlauf an der Wand in Abhängigkeit von x. 
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aufgetragen und in Abb. S über > Während der Druck p, bei dem ersten Modell noch 
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nahezu proportional n verläuft, ist dies bei dem zweiten und dritten Modell wegen des 
Anwachsens der Gleitschieht im Tal nicht mehr der Fall. Durch den Ablösungsvorgang hinter 


’, J. MH. Mitchell: On the highest waves in water, Phil. Mag. (5) XXXVI London. 
“») 6.6. Stokes: Theory of waves, ineluding tides. Brit. Ass. Rep. 1846, Tl. 1, S. 1. 77. 
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dem Kamm des vierten Modells nimmt der Winddruck hier einen gänzlich anderen Verlauf. 
Ähnliche Messungen wurden von Th. Stanton’) durchgeführt. Die von ihm untersuchten 
Wellen haben Sinusform, deren Amplitude groß ist, so daß die Strömung an der Leeseite der 
Welle abreißt. 


Der Druckwiderstand der Wellen wurde aus dem Druckverlauf an der Wand durch 
Integration bestimmt. Es ist für die Breite Eins der Druckwiderstand 


S 
W,=\psinads, 
0 


wenn 
S die Länge der Abwickelung der Welle, 
ds ein Element derselben und 


a der Winkel, den das Element ds mit der Hauptströmungsriehtung bildet. 


Da dssina=dn, so findet man den Druckwiderstand als den Flächeninhalt, der von der 
Kurve p (n) eingeschlossen wird. 


Wegen des Druckabfalles im Kanal sind die Inhalte nicht gleich, wenn man von Wellen- 
berg bis Wellenberg oder von Wellental bis Wellental integriert. Bei einem dem periodischen 
Drucekverlauf sich überlagernden linearen Druckabfall erhält man den wahren Druck widerstand 
dureh Mittelwertbildung aus diesen beiden Flächeninhalten. 


Für die beiden ersten Modelle wurde der Widerstand für drei verschiedene Geschwindig- 
keiten bestimmt, für das dritte und vierte Modell für je vier. Es wurde bei allen Modellen 
eine quadratische Abhängigkeit von der Reynoldsschen Zahl festgestellt. Die sich ergebenden 
Beiwerte des Druck widerstandes 


r 10 
h R d ) 
Cd: > 
OUWm 


mit 26,, = Maximalgeschwindigkeit sind aus der Tabelle I zu entnehmen. 


Tabelle 1. 








; Wı W, 
Modell Re See Re 0 Um? . ” 0 Um? - c=cqa+c; obere Wand c,® 
v ze. "ae 

1 330000 V.OOUSD 0.005375 V.VO4HU 0.003459 VOLVO 
2 330000 0.0024 0.0038 0.0062 0.00365 -0.00420 
d 165000 0.0028 0,0035 0.0063 0.0034 —0.00395 
! I5000 0,0195 0.0026 0.0221 0.00425 

| 174500 0,0194 0.0028 0.0222 0,00410 
4 209000 0,0194 0,0028 0,0222 0.00395 

| 267000 0.0195 0.00295 0.02245 0.00365 


Der Druckabfall im Kanal konnte wegen des kurzen Wellenstückes nicht genau ermittelt 
werden, so daß sich der Reibungsanteil des Widerstandes nur nahezu berechnen ließ. 

Ist der Druckabfall über eine Wellenlänge Ap, der Reibungswiderstand der oberen 
glatten Wand W,’, der der welligen W, und ist die Kanalhöhe kA, dann lautet der Impuls- 
satz für die Breite Eins 


W, Ip ” h W,„’ ar Wı . 
», Th. Stanton: Marsehall and Houghton: The growth of Waves in Water due of the Action of the Wind. 
Proc. Roy. Soc. Ser. A 137, 1932. 


| 10) Bei den Modellen I bis III ist unter “,, das Maximum des Geschwindigkeitsprofils über den Punkten 70 der 
Wellen verstanden. 
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Der Widerstand der glatten Wand wurde nach dem universellen Wandreibungsgesetz 
der turbulenten Strömung von Prandtl und v. Kärmän'') ermittelt. Hiernach gilt für die 
Verteilung der Geschwindigkeit an der Wand 


u—=9»,19r,logy-+ const mit v,= i 
- 5 O 
wenn 7 die Wanedschubspannung und o die Dichte der Luft bedeuten. Es wurden also die 
gemessenen Geschwindigkeiten «'*) über log y aufgetragen und aus der Richtung der dureh 
die Meßpunkte hindurchgelegten Geraden », bzw. r bestimmt. Die gefundenen Werte des 
Gesamtwiderstandes der Wellen sind ebenfalls aus der Tabelle 1 zu entnehmen. 
































4 Mod 2:3.Mod. 1.Mod Tabelle 2. 
T 
Un hi? co, = 
Re In 1 OU Ein 
| v > 
J | k | 
7 N 63000 0.0040 
| NN | 82000 0.00375 
; \ F 
MELLE FT | | | | 106000 0,00345 
| ‘| ı \ | | BD 110600 0.00345 
mo m oO EE . SE HT 90 Mm Br 3: 
70 129000 0.00330 
Fu 
Um 1652000 0.00315 


Abb. 8. Druekverlauf an der Wand in Abhängigkeit von n. 


Für einige Reynoldssche Zahlen wurden auch die Geschwindigkeitsverteilungen an 
der ebenen Kanalwand (ohne einzebaute Wellen) bestimmt und aus diesen der Widerstand 
der Wand ermittelt. Die gefundenen Beiwerte 


T 
0 Hy“ 


.) 


eibt Tabelle 2. 


b) Geschwindiekeitsverteilung. Die Verteilunz der Geschwindigkeit wurde bei den 
beiden ersten Wellenmodellen für drei, bei den übrigen für vier Reynoldssche Zahlen 
vemessen. Da sich die Kurven bei «imensionsloser Auftragung nur wenig voneinander 
unterscheiden, so sei jeweilig der Geschwindigkeitsverlauf bei einer mittleren Reynolds- 
schen Zahl ausgewählt. Abb. 9 zeigt diesen für das Nlache Sınusprofil. Die Strömung liegt 
überall an. Der Gradient der Geschwindiekeit dieht an der Wand ıst bei Druckabfall 
erößer als bei Druckanstieg. In Abb. 10 sehen wir dasselbe für das zweite Sinusprofil. Auch 
hier liegt «lie Strömung noch an. Im Wellental zeigt das Pitotrohr in Wandnähe, wie man 
aus den Geschwindigkeitsprofilen erkennt, zu große Mittelwerte, so daß zu vermuten ist, daß 
hier zeitlich Rückströmungen auftreten, die, da sie das Pitotrohr rückwärts anströmen, nicht 
mit ihrem wahren Wert in den Mittelwert eingehen. Auf den Wellenkämmen herrscht in 
Nähe der Wand ein starker Geschwindigkeitsanstieg, der mit wachsender Reynoldsscher 
Zahl noch zummmt. In Abb. IL ist die Geschwindiekeitsverteilung für das zweite Sinus- 
profil aufgetragen bei geänderter Geschwindiekeitsverteilung im Kanal. Die Änderung wurde 
dureh quadratische Leisten hervorgerufen, die auf der unteren Kanalwand über einer Strecke 
von 2 m vor dem Wellenstück aufgenagelt wurden. a 

Die gezeichneten Gesehwindigkeitsprofile sind die über der letzten Welle der Modelle 
gemessenen. Die vorhergehenden entsprechenden Kurven zeigen in Wandnähe etwas kleinere 
Geschwindigekeiten. Es ist nieht zu sagen, ob die bei einer größeren Anzahl von Wellen sich 
endgültig ausbildenden Geschwindigkeitsprofile von denen der Abb. 9 bis 11 noch stark 
abweichen. 


11) S, Fußnote S. 193 
Nr. 5, 1933 


17) Hierbei wurden die Gesehwindigkeiten aus den einzelnen Profilen über einem Wellenzug gemittelt. 


und auch L. Prandtl: Neuere Ergebnisse der Turbulenzforschung. Zeitsechr.d. V.D.1l. Bd. 77. 
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Abb. 10. Geschwindigkeitsverteilung 

über dem zweiten Wellenmodell 

(unten) und Verlauf des statischen 
Druckes (oben). 
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ei dem dritten Modell wurden, um einen endgültigen Zustand zu bekommen, vor dem 
Wellenstück auf dem Kanalboden Halbrundleisten im Abstand der Länge der Trochoidenwelle 
über eine Strecke von 2,7 m angebracht. Hierdurch wurde erreicht, daß, falls die sich ent- 
sprechenden Geschwindigkeitsverteilungen über den drei letzten Wellen des Modells noch 
Unterschiede aufwiesen, diese dann innerhalb der Meßgenauigkeit lagen. Abb. 12 zeigt diese 
Geschwindigkeitsverteilungen. Wegen der besseren Übersicht wurden einige gemessene 
Kurven ausgelassen. 
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Abb. 13. Geschwindigkeitsverteilung 

über dem vierten Wellenmodell 

(unten) und Verlauf des statischen 
Druckes (oben). 
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Die bei dem dritten Wellenmodell auf dem Kanalboden aufgenagelten Halbrundleisten 
wurden bei dem vierten Modell durch Dreieckleisten ersetzt, wodurch ebenfalls das Gewünschte 
erreicht wurde. Abb. 13 gibt die Geschwindigkeitsverteilung über den scharfkämmigen Wellen 
wieder. Die Strömung, die hinter dem Wellenkamm abreißt, liegt im tiefsten Punkt des Tales 
wieder an. Wegen des größeren Widerstandes der unteren Wand ist das Geschwindigkeits- 
maxımum ım Kanal zur glatten Wand hin verschoben. 

In den Abb. 9 bis 13 ist über den Geschwindigkeitsprofilen der Verlauf des statischen 
Druckes im Innern der Strömung aufgetragen. Die durch Extrapolation gefundenen Werte 
an der Wand stimmen gut mit denen der Anbohrungen überein. 
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e) Stromlinien, Gesamtdruck auf den Stromlinien und Verteilung der scheinbaren Schub- 
spannungen. Aus den Gesehwindigkeitsprofilen wurden als Kurven 


Y 

\udy=const 

0 
ie Stromlinien bestimmt, wobei « mit dem Betrag der Geschwindigkeit identifiziert wurde. 
ls ist hierin eine geringe Vernachlässigung, da man dann eigentlich normal zu den Strom- 
linien hätte integrieren müssen. Die Strömungsbilder der drei letzten Wellenmodelle sind 
in den Abb. 14 bis 16 wiedergegeben. Die darüberliegenden Kurven zeigen den Verlauf des 
(Gesamtdruckes auf den Stromlinien. Man sieht, daß hinter den Wellenkämmen auf den 
unteren Stromlinien ein Energieverlust eintritt, der dann stromabwärts weiter in das Innere 
der Strömung hineinreicht. Bei Abfall des statischen Druckes nimmt die Energie auf den 
unteren Stromlinien wieder zu. Dieser Energiegewinn stammt aus dem Innern der Strömung, 
da hier ein ständiger Gesamtdruckabfall auf den Stromlinien vorhanden ist. 
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Abb. 14. Stromlinienbild (unten) und Verlauf des Gesamt- Abb. 15. Stromlinienbild (unten) und Verlauf des Gesamt 

druekes auf den Stromlinien (oben) bei dem zweiten Wellen- druekes auf den Stromlinien (oben) bei dem dritten Wellen- 


modell. modell. 


n. - 
; d 
1 da )(ds +k, ds dn) 
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Abb. 16. Stromlinienbild (unten) und Verlauf des Gesamtdruckes Abb. 17. 


auf den Stromlinien (oben) bei dem vierten Wellenmodell. 


Die Verteilung der Schubspannung läßt sich bei turbulenter Strömung nicht aus dem 

Newtonschen Reibungsansatz = u e fy —Zähigkeit der Flüssigkeit) ermitteln, da durch 
r dy‘ „ ww | Fu; 

den Impulsaustausch der turbulenten Schwankungen noch eine scheinbare Reibung zwischen 
den einzelnen Schichten der Flüssigkeit hervorgerufen wird. Um die sogenannte scheinbare 
Schubspannung zu ermitteln, gehen wir aus von der ersten Bewegungsgleichung, die, wenn 
sie in jedem Punkt für die Richtung der durch ihn hindurchgehenden Stromlinie und der 
senkrechten Trajektorie zu den Stromlinien gilt, folgendermaßen lautet: 


09 or 
— +2tk,, 

Os on 

wenn 0 s ein Linienelement der Stromlinie, d n ein solches der Trajektorie und k, die Krümmung 

der Stromlinie bedeutet 
















202 Motzfeld, Die turbulente Strömung an welligen Wänden a a 





Zu dieser Gleichung gelangt man, indem man ein Flüssigkeitselement betrachtet, das 
von einer Stromlinie, einer Parallelkurve dazu und zwei senkrechten Trajektorien begrenzt 
wird und die Breite Eins hat. Die Komponente in der s-Richtung der Resultierenden aus 
dem statischen Druck und den Schubspannungen ist gleich dem Produkt aus der Masse und 
der Beschleunigung in der s-Riehtung des betrachteten Elementes. In Abb. 17 ist dieses 
Element mit den angreifenden Kräften dargestellt. Aus dem beigefügten Kräfteplan ergibt 
sich für die Komponente der Schubspannungen in der s-Richtung 


| Ye or OT 
th, dsdn-+ „|? +ds N .) k,dsdn+ [7 + dn ) ) (ds+k,dsdn)—rds. 
Aus dem statischen Druck resultiert in der s-Riehtung die Kraft - dsdn z m, ferner ist Masse 


A) 


u 
mal Beschleunigung gleich odsdnu Fr also gilt, wenn wir durch dnds kürzen, 
0 s | 


Io17 OT 07 OP, o du? 


sts + „dn 43 


.) 
aT k. r a wi 
on os 


— 
u“ 


II 


Os Os’ 


. ® » 0 7 ° 
lassen wir nun dn und ds>0O gehen und bedenken wir, daß p+Z «= ist, so erhalten 


wir obige Gleichung 


09 or ' 
z £ +2tk.. 
Os ON 


r . r} Ö q hd .. * 
Aus den Kurven der Abb. 14 bis 16 wurde \ bestimmt und über y aufgetragen und 
( Fr « \ 


dann hierüber integriert. Wo es nötig erschien, wurde noch eine Korrektur angebracht unter 
.) . 
. . . 1]: Thy r: . \ . i ug . 
Berücksichtigung des Gliedes — (a Neigungswinkel der Tangente an die Stromlinie). Es 
OS A ö 


sei aber «darauf hingewiesen, daß die Integralkurven wegen der Ungenauigkeiten bei der 


RR Or 
KErmittelune von 


F und A, mehr nur einen qualitativen Verlauf der Schubspannungsverteilung 
(8 g 


eeben können. Abb, 15 und 19 zeigen die Verteilung für das zweite und vierte Modell. 
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4. Vergleich der Ergebnisse mit denen der turbulenten Strömung an ebener Wand. a) Potenz- 
vesetze. Auf den Geraden „== const nimmt die Geschwindigkeit beim Fortschreiten in der 
‚-Riehtung verschiedene Werte an. Erst für größere y behält sie ihren Wert bei. Im folgen- 
den sei nun unter « der mittlere Wert der Geschwindigkeit auf den Geraden y == const ver- 
standen, Es läßt sich dann diese Geschwindigkeitsverteilung über den vier Modellen durch 
folgende Potenzen von y wiedergeben: 





Modell | 1l 111 IV 
: 48 1/2*) == 
n 1/9 . 116,5 16,5 1/2,5 


*, in Wandnähe. 
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TR Gr 
Für Reynoldssche Zahlen a 50000 fand Blasius'?) für die turbulente Strömung 


in glatten Rohr, daß « proportional zu der 1/7 Potenz des Wandabstandes ist. In Abb, 20 
sind die dimensionslosen Geschwindigkeiten «/t,„, über y/d aufgetragen (hierin bedeutet ö den 
Wert von 9, für den « sein Maximum annimmt), dazu die von G. Wüst'') 1920 über der 
()stsee gemessenen. 
b) Logarithmisches Gesetz'’). Nach den theoretischen Überlegungen von Prandtl und 
\ von v. Kärmän ist es möglich, die turbulente Geschwindigkeitsverteilung an ebenen Wänden 
wiederzugeben durch den Ausdruck 
9* 
. log y-+(), worin v, 


b} 7 “ 
U, X 0 


De 


T 


| die Schubspannungsgeschwindigkeit und z ein universeller Zahlenwert bedeuten. Die Inte- 
orationskonstante € ist noch von der Beschaffenheit der Wand und der Reynoldsschen 
g 


Zahl abhängig. 


. h 30 
Aus den Rohrmessungen von Nikuradse'*) ergab sich für z=0,4 und C=(,= log - 
[ 


| 3m, 
| für das rauhe Rohr, wenn k die Korngröße der Rauhigkeit ıst und = („= log a für das 


elatte Rohr, worin » die kinematische Zähigkeit bedeutet. 
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Abb. 20. Dimensionslose Gesehwindigkeitsverteilung. Abb. 21. Dimensionslose Geschwindigkeiten in Abhängig 
(Die Kurven sind jeweilig um 0,1 nach unten verschoben.) keit vom Logarithmus des Wandabstandes. (Die Kurven 
sind jeweilig um 0,1 nach unten verschoben ) 
Aus unseren Messungen an der ebenen Kanalwand (siehe Tabelle 2), die wohl nicht 
| vollkommen glatt angenommen werden darf, ergibt sich ein etwas anderer Wert. Es ist: 
33V, 
C,= log 2 


Bei rauher Wand bleibt die Frage noch offen, wo der Ursprung der y-Achse zwischen den 
„Höckern der Rauhigkeit“ genau zu liegen hat. In Abb. 21 sind die dimensionslosen Ge- 


> - id u .. . .. » 24 hj 
schwindigkeiten z über log (y — y,) aufgetragen, wobei ,, so gewählt wurde, daß die Gerade, 
m 


die sich durch die wandfernen Punkte legen läßt, eine Riehtungskonstante 
v, 


m—9,1d 
u m 


T 4% ) >| m N 
OU“ Um BD) 


.) 


hat, so daß 


gleich den in Tabelle 1 angegebenen Werten wird. 


13) H. Blasius: Das Ähnlichkeitsgesetz bei Reibungsvorgängen. Phys. Z. Bd. 12, 1911 oder Forseh.-Arb. 
Ing.-Wes. H. 131. 
a8) Ü. Wöüst: Verdunstung auf dem Meer. Veröff. Inst. Meereskunde. N.F.A. 6, Berlin 1920, 
15) Siehe Fußnote S. 193. 
16) Siehe Fußnote S. 193. 









































s ö N a En Ztschr. f. angew. 
204 Motzfeld, Die turbulente Strömung an welliren Wänden Math. und en 





Man erkennt, daß in Nähe der Wand das logarithmische Gesetz nicht mehr gilt, und 


I 
7 bei den 


. * . ® ..g . ( 
zwar sind die Abweichungen um so bedeutender, je größer das Maximum von 7 
ai A 


einzelnen Modellen ist. 
Für das vierte Wellenmodell, welches wir als rauhe Wand auffassen können, liegt der 
.) . 


Punkt „9, = y=V um ya höher als der Wellenkamm. Nehmen wir diesen Punkt als 


Ursprung und setzen für k=2a, so lautet die Integrationskonstante 


1.25 


C,=10g5 
Der Unterschied zu Nikuradse ist wohl darauf zurückzuführen, daß die Korngröße %k nicht 
mit den wirklichen Höhenunterschieden (2 a) identisch ist, da bei Nikuradse der aufgeklebte 
Sand noch mit Lack überzogen wurde. Weiter sind die Arten der Rauhigkeit wohl nicht 
initeinander zu vergleichen. 








“) 
’ , ' ’ a ,„ 2a 
Bei der ersten Sinuswelle liegt der Ursprung y=0 um Y= Fa bei der zweiten um 
r 2 Ad . H rg . r da .. r rm» m 
Y= , ‚bei der Trochoidenwelle um Y: g höher als das Wellental. Trägt man Y über 
dı 
"| auf, so erkennt man, daß 
wu MAX 
Vdaldn da e dn 
we tg a, oder nahezu Y=rza Ir 
IT (A max ZT max 














gesetzt werden kann, wenn a, der größte Neigungswinkel der Tangenten an die Wellen ist. 








4 T 
e ebene Platte 
o I. Modell 
N © 7. Modell 
ZI oe Made 
Gr | 





0 “ r 
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647.22 lg Km 























Abb. 22. Abhängigkeit der 
Integrationskonstante C,„ von tg Em: 
Weiter zeigt sich, daß auch die Integrationskonstanten für wellige Wände (,, die sich für 
die Wellen der 3 ersten Modelle ergeben, linear von tga, abhängen (Abb. 22), und zwar 
ereibt sıch 


) 
>)UV, 


Co» =log — IS.2tg an. 
v " 
In genügender Entfernung von der Wand lassen sich also die Geschwindigkeitsprofile 
wiederzeben dureh die Ausdrücke 


. 124 
u D1Dr,log 


2a 
für scharfkämmige Wellen, wobei die Frage nach dem Nullpunkt von y noch. offen bleibt und 


IV, 


9,19 0,108 18.2tg a, 


Y 
für Wellen, bei denen keine Ablösung des Windes eintritt. 
Aus diesen beiden Formeln findet man nun leicht Ausdrücke für die mittlere Wand- 
schubspannung. Da r=or,? ist, ergibt sieh für scharfkämmige Wellen 


\ 


ou 
T — - 


1.5-1\? > 
35 llo& 
> ») A 
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wenn 7; die Geschwindigkeit im Abstand y==! ist und für flache Wellen 


ou 
T . } > 
/ T er 
l \ 
2‘ 75] o\2 E IS.2te 
iO 1IOE\» y D,-lS Am 


woraus sieh 7 durch Iteration ermitteln läßt. 


5. Anwendung der Meßergebnisse auf das Problem der Entstehung und des Anwachsens 
von Wasserwellen bei Wind'’). Die Frage nach dem Anfachen von Wasserwellen durch den 
über die Wasseroberfläche streichenden Wind wurde zuerst von Lord Kelvin'*) unter der 
Annahme reibungsfreier Flüssigkeiten und unter Berücksichtigung der Oberflächenspannung 
behandelt. 


Setzt man für die Stromfunktion im Wasser 


Fan, ahreit+er tie rer erre 


dann verlangt die Kontinuitätsbedingung 
AP=0, 


daß 
p,=G, (l) ei** und Eier 5. er 
ist. 
Für ein Teilchen der Oberfläche ergibt sich aus Gl. (1) eine Geschwindigkeit in der 
y-Richtung 





od 2 
"y—Wu 2 yl»—o Pr (q 1 4 Br a (»). 
u a a . D ' ’ 
Weiter gilt für ein Teilchen an der Oberfläche »,—o gg (worin D die substantielle 


Differentiation ausdrücken soll), was sich auch angenähert schreiben läßt 


0 ur. 0 
" = +u u Er SA ee Ser 0 er 
YUV Vf Um dx (+), 
un ‚On a ee en 
wenn wir das Glied v \ 5 vernachlässigen und für «” die Strömungsgeschwindigkeit «„, setzen. 
( 
Hat das Wasser an der Stelle y=—h eine feste Begrenzung, so gilt hier #,—n =V 
und somit 
e"ho,—e'ho, 3:9 


Diese Beziehung liefert im Hinblick auf (2) 
ER 5 0MCE 2) 1. WR EEE EEE EEE EEE  ° 


Nun läßt sich die Potentialfunktion unter Verwendung der Gl. (3) bis (5) durch die 
Oberflächenform 7 (x, t) ausdrücken, und zwar ist 


cosh # (y-+h) 


D' Wnd- z 
a. sinh x h 


A Umyy ae 





097 ‚eo "| 


Wegen der Kontinuitätsbedingung muß darin n= H(t) e=-** gesetzt werden. 
Entsprechend ergibt sich für das Potential der Luftbewegung, wenn wir uns die Luft 
über dem Wasserspiegel bis ins Unendliche erstreckt denken: 


e—-*"y 
D- Um X Re 





0 ON 
Ö t Un A) 1 


17) H. Thorade: Probleme der Wasserwellen. Hamburg 1931, S. 37 bis 54. 
18) W. Thomson: Hydrokinetie Solutions and Observations. Phil. Mag. 4, 1871. 
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Nach Bernoulli folgt aus (6), bei Beschränkung auf lineare Glieder für den Druck an 


der Trennungslläche 
E_. 40 n.7,St nad „,. Ant d ‚nl, | 
og eonst — eteh # I r | NW + 2uU a Ardt Um Ir: +9 at), . (8 
und aus (7) 
p | ı /0? ir, f) a 0°’ (z,8) ol (8, t)\ : 
” eonst + y | AWE Um a..: Un Ya: gynıe, 2 re VE (9) | 
oder, wenn wir n = H (te ‘** einführen, 
p’ ”. en, | u IH (t) 2: +, | En 
ons cighahi | ge ZU mirgg e m u) +4 ZAUR (10), 
I /d? Hit ‚dHdıt Sol ar 
Er . ( Zn | um xi- di Ä Um” ZA) g ZA BIER. oe A 
‘) 4 / 





Der Druckunterschied, hervorgerufen durch die Oberflächenspannung T, läßt sich an- 
genähert ausdrücken durch 


[4 ar 7 


p—-p= ir r oder p-p =—Tx?’ H{t)e-i** 


und es ergibt sich für /(t) die Differentialgleichung, wenn wir noch durch e=-‘** kürzen: 


d®H dH\ 


N . ’ ’ 
(o + 0’ eteh x h) r 2ir luno+ wo’ eteh x h 
C ( AR (Hm & 1 mx > ) At 


(ot, +0’ WW, etghzh)H - x glo — o' ctgh xh) 


“"T-H 


(12). 





Als Lösung setzen wir an 
Hau get" 


darin ist a die halbe Höhe und o die Phasengeschwindigkeit der Oberflächenwellen, die jetzt 
die Form haben 


DR, 5.5 ee ar rt Fi 
und es ergibt sich als charakteristische Gleichung: 
(0 + 0’ etgh # h)2’ 0’ - 2(u„no+t wm oetghzh)z’o + louyn to wm etghxh) 
+xg(o — o’ctghzh)=»"T. 


Nach der Auflösung nach o folgt daraus für die Phasengeschwindigkeit 





40 + um 0’ eteh x h 
== “ u 
0-0 eteh x h 
(14). 
q0o- o'ctehzh lx 00 eteh x h Ä z NE 
+ ee 4 ; Zu um — U m)“ 
\ zo+o'ctehzh o-+o'ctghzh (o+o’ctgeh “hy? " . 


Die Trennungsfläche ist instabil, wenn 


go "cteh’(zh)- 0° o+o’ctghxh 
00’ etgh x h 00 etgh x h 


(un — Um) 5 U 
Ist. Um zu der kleinsten Windgeschwindigkeit zu kommen, die noch Wellen zu erzeugen 
vermag, fragen wir nach dem Minimum der rechten Seite dieser Ungleichung. Ein solches 
trıtt ein, wenn 
0 eteh “h o : M 
AEEN  2_ 1) 
ist und dann gilt 


‚0 etehx.h+o 


1 
(Um -un)>2 -[Tg(octghh—o) - . .». .»... MM. 


O 0" etgh H h 
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Pür tiefes Wasser liefert diese Formel für die minimale Windstärke ı, 1 — u = 640 em/sec 
nd für die Länge der ersten Wellen 1,5 em. Dieses Ergebnis steht aber zu den Beob- 
achtungen (s. unten) im Widerspruch. 


Einen anderen Weg zur Bestimmung der kritischen Windgeschwindigkeit beschreitet 
Teffrevs'®). Er geht von der Tatsache aus, daß der Wind, da die Luft nicht reibungsfrei 
ist, auf der Luv- und Leeseite der Welle eine verschiedene Druckverteilung hervorruft, wo- 
durch der Welle, wie gleich gezeigt werden soll, aus dem Wind ständig Energie zugeführt 
wird. Ein Anwachsen der Wellen tritt nun ein, wenn diese zugeführte Energie größer ist 
als die dureh Reibung im Wasser verzehrte. 


Wir wollen annehmen, daß das Wasser am Boden h in Ruhe ist, d.h., daß «’ „= 0 
ist. Dann lautet Gl. (6) 


I coshı z(y+h)0ı 


„ . S 
d x sinh#zh 01 (15) 
oder mit n = acos# (ot —- x) (vgl. Gl. (13)) 
coshz(y-+h) _. 
D' - J asın lot Bun 2. BE Fi an SEES 


sınh # h 


Die dureh innere Reibung in der Zeiteinheit im Wasser verzehrte Energie berechnet 
sich nach Lamb°’) aus 


Du \\ 2 rn dF, 


worin a die Zähigkeit des Wassers bedeutet und 


vu ( D' } . ? 2 


Or Ri Y 


ist. Betrachten wir ein Volumen Wasser, welches sıch senkrecht unter einem Stück der 
Oberfläche von der Breite b und der Länge / befindet und bis zum Boden reicht, so ist hierin 
(die Dissipation 


EI ': WG 9  ; \  } 


Für den Teil der Druckverteilung an der Oberfläche, der für seine Rechnung von 
Bedeutung ist, setzt Jeffreys an 


(20). 


p* = 5 $ ( 


Die von diesem Druck in der Zeiteinheit geleistete Arbeit für ein Oberflächenstück der Breite b 
und der Länge A ist dann, da die Geschwindigkeit eines Oberflächenteilchens in der y-Riehtung 


uch EEE u 
gleich \} ist und dssina=d.x gesetzt werden kann, worin ds ein Längenelement der Ober- 
( 7 


Näche und « der Neigungswinkel desselben gegen die «-Richtung ist: 


7 
(0909 
V=Ss0 (tm ob \ ws dx 
0 
2 
sollt. oFazob\sin’z(ot a)dı 
1 
BEER N EI N rar 


»), H.Jeffreys:On the Formation of Water Wavesby Wind. Proe. Roy. Soe. Ser. A. 107 und Ser. A. 110, 1925/26. 


-+), H. Lamb: Hydrodynamies, $ 317. 
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Aus Gl. (19) und (21) ergibt sieh nun nach einer kleinen Umformung für die Anfachung von 


i 3 u NIT ia 
Wellen die Bedineune. wenn wir noch dureh b-/ kürzen und — =rv’ einführen 
ı) 


De 


(mo) tr’ o’ 


DIEBE N ee 
OÖ 50 a 


Vernachlässigen wir in Gl. (14) die Oberflächenspannung gegenüber der Schwere, was wir 
können, sobald 4> Sem bzw. 2<0,8 ist, und beschränken wir uns auf kleine Wind- 
veeschwindiekeiten, dann erhalten wir für die Phasengeschwindigkeit, da «’„=0 angenommen 
wurde, 

4 0" etgh (z h) 2) (23) 
X 0 etgh (2 h\—+ O . 


Unter Verwendung dieses Ausdruckes läßt sich (22) dann schreiben 


vg 0’ 0 eteh (# h) — 0 
etsh x I eK 


u 0)” 0 
vr s oo’ ectgh (zh)—+o 


Um zu der Windstärke zu gelangen, die gerade noch Wellen zu erzeugen vermag, suchen 


wir das Maximum der linken Seite dieser Ungleichung. Es tritt ein, wenn 0 = o, 3 Uom 
ist, und es gilt für diesen Fall 
i _v’g0' 0 eteh (2x h)— 0 E 
ER = = -Ushah . .:. .» : . 5 - Ä. 


s 0 o’ etgh (#h)-—+o 
Unter Berücksiehtireune von Gl. (23) folet hieraus 


5) ’ 
i „4 0 o' eteh(zh) - 0 4 \ 
> Aa 5a ro PN ’ tech? % h . ö ä . ö . . 26). 
PTrao’o eteh (#h)-+o . 


Kine genaue Nachprüfung dieser Formeln war bisher nicht möglich, da der Faktor s 
nicht bekannt war. Jeffreys ermittelte umgekehrt aus Beobachtungen der ersten wind- 
erzeugten Wellen, für die er eine Wellenlänge von S bis 10 em bei einer Mindestwindstärke 
von 104 bis 110 em/see feststellte, sowohl aus (25) als auch aus (26) s zu 0318 - 0,229 und 
sah hierin eine Bestätigung seiner Theorie. Gleichzeitig wird hierdurch die Vernachlässigung 
der Oberflächenspannung gerechtfertigt, so daß also die winderzeugten Wellen keine „Rippel“ 
(Thomsoni, sondern Schwerewellen sind. Bei Laboratoriumsversuchen von Stanton’') er- 
cab sich dagegen u 250 em/see und /, = 6em. Diese Werte liefern aus (25) s = 0,025 und 
aus (26) Ss = O,48. 


om 


Berechnen wir nun den Druckwiderstand einer Welle „=asinz.x für die Breite Eins 
unter der von Jeffreys gemachten Annahme der Druckverteilung 


& ( an 
"—soln 0)" f 
| "U me I dr 
2 d 7 
so Ist weren dı, Äx 
dx 
/ 
Nur soltuu- oa z\cos’zar da 
0 
oder 
Kerr, iR... ira er } 
rn W 
Für den Beiwert e,= e folgt daraus 
lt 0 A 
CA ee. ; . ; . ; R : ! N : ; . i (28). 


21) s, Fußnote S. 197. 









































sand 17, Heft 4 


uzust 1937 Motzfeld, Die turbulente Strömung an welligen Wänden 209 








Setzen wir einen mittleren der von Jeffreys gefundenen Werte von s hierin ein, so erhalten 
wir als Widerstandsbeiwert für die Sinuswellen unserer beiden ersten Modelle e4 = 0,0068 
und e7 = 0,027, also Werte, die um eine Zehnerpotenz größer sind als die gemessenen. Wir 
schen also, daß die Jeffreysschen Ungleichungen (25) und (26) leider keine befriedigende 
Lösung des Problems darstellen. 


Für die weitere Untersuchung wollen wir jetzt die kinetische Energie der Welle 
bereehnen. Nach Lamb’) ergibt sich diese aus 


| 


») 


2 N\p’ 
, ( ’ 
) \ \ «D' dAF, 
* “ Ö n 


es 
nm 


worin n ein Element der äußeren Normalen der Kontrollläche F ist. Betrachten wir als 

Kontrollfläche einen Quader, der senkrecht unter einem Stück der Oberfläche von der Breite 

» und der Länge 4 liegt und bis zum Boden reicht, so ist das Integral nur über das Stück 

der Oberfläche zu erstrecken, da sich die Beiträge der Seitenwände des Quaders herausheben 
N\cb’ 

und am Boden po V ist. Es ergibt sich also für die kinetische Energie, wenn wir PD’ aus 


(+]. (1Sa) einführen, 


| 
E vo @o’zb-.4ctehzh. 


l 


Da wie bei allen Schwingungen die potentielle Energie gleich der kinetischen ist, so ist die 
totale Energie 
| 


.) 0 da? ozhb-, eteh z . ‘ z - - : { f Ä ' j (29). 


Die von dem Wina der Welle zugeführte Energie ergibt sich aus dem Produkt von 
Kraft mal Weg. Für die Zeiteinheit und einem Stück der Oberfläche der Breite b und der 
Länge / lautet dies 


| 
() = (450 (tt, 0) 0- h».% für Mn oO 
(30). 
() Ct5oltm oo-b-/ für u.<o 





Für die weitere Rechnung wollen wir uns auf den Bereich «,„ > 6 beschränken, da wegen 
(les Energieentzugs für u, <o stets eine Dämpfung der Wellen eintritt. 

Der Widerstandsbeiwert ey ist nun, wie unsere Versuche gezeigt haben (Tabelle I), noch 
von dem Verhältnis : abhängig. Da sich diese Abhängigkeit nur empirisch ermitteln läßt, 
wollen wir angenähert setzen 

ee a Ya ee ze A a 


worin der Proportionalitätsfaktor s und der Exponent » noch aus unseren Messungen zu er- 
mittelnde Größen sind. 


Wegen der Erhaltung der Energie muß nun die zeitliche Anderung der totalen Energie 
gleich der Differenz aus der zugeführten und der verzehrten Energie sein. Dies führt zu der 
Gleichung 


da | 
’ 5) 7 4 0) . 7 °), 
era ms "a" olum—o”-otghrh—-2ua"o’ . . . . . .. (82), 
( Zu l 


wenn wir nur a mit der Zeit veränderlich annehmen und durch b-/ kürzen. 


-?) Lamb: Hydrodynamies, $ 44. 
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Als erste Lösung dieser Gleichung erhalten wir @a=0. Für a==0 können wir (32) auch 
schreiben 


da f' 4 (' 
Ad =. 
dt 1 r- 
wenn wır für 
| >) n LET o)” ‚ Dr ) ’ 5 1] 
un tehch=C, und für 2v«’=[, 
Z 7) Ö 2 
setzen. Die zweite Lösung lautet dann 
gt da | 
\ s nl + eonst 


Car? _C, f 
(n — 2) In 0 u 


(In 2)\C,-t=In -—- r 
3 BE Pak + a 


Für den Grenzfall n=2 zeht diese Gleichung über in 


und für den Grenzfall verschwindender Wasserzähigkeit, d.h. für (,=-0 aber n ==2 in 


D) „ 


2 —n)C,t=ua"—a/ ee 


0 


Aus der letzten Gleichung erkennt man, daß bei verschwindender Wasserzähigkeit die 
Oberllächenwellen angefacht werden, wenn u„>6 ist. Aus Gl. (33) ergibt sich für die An- 
fachung die Bedingung für 0-0 t,. 


C,a” 2 (' 


, * 


Daraus folet 


l. für n<2 
) 
5 My 6)” ) n Y32 
; ; tech 2 h ee ee 
Oo KV I 
und 
ya für 7) 2 
Ss Oolltıy o)” : > ad 
| re tsh x I) Re. 
Eee 4 AR, 


Die erste dieser beiden Ungleichumgen besagt, daß die Wellen dann angefacht sind, wenn 
EEE 5 a ar ,__: er = 

das Verhältnis von Amplitude zur Wellenlänge | 7 unter einer von der jeweiligen Wind- 

| - 


stärke abhängigen Grenze liegt, oder anders ausgedrückt, zu jeder Windstärke und zu jeder 
Wellenlänge gibt es eine bestimmte stationäre Amplitude. Dagegen folgt aus der zweiten 
Ungleichung, daß der Wind nur dann die Wellen anzufachen vermag, wenn ihr Verhältnis 
von Amplitude zur Wellenlänge schon einen gewissen von der jeweiligen Windstärke ab- 
hängeigen Wert überschritten hat. 


Aus Gl. (34) folgt schließlich die Instabilitätsbedineung 


3. für n > 


was natürlich zu den Jeffreysschen Ungleiehungen (25) und (26) führt, da unter der 
leffreysschen Annahme der Druckverterlung, wie Gl. (25) zeigt, der Exponent » sich zu 2 
ergibt. 
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Zur Bestimmung von n und s wollen wir nun unsere Ergebnisse der Druck widerstands- 
Ad u i d . 2 E 
messungen am ersten . 0.025) und zweiten ; 0,05) Modell zugrunde legen. Es sei aber 


darauf hingewiesen, daß die sich ergebenden Werte nicht sehr genau sein können, da zu ihrer 
Bestimmung eine größere Anzahl von Wellenformen untersucht sein müßte. Aus den zu- 
vehörieen Beiwerten e4 = 0, W)S5 bzw. ea == 0,0024 ergibt sich » zu 1,5 und s zu 0,014 und es 


_ 


eilt somit die Instabilitätsbedingung (26). Nach Einführung obiger Werte lautet diese 


= 4 l (u 0)” ’ 
(0,0035 S . teh = h) ea an er 
0 „_ Ö 


Bei vorgegebener Windstärke gehört also zu jeder Wellenlänge eine bestimmte stationäre 


dad. % { 
3 berechnet. Die 


(, 


Wellenlänge a,. Für einige Windstärken wurde das Verhältnis von 


Werte sind in der Abb. 23 über der Wellenlänge aufgetragen. Die Kurven zeigen, daß es 























0 7 30 120 7) 240 TE 77 1200 
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Abb. ?3a und >3b. Abhängigkeit der stationären Amplitude von der Wellenlänge bei verschiedenen Windstärken. 


. . r* nie . . A . r . ® *n. Y 
bei jeder Windstärke ein Maximum von Fr eibt. Um dieses zu bestimmen, führen wir o 


aus (23) ın (38) ein 


00023 0 teh Hr h 0 eteh (2 h) 4 0 r 9 Fr 
£ )») 0" y’ g 0" eteh (x h) Ö m 6) tx | | | | (39) 
| 


Hierin wird die linke Seite ein Maximum, wenn 0=0,° 5 tt, wird, d.h. das Verhältnis 
> 


wird am größten für die Welle der Länge 


(40). 


Von allen geometrisch ähnlichen Wellen ist die mit der Länge /, am stärksten angefacht. 
Um die Wellenlänge der Welle mit der größten Anfachung zu bekommen, wobei für alle 
Wellen die gleiche Amplitude vorausgesetzt wird, schreiben wir (39) folgendermaßen: 


/ „oteh x h /o’ eteh (x h\—+ o\'» 3 
0,0035 7, | Ei ade, or) rn VER 
| o vg» \octgh(zh) - o, 


Hierin wird die linke Seite ein Maximum, wenn 6=0,--, tt, wird. Es wird also in diesem 


Fall die Welle der Länge 


- T Um 
4, = (42) 
49 


am stärksten angefacht. 


Für tiefes Wasser ergibt sich bei einer Windstärke von 104 em/see für 4, 7,7 em und für 
/, 17,3 em. Bei einer Windstärke von 250 em/seec ist 4, 44,5 cm und 4, = 100 em. 


15” 
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/um Schluß wollen wir noch die eefundenen Werte von n und s in (33) einführen. Es 
ereibt sieh dann für die Zeit des Anwachsens der Wellen 


R 
„oteh(zh) (u ,,— 0) = 
0.003575, (a,x) °—1 a 
| ovxı Oo 0.» A j 
f et 7 De N 
ni „.otehzh (u, — 0) m ” U 
VO, (ax) “—| 


Te Ö 


er a , er A n 
In Tabelle 3 sınd einige Zeiten anzereben, die die Wellen benötigen, um von ; 0015 auf 


dd 





- 0,08 anzuwachsen. 
/ 
Tabelle 3. 

Un: 250 em/ Sec 
/. 25 >0 1» 100 em 
/ 0.17 0.555 1.56 4.15 Std. 
n >) I00 500 S00 1000 em 
/ 0.465 2.17 7.09 20 66 Std. 


Zusammenfassung. Es ist die turbulente Strömung längs welligen Wänden untersucht 
worden und zwar wurde der Druckverlauf und die Geschwindigkeitsverteilung an vier Wänden 
nit verschiedenen Wellenformen bestimmt. Aus den Messungen wurde sowohl der Druck- 
als auch der Reibuneswiderstand der Wände bestimmt, weiter das Stromlinienbild und die 
Verteilung der Schubspannung. Ein Vergleich der Ergebnisse mit denen der turbulenten 
Strömung an glatten und rauhen Wänden ergab, daß in Wandnähe die Geschwindigkeits- 
verteilune dem logarithmischen Gesetz nicht mehr geborcht. In größerer Entfernung von der 
Wand ließen sie sich jedoch wiedergeben durch die Ausdrücke 


be Sale ı»y 
7 DD r,log Fr 
. Ad 
für Wellen mit scharfen Kämmen und 
+.) = 
a rt Pr „7 ® 
a =5,(5 v, log —- IS.2 te a, 
x r 


für Wellen mit runden Kämmen. 


Im zweiten Teil wurden die Ergebnisse der Druckwiderstandsmessungen auf die Frage 
nach dem Entstehen und dem Anwachsen von Wasserwellen bei Wind angewandt. Es ergab 
sich aus einer Energiebetrachtung, daß die ebene Wasseroberfläche für jede Windstärke 
instabil ist und daß es für jede Windstärke und jede Wellenlänge eine bestimmte maximale 
Amplitude gibt, «die nieht überschritten werden kann. Für einige Windstärken wurde diese 
Amplitude berechnet. Es zeigte sich, daß dieselben bei Windstärken unter 100 em/see sehr 
klein sind. Die stärkste Anfachung von allen geometrisch ähnlichen Wellen hat die mit der 
Länge 


und von Wellen mit der gleichen Amplitude die mit der Länge 


ü x Mu 
IR 
a 


Zum Schluß wurden noch die Anwachszeiten für einige Wellen verschiedener Wellenlänge 
bei zwei verschiedenen Windstärken berechnet. 647 
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Fin elektrisches Gerät zum Lösen von Systemen 
linearer Gleichungen. 
Von Helmuth Bode VDE, Berlin. 


(Mitteilung aus dem Elektrotechnischen Laboratorium der Technischen Hochschule Berlin.) 


n vielen Stellen der Technik besteht die Notwendigkeit, Systeme von linearen Gleichungen 
A mit einer größeren Anzahl von Unbekannten zu lösen. Es wurde ein Gerät gebaut, das 
die Abbildung und Lösung von symmetrischen, linearen Gleichungssystemen dureh ein elek- 
trisches Netzwerk in einfacher Weise und mit beliebig hoher Genauigkeit gestattet. 


1. Bekannte Verfahren zur Auflösung linearer Gleichungssysteme. Die Auflösung eines 
Systems linearer Gleichungen mit mehreren Unbekannten bereitet mathematisch keine Schwierig- 
keiten: sie erfordert jedoch, besonders wenn es sich um eine größere Anzahl von Gleichungen 
und Unbekannten handelt, eine umfangreiche und zeitraubende Rechenarbeit. Auch beı Be- 
nutzung der bekannten Näherungsmethoden') gestaltet sich die Auflösung noch außerordentlich 
umständlich; um die Näherungsverfahren anwenden zu können, sind häufig noch Umformungen 
der Gleichungen erforderlich. Es sind daher zahlreiche Versuche gemacht worden, lineare 
Gleiehungssysteme maschinell zu lösen. Der erste Versuch scheint von Lord Kelvin zu 
stammen: er baute ein mechanisches Gerät, bestehend aus Waagebalken mit Rollen und darüber- 
laufenden Sehnüren?). Die Zahl der Wagebalken und Schnüre ist gleich der Zahl der 
Gleichungen und Unbekannten. Eine Maschine von J. Nowak’) stellt ein zwangsläufiges 
Getriebe dar; Zahnradübersetzungen bilden die Gleichungskoeffizienten ah. Ganz ähnlich 
arbeitet ein Gerät von M. Nähbauer?). Ein hydraulischer Apparat wurde von K. Fuchs 
gebaut’). 

Neuere Geräte und Vorschläge beschreibt W. Cauer in „Elektrische Nachrichtenteehnik* 
1935, Bd. 12, Heft 5*); es handelt sich dabei um elektrisch arbeitende Apparate, wie die 
Mallocksche Maschine, die Abbildung eines Fachwerks durch eine elektrische Schaltung, 
wobei der Bedingung der Statik, daß die Summe aller auf einen Knotenpunkt wirkenden 
Kräfte gleich O0 sein muß, die elektrische Bedingung, daß die Summe aller einem Knoten- 
punkt zufließenden Ströme O0 sein muß, entspricht, u. a. 

Cauer selbst schlägt ein Gerät vor, mit dem sich nach einem Kompensationsverfahren 
die Elimination von Unbekannten elektrisch durehführen läßt: man erhält so ein gestaffeltes 
Gleichungssystem, dessen letzte Gleichung nur noch eine Unbekannte besitzt; damit ist das 
Gleichungssystem praktisch gelöst. Durch Verwendung von Schrittschaltwerken und Relais 
kann der Kompensationsvorgang weitgehend automatisiert werden. 

Eine Umgehung der Auflösung linearer Gleichungssysteme bedeuten sog. Netzmodelle, 
durch die zu berechnende elektrische Leitungsnetze nachgebildet werden, und an denen dann 
dureh Spannungs- und Strommessungen die erforderlichen Daten eines zu bauenden Netzes 
oder das Verhalten eines vorhandenen Netzes bei bestimmten Betriebsfällen ermittelt werden 
können. Es handelt sich dabei um die experimentelle Lösung der ursprünglichen Aufgabe 
bei geänderten Maßstäben. 


Ziel der vorliegenden Arbeit war es, die Möglichkeit zu untersuchen, beliebige lineare 
Gleichungssysteme, auch solche, die nieht aus einem Netzproblem entstanden sind, durch 
Abbildung in ein Netzwerk zu lösen, und eine praktische Ausführungsform für ein auf dieser 
Grundlage arbeitendes Gerät zu finden. 


2. Die Abbildung elektrischer Netzwerke durch lineare Gleichungssysteme. Die Berechnung 
eines Netzes von Ohmschen Widerständen führt auf ein System linearer Gleichungen mit 
reellen Koeffizienten; umgekehrt kann ein lineares Gleichungssystem unter gewissen Voraus- 
setzungen in ein Netz von Ohmscehen Widerständen abgebildet werden und durch Strom- 
und Spannungsmessungen gelöst werden. Die Aufstellung der Gleichungen eines Netzwerkes 
kann mit Hilfe der Kirchhoffschen Stromgleichungen oder Spannungsgleichungen erfolgen, 
wobei der Ausgangspunkt der Berechnung einmal die Knotenpunkte, das andere Mal die 
Maschen des Netzwerkes sind; man benutzt die Stromgleichungen, wenn die Spannungen in 
den Knotenpunkten gesucht sind und die den Knotenpunkten zufließenden Ströme bekannt 
sind. Die Spannungsgleiehungen werden verwendet, wenn die Maschenströme unbekannt sind 
und die an den einzelnen Maschen liegenden Spannungen gegeben sind. An zwei Beispielen 
soll der Unterschied zwischen diesen beiden Formen der Netzgleichungen gezeigt werden. 


‘) Bei Erscheinen der genannten Arbeit waren die hier beschriebenen Untersuchungen bereits abgeschlossen. 
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Zunächst seien für ein Netz mit 4 Knotenpunkten (Abb. 1), zwischen die die Leitwerte 
OO: TODE VOBEREnGee oeschaltet sind, und deren Spannungen gesucht sind, die Netzgleichungen auf- 


’ alt, u,)—+ 93 lu, — Uu,) + 90: lu, ") 
0, = (a1 (U, — U,) + Gas ' (U, — U,) + Geo‘ (Us — U,) NM 


b; Isı' (1, ,) T(Iss° (1, u,) + Jz0 ' (tt, U.) 


o = In (U — U) FT Ion’ (Un — Ur) F Goa‘ (Un — Us) 





Durch Umformunge erhält man hieraus: 


1 


7 2 tr Iıs TI) U, I. '%U Jız "U; Jo’ U, 
l, Ja Ur (Ist I2s t I20)° tt, Ja: U; U FIAT 
7 Sa U, Ja "U, r (Ysı F Ja» + 30) U; Io U 
t, Jo U, "Jo2 ' Ur Jos "Us + (Yoı + Joa + Joa) Up: 








Abb. 1. Abb. 2. 


In «diesem Gleiehungssystem ist eine Gleichung die Summe aller übrigen. Für die 
Ströme in den Zweigen des Netzes sind nämlich nieht die Spannungen in den Knotenpunkten 
selbst, sondern die Spannungsunterschiede zwischen den Knotenpunkten bestimmend. Die 
Spannung eines Knotenpunktes kann daher beliebig angenommen werden und eine Gleichung 
gestrichen werden. Man erhält dann ein System von 3 voneinander unabhängigen Gleichungen. 


nV 
', IYıs + ist 940) u, Jıa ' U, Sıs "Us | 
R .) 
l, Js, U, Tr Is WER + 920) "Us Jo; U; (=). 
l, : I; u, Ja "Mg + (Is: + MER) -7 0) "U; 


Entsprechend ergibt sich für das zu obigem Netzwerk duale Netz nach Abb. 2 mit den 
Widerständen 9 Fa Far... In den Maschen bei einer ähnlichen willkürlichen Annahme 


„0 


zur Eliminierung der überzähligen 4. Gleichung: 


"; UF u SP TI BEZ Yı2 't, Ya 1; | 

A . . >) 
", "nt, mit, + No; Fr Yo), Na; tz nn. A 
nes Na 7 Nana’ i; nu (r;, + "39 + No) b; | 


Die Gleichungssysteme (2) und (3) bestehen beide aus 3 voneinander unabhängigen 
Gleichungen; sie stimmen in ihrem Aufbau vollkommen miteinander überein, da sie zu zwei 
zueinander dualen Netzwerken gehören: man erhält das eine Netzwerk aus dem anderen, 
wenn man Widerstände durch Leitwerte, Spannungen durch Ströme, Netzmaschen durch Netz- 
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knoten und umgekehrt ersetzt. Die beiden Netze haben 4 Knotenpunkte bzw. 4 Maschen 
und führen auf Gleichungssysteme mit 3 voneinander unabhängigen Gleichungen; in gleicher 
Weise führt die Berechnung von Netzwerken mit » Knotenpunkten oder Maschen auf 
Gleichungssvsteme mit (n — 1) Gleichungen. | 


3. Die Abbildung linearer Gleichungssysteme durch elektrische Netzwerke. Um ein lineares 
Gleiehungssystem mit » Gleichungen elektrisch abzubilden, braucht man nach dem Vorher- 
oehenden ein Netzwerk mit (rn -+1) Knotenpunkten bzw. ein solehes mit (n—+1) Maschen; 
in ersten Fall sind die Koeffizienten der Gleichungen durch Leitwerte und die absoluten 
Glieder durch Ströme, die in die Knotenpunkte hineingeschickt werden, darzustellen; die 
Unbekannten werden durch Spannungsmessungen zwischen den Knotenpunkten erhalten. Im 
zweiten Fall werden die Koeffizienten durch Widerstände und die absoluten Glieder dureh 
Spannungen abgebildet; die Unbekannten erhält man durch Messung der Ströme in bestimmten 
/weizen des Netzwerkes. 

An dem Beispiel eines Gleichungssystems mit 3 linearen Gleichungen (4) werde nun 
eezeigt, unter welchen Bedingungen sich ein lineares Gleichungssystem durch ein Netzwerk 
aus Ohmsechen Widerständen oder Leitwerten abbilden läßt. 


N; d,'K4, Tr Aa Karl, KH | 


Y,=d,,'% — Ooo* Ko 7” UÜgg *' Ka . . . . . . . . . | N). 
x 
3 


21 l ‚ 22 p. 23 
Ns Ay, KT Ag d, + Ay 
Der Vergleich dieses Gleichungssystems mit den Systemen (2) und (3) ‚ergibt die Ströme 
und Leitwerte bzw. die Spannungen und Widerstände des abbildenden Netzes; man erhält 
im ersten Fall: 


t| N; ls ur Ur = N; 
d,, (Jia + Iıs tg) As —(l,, ir 
(yo FE (Js; + Jos ER Jon} (Ag; Ay, an: Js; 
d,, — (Isı + MER + zo) d;,; Zen zn Js» 


/wei Bedingungen müssen erfüllt sein, damit eine Abbildung durch ein Netzwerk Ohm.- 
scher Widerstände möglich ist: 


I. Da der Leitwert zwischen zwei Knotenpunkten unabhängig von der Stromriehtung 
ist, da also z. B. 9, = 9, Ist, muß auch für die Koeffizienten des abzubildenden Gleichungs- 
systems die Bedingung erfüllt sein, daß a;,;==a7;; ist; d. h. das abzubildende Gleichungs- 
system muß eine symmetrische Hauptdeterminante haben. 


2. Da zur Abbildung nur positive Leitwerte bzw. Widerstände zur Verfügung stehen, 


müssen die Koeffizienten «a,,, a,, und a,, negativ und die Koeffizienten «a, ,,«,, und a,, positiv 
sein; außerdem muß 


11? 


d;, 0, +0] 
(Ay, | Ay; + das] 


Ass d;, | Ay, | 


sein. Allgemein bestehen also die Bedingungen: 


Gl; dp <V a4;>|2axr 


Rust k=jzi 


Entsprechende Umrecehnungsformeln und die gleichen Bedingungen gelten auch für die 
Abbildung mit Widerständen und Maschen: außerdem besteht aber bei dieser Art der Ab- 
bildung noch eine weitere Einschränkung. Die Indizes der Widerstände geben an, zwischen 
welchen Maschen diese als Kopplungsglieder liegen; sind bei einem Gleichungssystem alle 
Koeffizienten von O0 verschieden, so müssen bei dem abbildenden Netzwerk alle Maschen 
unmittelbar miteinander durch Widerstände gekoppelt sein. Nur 4 Maschen jedoch können 
praktisch unmittelbar miteinander gekoppelt sein, d. h. einander direkt benachbart sein 
(vgl. Abb. 2); bei einer größeren Anzahl von Maschen ist nur eine mittelbare Kopplung aller 
Maschen untereinander möglich. Die Abbildung eines Gleichungssystems mit mehr als 
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3 Gleichungen nach dieser zweiten Art erfordert also die Verwendung von Hilfsmaschen und 
eine komplizierte Berechnungsart für die Widerstände des abbildenden Netzwerkes; sie kommt 
aus diesem Grunde für den Zweck der Gleichungslösung nicht in Frage. Sie würde den 
weiteren Nachteil haben, daß für »die Darstellung der absoluten Glieder der Gleichungen 
ebensoviel voneinander unabhängige Spannungsquellen nötig wären, wie das abzubildende 
(Gleichungssystem Gleichungen hat. 

Bei der ersten Art der Abbildung, bei der die Koeffizienten des Gleichungssystems durch 
Kopplungsleitwerte zwischen den Knotenpunkten dargestellt werden, besteht die vorstehende 
Einschränkung nicht, da es ohne Schwierigkeiten möglich ist, jeden Knotenpunkt mit jedem 
anderen durch einen Leitwert direkt zu koppeln; es können also Gleichungssysteme mit 
beliebige vielen Gleiehungen und Unbekannten durch Netzwerke mit einer um 1 größeren 
Zahl von Knotenpunkten abgebildet werden, wenn sie die beiden vorher genannten Bedingungen 
erfüllen. Die in das Netzwerk zu schiekenden Ströme können einer gemeinsamen Spannungs- 
quelle entnommen werden. Die Unbekannten werden in einfacher Weise durch Messung 
der Spannungen zwischen dem (n» +1).Knotenpunkt einerseits und den übrigen n Knoten- 
punkten andererseits erhalten. 


4. Erweiterung auf Koeffizienten beliebigen Vorzeichens. Die Zahl der Gleichungssysteme, 
die auf diese Art durch Abbildung ın ein Netz von Ohmschen Leitwerten gelöst werden 
können, ist besonders durch die zweite Bedingung eingeschränkt, während die erste Bedingung, 
die Symmetrie der Hauptdeterminante fordert, bei den meisten in der Technik vorkommenden 
Gleichungssystemen erfüllt ıst. 

Um jedes symmetrische Gleichungsysstem elektrisch abbilden und lösen zu können, 
müßte man für den Aufbau des Netzwerkes positive und negative Leitwerte zur Verfügung 
haben. Negative Leitwerte können durch Anordnungen von Elektronenröhren dargestellt 
werden; es ist sogar möglich, mit Hilfe von Röhren richtungsabhängige Leitwerte zu bauen, 
so daß auch die Bedingung der Symmetrie beseitigt werden könnte. Es ist jedoch sehr 
schwierig, derartige Leitwerte in einem großen Bereich veränderlieh und einstellbar zu machen; 
ihre zeitliche Konstanz und ihre Linearität in weiten Bereichen ist wahrscheinlich ungenügend; 
außerdem erfordern sie die Verwendung von zahlreichen besonderen Spannungsquellen, die 
den Bau und den Betrieb eines solchen Geräts zu kompliziert machen würden. 

Die vorstehenden Schwierigkeiten werden vermieden bei Betrieb des Gerätes mit Wechsel- 
strom und Benutzung von Blindleitwerten; als positive Blindleitwerte werden Kapazitäten, 
als negative Blindleitwerte Induktivitäten verwendet. Mit einem solchen Netzwerk kann 
jedes symmetrische, lineare Gleichungssystem mit reellen oder auch komplexen Koeffizienten 
aleichen Phasenwinkels abgebildet werden. Ein für den praktischen Gebrauch geeignetes 
Gerät muß eine größere Anzahl kapazitiver (positiver) und induktiver (negativer) Leitwerte 
besitzen, «die einem zu lösenden Gleiehungssystem entsprechend zu einem Netz zusammen- 
geschaltet werden können und in weiten Grenzen veränderlich sind. Die positiven und nega- 
tiven Leitwerte sollen rein imaginär sein, damit sie die Voraussetzung ihrer Verwendbarkeit, 
entzegengesetztes Vorzeichen zu besitzen, erfüllen. Die in das abbildende Netzwerk hinein- 
zuschiekenden Ströme müssen gleichphasig sein, damit die von ihnen an den einzelnen Leit- 
werten hervorgerufenen Spannungen sich als Beträge addieren oder subtrahieren. Die Frequenz 
der zur Speisung des Netzwerkes verwendeten Spannungsquelle muß konstant sein, damit 
positive und negative Leitwerte eine feste Eiechung in Leitwerteinheiten erhalten können. 
Die Messung der Spannungen im Netzwerk zur Ermittelung der Unbekannten muß mit ver- 
brauchslosen Spannungsmessern erfolgen, damit durch ihre Einschaltung keine Änderung der 
Leitwerte des Netzes verursacht wird. 


5. Die Ausführung der abbildenden Blindleitwerte. Als kapazitive Leitwerte sind Dreh- 
kondensatoren mit Luftdielektrikum besonders gut geeignet; sie sind in weiten Grenzen in ihrem 
leitwert veränderlich und erhalten bei halbkreisförmigen Platten eine glejehmäßig geteilte 
Skala; sie erfüllen die Forderung, einen rein imaginären Leitwert zu besitzen, in sehr hohem 
Maße. 

Als günstieste Lösung für die Darstellung eines veränderlichen, induktiven Leitwertes 
wurde die Parallelsehaltung einer festen Induktivität, die einen Leitwert von der Größe des 
Maxıimalwertes der verwendeten Drehkondensatoren hat, mit einem dieser veränderlichen, 
kapazitiven Leitwerte gefunden. Durch Änderung der Kapazität des Drehkondensators von 
etwa 0 bis zum Höchstwert läßt sich der induktive Leitwert der Parallelschaltung von seinem 
Maximalwert bis 0 verändern. Die Drehkondensatoren erhalten eine Eichung in positiven 
leitwerteinheiten, die dann gilt, wenn sie allein benutzt werden, und eine zweite in nega- 
tiven Leitwerten für den Fall, daß sie in Parallelschaltung mit einer festen Induktivität zur 
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Darstellung eines veränderlichen, induktiven Leitwertes verwendet werden. Auch im zweiten 
Fall ist der Leitwert linear von der Kapazität abhängig, die Eichkurve also eine gerade Linie: 


| 

A: * -(] "+ L- 0). 
Q J o ? 4 
Dabei ist zunächst anzenommen, daß auch die benutzten Spulen als reine Blindwiderstände 
anzusehen sind; diese Voraussetzung ist nur angenähert erfüllt; damit das Gerät schnell 
cenaue Ergebnisse liefert, ist jedoch die Verwendung weitgehend rein imaginärer Leitwerte 
unbedingt erforderlich. Es ist daher noch zu prüfen, wie sich die sog. Güte y Mr der 

. i 
Parallelschaltung mit der parallelgeschalteten Kapazität ändert. 
» sei die Güte der Spulen ohne Parallelkapazität bei der später zu benutzenden Meß- 
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frequenz, y; die Güte der Parallelschaltung, q, und g, die entsprechenden Leitwerte. 
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Da die Güte der Spulen selbst sehr hoch gewählt wird, ist: 
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Aus dem letzten Ausdruck für q, ist zu ersehen, daß der imaginäre Leitwert der Parallel- 
schaltung ungefähr linear mit der Kapazität abnimmt wie bei Benutzung einer reinen Induk- 
tivität, daß zu dem imaginären Anteil des Leitwertes aber ein ungefähr konstant bleibender 
reeller Anteil hinzukommt. Für die Güte der Parallelschaltung erhält man folgenden Ausdruck: 

2 y,.(1 o°- L-C), 
d. h. die Güte nimmt mit dem Leitwert der Parallelschaltung linear ab. Soll daher der 
induktive Leitwert im Verhältnis 1:100 geändert werden können und dabei die Güte 5 nicht 
unterschritten werden, so müßte die Spule selbst eine Eigengüte 500 haben; dieser Wert wird 
in der Hochfrequenztechnik ungefähr erreicht. 

Es standen Drehkondensatoren mit einer Maximalkapazität von ca. 4500 pF in größerer 
Anzahl zur Verfügune. Für die induktiven Leitwerte wurden aus Gründen der einfachen : 
Herstellung und des geringen Raumbedarfs Induktivitäten von 10 mHl gewählt. Nimmt man 
den Bereich, in dem die Leitwerte veränderlich sein sollen, mit etwa 1:100 an mit den 
vorhandenen Kondensatoren wurde nur ein Bereich 1:67,5 erzielt —, so erhält man einen 
Kleinstwert der Kapazitäten von 45 pF; da positive und negative Leitwerte natürlich den 
gleichen Größenbereich haben sollen, ergibt sich die Frequenz der zu verwendenden Spannungs- 
quelle von selbst: 
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Man erhält in diesem Falle eine Frequenz von 239 kHz (» = 1,5: 10% s '). Die verwendeten 
leitwerte hatten bei dieser Frequenz einen Bereich (gegeben durch die Anfangs- und End- 
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kapazität der Drehkondensatoren) von 0,1: 10° bis 6,75 - 10 (> ; als Leitwerteinheit wurde 


ii 
10-3. 7 gewählt. 


Als Spannungsqueile kommt bei der erforderlichen, hohen Frequenz von f-- 239 kHz nur 
ein Röhrengenerator in Frage. Es wurde zunächst ein selbsterregter Röhrensender mit einer 
Telefunkenröhre RS 213 benutzt, dem die Hochfrequenzspannung mit einer verschiebbaren 
Ankopplungsspule entnommen wurde; sie konnte durch Änderung der Kopplung dieser Spule 
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mit dem Sender und durch Ändern der Heiz- und Anodenspannung des Senders geregelt 
werden. Da die Frequenz dieses Senders infolge des Auftretens von Kopplungswellen bei 
Resonanznähe in irgendeinem Teil des Netzwerkes nicht die erforderliche Konstanz aufwies, 
wurde später ein Sender mit Fremdsteuerung benutzt (s. Abb. 3). Eine praktisch absolute 
requenzkonstanz könnte durch Quarzsteuerung erzielt werden; jedoch genügte die Fremd- 
steuerung des Senders schon vollkommen. Außerdem wurde dureh die Fremdsteuerung eine 
sehr große Oberwellenfreiheit erreicht, da der Sender jetzt 3 auf die Grundwelle abgestimmte 
Kreise besaß. 
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Abb. 3. 


Das eine Ende der Ankopplungsspule wird fest mit dem Nullpunkt des Netzwerkes, 
d.ı. der (n-+-1). Knotenpunkt, verbunden, während das andere Ende über je einen Strommesser 
und Regelwiderstand mit allen übrigen Knotenpunkten verbunden wird. Die einzelnen Ströme 
sollen gleiehphasig sein; sie sind das schon insoweit, als sie nur gleichgerichtet oder entgegen- 
gerichtet sein können, da das Netzwerk ja nur positiv oder negativ imaginäre Leitwerte, 
nicht jedoch komplexe Leitwerte besitzt. Als Regelwiderstände kommen daher auch nur 
imaginäre Widerstände in Betracht; besonders gut eignen sich hier wieder Drehkondensatoren. 


6. Messung der Spannungen und Ströme. Als Strommesser wurden zunächst Hitzdraht- 
strommesser benutzt: da der ausnutzbare Meßbereich dieser Instrumente sehr klein ist, und 
da sie außerdem durch ihren Ohmschen Widerstand die Gleichphasigkeit der Ströme im 
Netzwerk beeinträchtigen, wurden später die Ströme durch Messung der Spannung, die sie an kapa- 
zitiven Widerständen hervorriefen, gemessen; es wurden also in die Zuleitung zu jedem 
Knotenpunkt mehrere Kapazitäten verschiedener Größenordnung eingeschaltet, an denen mit 
einem Röhrenvoltmeter die Spannung gemessen werden konnte; die Röhrenvoltmeter erhielten 
Kiehungen in Stromeinheiten für die einzelnen Kapazitäten. Diese Strommessung hat den 
Vorteil, daß der Meßbereich sehr groß ist, daß die Gleichphasigkeit der Ströme nicht gestört 
wird, und daß eine Überschreitung des Meßbereiches die Instrumente nicht gefährdet. 

Als Röhrenvoltmeter für die Strommessung wie auch für die eigentlichen Spannungs- 
messungen im Netzwerk zur Bestimmung der Unbekannten wurden Röhrenvoltmeter mit 
Gittergleichriehtung verwendet, da sich alle anderen Arten von Röhrenvoltmetern, die mit 
negativer Gittervorspannung arbeiten, im vorliegenden Falle als ungeeignet erwiesen. Da im 
Netzwerk wie bei der Strommessung Spannungen an rein kapazitiven Widerständen zu messen 
sind, würde eine negative Gittervorspannung statische Aufladungen des Netzes zur Folge 
haben: der Ausschlag des Instruments würde nicht nur von der zu messenden Wechsel: 
spannung, sondern auch von der statischen Spannungsverteilung im Netz abhängen. 

Das Röhrenvoltmeter mit Gittergleichriehtung, das ohne eine äußere negative Vor- 
spannung arbeitet, besitzt diesen Nachteil nicht; seine Schaltung zeigt Abb. 4; der Meßbereich 
dieses Instruments reichte von O,1 bis etwa 5 Volt; da aber im Netz Spannungen bis etwa 
O0 Volt auftraten, genügte das Instrument in dieser Form noch nicht; durch Umkehrung 
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Abb. 4 und 5. Abb. 6. 
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der Funktionen der Elektroden einer Röhre kann man jedoch ein Röhrenvoltmeter erhalten, 
das je nach der verwendeten Röhre einen Meßbereich bis 150 Volt und mehr besitzt. Die 
zu messende Spannung wird nieht an das Gitter und die Kathode der Röhre, sondern an 
Anode und Kathode gelegt, während im Gitterkreis eine Batterie und ein Strommesser liegen. 
Die Erhöhung des Meßbereiches bei dem umgekehrten Röhrenvoltmeter hängt von dem Durch- 
ori! D der verwendeten Röhre ab; sie beträgt beinahe das 1/D°’-fache; Abb. 5 zeigt die 
Schaltung des umgekehrten Röhrenvoltmeters. 


7. Bestimmung des Vorzeichens der Unbekannten. Zur Ermittlung der Unbekannten eines 
abeebildeten Gleichungssystems ist es nicht nur erforderlich, die Größe der Spannungen 
zwischen den Knotenpunkten des Netzes und dem Nullpunkt zu messen, sondern diese Spannungen 
uch dem Vorzeichen nach zu bestimmen. Diese Vorzeichenbestimmung kann nur indirekt 
erfoleen. Jede unbekannte Spannung wird zunächst der Größe nach gemessen; dann wird 
lie Differenz zwischen dieser Spannung und einer dem Vorzeichen nach bekannten Spannung 
oemessen; ist diese Differenz kleiner als die zuerst gemessene Spannung, so hat die letztere 
das gleiche Vorzeichen wie die Hilfsspannung; Ist die Differenz größer, so hat die zu bestimmende 
Spannung das entgegengesetzte Vorzeichen wie die Hıilfsspannung. Die praktische Ausführung 
dieser Vorzeichenbestimmung erfolgte mit einer Meßanordnung nach Abb. 6. 

Zur Bestimmung der Unbekannten #, (Gleichungssystem 4) soll in dem abbildenden 
Netzwerk die Spannung «, zwischen dem Knotenpunkt I und dem Nullpunkt der Größe und 
dem Vorzeichen nach ermittelt werden. €,, ist der Regelkondensator, mit dem der Strom, 
der in den Knotenpunkt 1 zu schicken ist, auf die erforderliche Größe eingestellt wird; an 
die hochfrequente Spannungsquelle ist außer dem eigentlichen Netz noch eine Reihenschaltung 
einer festen Kapazität e, und einer veränderlichen Kapazität e, angeschlossen; an der Kapa- 
zität e, liegt also unter allen Umständen eine der Gesamtspannung gleichgerichtete Spannung 
und, da das gesamte Netz nur imaginäre Leitwerte besitzt und in ihm daher nur gleich- 
verichtete oder entgegengerichtete Spannungen auftreten können, eine positiv zu nennende 
Spannung. Die Kapazität e, ist zunächst bei der Messung der Größe der Spannung «u, kurz- 
eeschlossen; die Kathode der Meßröhre liegt am Nullpunkt des Netzes. Zur Vorzeichen- 
bestimmung wird die Taste T gedrückt und damit der Kurzschluß des Kondensators e, auf- 
gehoben; die Kathode der Meßröhre erhält dadurch gegenüber dem Nullpunkt des Netzes eine 
positive Spannung. Ist die dem Vorzeichen nach zu bestimmende Spannung «, nun auch 
positiv, so liegt an der Röhre jetzt eine geringere Spannung, und das Meßinstrument zeigt 
eine geringere Spannung an; ist «, negativ, so wird von dem Röhrenvoltmeter eine größere 
Spannung gemessen. Die Drosseln Z von je 0,1 MH stellen für die Hochfrequenz praktisch 
eine Isolierung dar und machen den Einfluß von Kapazität der Batterie und deren Zuleitungen 
vegen Erde auf den Meßvorgang unwirksam; das gesamte Röhrenvoltmeter soll möglichst 
weitgehend von allen anderen Teilen des Gerätes für Hochfrequenz isoliert sein. Bedingung 
für die Richtigkeit des Ergebnisses der Vorzeichenbestimmung ist, daß die Hilfsspannung 
an der Kapazität e, kleiner ist als das Zweifache der zu messenden Spannung «,; durch 
Anderung der Kapazität e, kann diese Bedingung auch in kritischen Fällen immer erfüllt werden. 

Zur Darstellung der absoluten Glieder eines Gleichungssystems sollen in die Knoten- 
punkte des abbildenden Netzwerkes Ströme von der Größe dieser Glieder geschickt werden; 
die gleichzeitige Speisung des Netzes mit allen diesen so eingestellten Strömen bereitet 
Schwierigkeiten, da bei der Einstellung eines Stromes die bereits eingestellten Ströme sich 
mehr oder weniger ändern; eine Verringerung dieser Rückwirkung ist durch Wahl sehr kleiner 
Itegelkondensatoren möglich, erfordert jedoch die Verwendung einer erheblich höheren, hoch- 
[requenten Spannung. 


8. Das Superpositionsverfahren. Hier hilft die Anwendung der Methode der Superposition: 
Die Gesamtlösungen eines Gleichungssystems mit » Gleichungen sind die Summen von Teil- 
lösungen, die von den absoluten Gliedern der einzelnen Gleichungen herrühren. Setzt man 
in einem Gleichungssystem alle absoluten Glieder bis auf das einer Gleichung gleich O0 und 
löst dieses vereinfachte Gleichungssystem auf, so erhält man für die einzelnen Unbekannten 
die Teillösungen, die von dem berücksichtigten, absoluten Glied herrühren; löst man in dieser 
Art die » verschiedenen Gleichungssysteme, die entstehen, wenn immer ein anderes absolutes 
Glied berücksichtigt wird, während alle anderen gleich 0 gesetzt werden, und bildet die 
Summen der Einzellösungen für die verschiedenen Unbekannten, so erhält man die Lösungen 
des vollständigen Gleichungssystems. Bei der elektrischen Auflösung von linearen Gleichungs- 
systemen bedeutet das Folgendes: In das abbildende Netzwerk werden nieht alle Ströme 
gleichzeitig hineingeschiekt und durch Messung der Spannungen zwischen dem Nullpunkt des 
Netzes und den » Knotenpunkten die Gesamtlösungen des aufzulösenden Systems erhalten, 
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sondern nacheinander werden die einzelnen Ströme in die entsprechenden Knotenpunkte 
geschickt, und jedesmal werden durch Spannungsmessungen Teillösungen ermittelt. Die Summen 
der » gemessenen Spannungen der gleichen Knotenpunkte ergeben die Lösungen des gegebenen 
Gleichungssystems. Die Anwendung der Überlagerung bringt außer der einfachen Einstellung 
der Ströme noch weitere Vorteile: Es wird nur ein Strommesser gebraucht, d. h. ein Röhren- 
voltmeter für die Strommessung; die Bedingung der Gleichphasigkeit der Ströme als Voraus- 
setzung für die Richtigkeit der Ergebnisse ist von selbst erfüllt, da statt der elektrischen 
Summierung eine rechnerische Addition der Beträge der Teillösungen erfolgt; zur Speisung 
des Netzwerkes wird eine weniger große Stromquelle benötigt. Gegenüber diesen Vorteilen 
fällt der Nachteil der größeren Zahl der Messungen meist nicht ins Gewicht. 
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Abb. 7. 


In Abb. 7 ist die Prinzipschaltung des Gerätes für Gleiehungslösung dargestellt, das 
nach den im Vorhergehenden erläuterten Gesichtspunkten gebaut wurde. 

RA ist die Strommeßeinriehtung, die aus den Kapazitäten X, und A, und einem Röhren- 
voltmeter besteht: A, hat die Größe 1/9: K, und ist bei Messung größerer Ströme durch einen 
Schalter kurzzeschlossen; durch Öffnung des Schalters S kann der Meßbereich des Röhren- 
strommessers auf '/,, des ersten Bereiches verkleinert werden. €, ist ein Drehkondensator 
zur Einstellung der Ströme auf die vorgeschriebene Größe; mit dem Wahlschalter W kann 
der Strommesser und der Regelkondensator mit den einzelnen Knotenpunkten des Netzes 
verbunden werden. Mit dem Röhrenvoltmeter RV werden jedesmal die Spannungen der 
Knotenpunkte 1,2,3,........ n gegen den Nullpunkt 0 des Netzes bei der Speisung des Netz- 
werkes in den einzelnen Knotenpunkten gemessen; die Reihenschaltung des Festkondensators e, 
und des Drehkondensators e, dient zur Vorzeichenbestimmung (s. vorher). 


9. Die Genauigkeit der gefundenen Lösungen. Nachdem ein praktisch brauchbarer Aufbau 
für das Gerät gefunden war, mußte noch untersucht werden, mit welcher Genauigkeit Lösungen 
von Gleiehungssystemen mit ihm zu erhalten sind. Eine Angabe der Ungenauigkeit der 
Ergebnisse auf Grund bekannter Fehler des Gerätes, z. B. infolge von Wirkanteilen der Leit- 
werte ist nicht möglich, da der Fehlwinkel der Leitwerte nicht konstant ist, und da die Fehler 
der Ergebnisse außer von der ungenauen Abbildung eines Gleichungssystems durch ein Netzwerk 
auch von dem Betriebszustand dieses Netzwerks, z. B. davon, in welchem Punkte gerade 
die Speisung stattfindet, abhängt. Trotz der ungenauen Abbildung und trotz fehlerhafter 
Spannungsmessungen lassen sich jedoch mit dem Gerät beliebig genaue Lösungen ermitteln, 
wenn man das Lösungsverfahren wiederholt anwendet, eine Möglichkeit, auf die schon Lord 
Kelvin in der Beschreibung seiner Lösungsmaschine hingewiesen hat. 

Setzt man die zuerst erhaltenen Näherungslösungen für die Unbekannten in die zu 
lösenden Gleichungen ein, so befriedigen sie diese nicht vollkommen: die gefyndenen Diffe- 
renzen können als die absoluten Glieder eines neuen Gleichungssystems angesehen werden, 
das sich im übrigen von dem ursprünglichen System nicht unterscheidet. Die Abbildung 
durch ein Netz ist also die gleiche wie die des ersten Systems; nur die Speisung des Netzes 
hat mit anderen Strömen zu erfolgen. Die durch Spannungsmessungen erhaltenen Lösungen 
des neuen Gleichungssystems stellen die ersten Verbesserungen der zuerst gefundenen Näherungs- 
lösungen des ursprünglichen Systems dar. Man erhält also verbesserte Lösungen für die 
Unbekannten durch Addition der ersten Näherungslösungen und der ersten Verbesserungen; 
setzt man diese verbesserten Lösungen in die zu lösenden Gleichungen ein, so werden diese 
zwar besser, aber noch nicht vollkommen erfüllt; die jetzt gefundenen Differenzen bilden 
wiederum die absoluten Glieder eines neuen Gleichungssystems, dessen Lösung neue Ver- 














7, Heft 4 be si r nn 
Band 11. His Bode, Ein elektrisches Gerät zum Lösen von Systemen linearer Gleichungen 221 





hesserungen ergibt; durch Wiederholung dieses Verfahrens kann die Genauigkeit der Lösungen 
für die Unbekannten des aufzulösenden Gleiehungssystems beliebig weit getrieben werden. 
Eine Kontrolle für die jeweils in den einzelnen Stufen erreichte Genauigkeit ermöglicht die 
in jeder Stufe erfolgende Bestimmung neuer absoluter Glieder; aus dem Größenverhältnis 
der neuen absoluten Glieder zu den ursprünglichen absoluten Gliedern ersieht man am ein- 
fachsten die wahrscheinlichen Fehler der absolut größten Unbekannten. 

Die vorstehend beschriebene Art der Bestimmung von Verbesserungen für die Lösungen 
setzt nicht voraus, daß die Abbildung des zu lösenden Systems vollkommen fehlerlos erfolgt 
ist, oder daß die das Netz speisenden Ströme genau eingestellt werden; sie erlaubt vielmehr 
die Bestimmung beliebig genauer Lösungen in jedem Falle, bei ungenauer Abbildung der 
Koeffizienten dureh elektrische Leitwerte, bei Verwendung nicht rein imaginärer Leitwerte, 
bei ungenauer Einstellung der das Netz speisenden Ströme und bei ungenauer Messung der 
Spannungen zwischen den Knotenpunkten. Der Grund für diese Möglichkeit der Bestimmung 
beliebig genauer Lösungen ist der, daß bei jeder Verbesserung auf die ursprünglichen 
Gleichungen und nicht auf die abgebildeten Gleichungen zurückgegangen wird, und so dauernd 
die Fehler, die durch ungenaue Abbildung und Messung entstehen, durch rechnerische Zwischen- 
kontrolle ausgeglichen werden. Die schließlich erhaltenen Endergebnisse brauchen durchaus 
nieht mehr eine auch nur angenäherte Lösung des Abbildungssystems darzustellen; sie sind 
jedoch eine Näherungslösung für das aufzulösende System. 


10. Das Einheitsstromverfahren. Die Bestimmung beliebig genauer Lösungen für ein 
Gleichungssystem bei nur einmaliger Benutzung des elektrischen Lösungsgeräts ermöglicht 
ein Verfahren, das als Einheitsstromverfahren zu bezeichnen wäre. Bei Anwendung der 
Superposition, d. h. bei Speisung des abbildenden Netzes jeweils nur in einem Knotenpunkt 
und Addition der 'Teillösungen für die Speisung in den verschiedenen Punkten sind diese 
Teillösungen den speisenden Strom proportional; es ist daher nicht erforderlich, in einen 
Knotenpunkt den Strom zu schicken, den das absolute Glied der entsprechenden Gleichung 
vorschreibt, sondern man kann einen beliebigen Strom zur Speisung benutzen und die erhaltenen 
Spannungen mit einem dem absoluten Glied entsprechenden Faktor multiplizieren; zweck- 
mäßigerweise wird man diesen Strom so auswählen, daß nur mit dem Zahlenwert des absoluten 
Gliedes zu multiplizieren ist. Die Gleichungen, die zur Ermittlung von Verbesserungen für 
die ersten Näherungslösungen zu lösen sind, unterscheiden sich von den ursprünglichen 
Gleichungen nur durch die absoluten Glieder; statt nun das abbildende Netz mit Strömen, 
‘lie gleich den neuen absoluten Gliedern sind, nacheinander zu speisen und die Teillösungen 
zu bestimmen, kann man die letzteren einfacher durch Multiplikation der Spannungen, die 
schon früher bei Speisung mit der Stromeinheit gemessen wurden und zur Ermittlung der 
ersten Näherungen dienten, mit neuen Faktoren erhalten. In gleicher Weise können alle 
weiteren Verbesserungen aus den Spannungsmessungen bei Speisung mit der Stromeinheit 
bestimmt werden. 


11. Allgemeiner Lösungsgang. Der Lösungsgang eines Gleichungssystems (5) mit n 
Gleichungen gestaltet sich jetzt folgendermaßen: 


ya, '%, +0,’ % +4, '’%t:.....t7 9°’ En 
ur Ad, 8, 7 Ayg da + Aa '%Tt.-.-- + 4G,n Con 
} B | 1»). 
N; Az, 7 U FFErEEn ET Zee ul U Free 2 Passau rAzn Ku 
Ya Ay Ct las Kt las ten. + Ayn Con 





Bedingung für die Möglichkeit der elektrischen Auflösung überhaupt ist die Symmetrie 
der Hauptdeterminante, d. h. es muß a;;, = a;; Sein. 

Zur Abbildung des Systems mit » Gleichungen ist ein Netzwerk mit (» +1) Knoten- 
punkten erforderlich; die Größe der Kopplungsleitwerte zwischen den Knotenpunkten (die 
Indizes geben an, welche beiden Knotenpunkte von dem betreffenden Leitwert verbunden 
werden) erhält man aus den Koeffizienten des Gleichungssystems nach folgenden einfachen 
Formeln: 


YA, + (As +0, Piss fi An Je 4 
Iso — Qss + ds; + Ay; ee + Agy Js3 Ayg 
Io — Ag r d;; + (Ay rar + Ay Jıs Az 


Ino Aynt Anı eu Angst 0. + Ayumn-v Imn” Imn- 
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leitwerte, für die sich positive Werte ergeben, werden durch Einschaltung von Dreh- 
kondensatoren allein zwischen den entsprechenden Knotenpunkten dargestellt, Leitwerte, für 
die negative Werte erhalten werden, durch Drehkondensatoren in Parallelschaltung mit festen 
Induktivitäten:; die Einstellung der Leitwerte auf die geforderten Werte erfolgt in beiden 
Fällen nur mit den Drehkondensatoren, die zu diesem Zweck eine Eichung in positiven und 
negativen Leitwerten besitzen. 

Nun wird das abbildende Netz nacheinander in den Knotenpunkten 1 bis » mit einem 
Strom von der Größe der Stromeinheit gespeist; bei der Speisung im Knotenpunkt 1 werden 
zwischen dem Nullpunkt des Netzes einerseits und den Knotenpunkten I bis » andererseits 
die Spannungen 


U, u u 


11 « 


129 13° 7 RE | m ZI 


nach Größe und Vorzeichen gemessen; bei Speisung im Knotenpunkt 2 werden die Spannungen 


I, Wa: Mare er Usn 


219 22) 23 
gemessen, bei Speisung im Punkte » die Spannungen 


Hynıs Une; Uns Unas a 2 Hun: 


Die Teillösungen für das Gleiehungssystem, bei dem alle absoluten Glieder bis auf eines 
oleich 0 gesetzt sind, erhält man dureh Multiplikation der entsprechenden, bei der Stromeinheit 
gemessenen Spannungen mit diesem absoluten Glied; die Teillösungen für den Fall, daß das 
Gleichungssystem nur das absolute Glied „, besitzt, sind also die folgenden: 


d“ Y U,» Jia Yı' Has Kg NY Us, 
Kun Y,' Ua; a er ee Cın W* U,n- 


Bei alleinigem Vorhandensein des absoluten Gliedes 7, erhält man die Teillösungen: 


| Ur s Jar ’ ya - N: z Uys ‚ Uyz UR P Us; ’ 
ET Vans een eh in 


Die ersten Näherungslösungen für das gesamte Gleichungssystem sind dann die folgenden: 


. » . ! s . Bern | A F x 
er ra FE Pr Pas ya U Pr Fr er F Ya’ Uns 
a. pi Fa Pin anns Pas ut u FU a nn + Ya’ Uns 
® : R ’ | | 5 I | Be. 1m | i 
un trin t un t: Hr I Lan = Yı Yin FY’ Yan ty Want + +++. + Yı Un: 


Die ersten Näherungslösungen werden in die zu lösenden Gleichungen eingesetzt; sie erfüllen 
diese nieht vollkommen, sondern es ergeben sich die Differenzen 


ea #n: 


Diese Differenzen sind die absoluten Glieder eines neuen Gleichungssystems, durch dessen 


l,ösune man die ersten Verbesserungen der ersten Näherungslösungen erhält: die ersten Ver- 
besserunzen haben folgende Werte: 


’ ’ ’ | ’ | ' 
Le, Y, u, Yo’ Us TYs'’ Ust: Hr, 11. Yn’ Uns 
‚ ‚ ‚ D 
1 Zi 1 u, au „ 2 U, „ 3 Use ZEIT EEG 5 „ N Us % 
R ‚ ' ’ ’ N | r 
vn YıdıınT Ya’'lon TYs’UsnT-:::+: +... TYn' Uun- 


Die verbesserten Lösungen des Gleiehungssystems sind die Summen der ersten Näherungs- 
lösungen und der ersten Verbesserungen; werden diese in die zu lösenden Gleichungen ein- 
gesetzt, so werden neue Differenzen erhalten, aus denen wie oben neue Verbesserungen 
bestimmt werden können; durch wiederholte Anwendung dieses Verfahrens kann eine beliebige 
Genauigkeit der Lösungen erreicht werden. 
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Maßgebend für die Schnelligkeit, mit der die gewünschte Genauigkeit erreicht wird, 
ist natürlich die Genauigkeit der bestimmten Teillösungen und damit der ersten Näherungs- 
lösungen; mit dem gebauten Versuchsgerät wurde bei den ersten Näherungslösungen eine 
Genauiekeit von etwa 2,5°/, für die größte Unbekannte erzielt. 


12. Anwendungsbeispiel. Es sei das folgende Gleichungssystem mit 3 Gleichungen und 
3 Unbekannten zu lösen: 


104 61x, + 19 x, +0, | 
31 190, 13x, Me, (6). 
>| Va, — 21 27 3) X, 


Für die Leitwerte des abbildenden Netzes erhält man folgende Werte: 


Jo 12, (so vo 19, J30 + IH, 
Jı2 19, Jıs 30, (s=+t21. 


Die Leitwerte Yo; Foo: 91, und 9,, Sind hiernach durch Parallelschaltungen von Induktivitäten 
und Drehkondensatoren und die Leitwerte 9g,, und 9,, durch Drehkondensatoren allein ein- 
zustellen. Man kann jedoch auch für die negativ berechneten Leitwerte Drehkondensatoren 
und für die positiv berechneten Parallelschaltungen wählen; im vorliegenden Fall wurden 
z. B. die Leitwerte q,, und 4,, durch Parallelschaltungen von Induktivitäten und Drehkonden- 
satoren und alle übrigen Leitwerte durch Drehkondensatoren allein eingestellt. 


Die rechnerische Lösung ergibt für die Unbekannten folgende Werte: 
ge 1.6228 A 2048 x VI0O33.... 


Das Lösungesgeerät lieferte unter Anwendung des Einheitsstromverfahrens die ersten 
Näherungen 


2, 1.855, x,—= —940, P3 0,02, 
die zweiten Näherungen 
1.204, 3169, VON 
und die dritten Näherungen 
1.6233 . 12032, VOII. 


Die vorstehenden Ergebnisse zeigen die Brauchbarkeit des gebauten Gerätes für die 
lösung symmetrischer, linearer Gleichungssysteme. Die Vorteile dieses Verfahrens gegenüber 
der reehnerischen Lösung treten besonders bei Gleichungssystemen mit einer größeren Anzahl 
von Unbekannten und Gleichungen hervor. 


Für die Anregung zu der vorliegenden Arbeit und für viele wertvolle Ratschläge bei 
ihrer Durchführung danke ich Herrn Professor Schwenkhagen, jetzt an der Technischen 
Hochschule in Danzie. 685 
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Über den Träger auf elastischer Unterlage. 
Von K. Marquwerre in Berlin-Adlershof. 


1. Einleitung. Der Träger auf elastischer (nachgiebiger) Unterlage ist wegen seiner 
eroßen technischen Bedeutung schon seit langer Zeit der Gegenstand wissenschaftlicher Unter- 
suehunzen gewesen. Kine besonders eingehende Behandlung hat das Problem in dem Buche 
eleichen Titels von K. Havashi gefunden'). Die Voraussetzung der Theorie ist die von 
Zimmermann?) benutzte Annahme der Proportionalität von Normaldruck und Normaldurch- 
senkung der nachgiebigen („bettenden“*) Unterlage; diese Annahme ist streng genommen nur 
erfüllt für eine Unterlage, die (wie die Flüssigkeiten) keinerlei Schubspannungen zu übertragen 
vermag man könnte also vom schwimmenden Träger sprechen. 


Die vorliegende Untersuchung betrachtet den Träger auf einer „wirklich“ (d.h. im 
HHookeschen Sinne) elastischen Unterlage. Sie erörtert insbesondere auch die Frage, ob ein 
auf diese Weise gelagerter, durch konzentrierte Lasten beanspruchter Träger sieh in einiger 
Entfernung vom Lastangriffspunkt von der Unterlage abhebt, wie dies aus der Theorie des 
schwimmenden Trägers folgt. 


Um den Unterschied (und auch die Verwandtschaft) der beiden Hypothesen über das 
Verhalten des Baugrundes deutlich hervortreten zu lassen, sollen beide Fälle in gleicher 
Weise nebeneinanderher behandelt werden. Als Beispiel diene der einfachste Fall: der 
unendlich lange Balken unter Einzellast, gebettet auf dem „elastischen“ bzw. „nachgiebigen“ 
Halbraum. Um von diesem ebenen Problem (was hier nur in den Resultaten geschehen soll) 
leicht auf das analoge achsensymmetrische Problem (die allseits unendlich ausgedehnte Platte 
unter Einzellast) übergehen zu können, soll aber nicht der Balken (ebener Spannungszustand), 
sondern die unendliche Platte unter Schneidenlast senkrecht zur betrachteten Ebene durch- 
eereehnet werden (ebener Formänderungszustand). Diese Aufgabe unterscheidet sich im Er- 
vebnis von der eigentlichen Balkenaufgabe nur in Größen, deren Zahlenwert (in engen Grenzen) 


u | 
von der Wahl der Querkontraktionsziffer m beeinflußt wird. 
* 


2. Koordinaten und Bezeichnungen. Die elastische Platte liege waagerecht auf der Unter- 
lage; die „Achse falle in die Berührungsebene, senkrecht zu der Schneidenlast (y-Achse); 
die Schneidenlast P’ (als verteilte Last p=pd(r)) wirke senkrecht zur Platte in Richtung 





’ 
P(P) 
d x(r) 
+ ZZ 
u 
A 655.1] -z 1 -mw 
Abh 1. 
fallender z (Abb. I). Die Dicke der Platte sei ö, ihr Trägheitsmoment ‚= 0*/12; ihre 
2 .) \ 
no ‚m 2m 
elastischen Konstanten E,G,m|E'=E =G 
m’ — | m l, 
Die Konstanten des elastischen Baugrundes seien E,, @,, m, ohne die absolute 


Größe der Querdehnungszahl zunächst festzulegen, sei zur Vereinfachung angenommen, daß 
m m,. also auch E/E, = @/@G,. Die Spannung zwischen Platte und Unterlage sei g=q (x). 
N sei die Bettungsziffer der nachgiebigen Unterlage. Im achsensymmetrischen Fall 
tritt an Stelle von .r der Radiusvektor r, an Stelle der Streckenlast P’ die Einzellast P. 


«, w sind in beiden Fällen die Verschiebungen in den Achsenrichtungen. 


3. Die Normalspannung zwischen Platte und Baugrund. Mit Hilfe des Fourierschen 
Integraltheorems definieren wir die den Funktionen 8 (r), p (#), q (#) zugeordneten Funktionen 
VW, 7,0 dureh die Gleichungen: 


IK Hayashi: Theorie des Trägers auf elastischer Unterlage, Berlin 1921. 


?) Hierüber und über weitere Literatur vgl. Handbuch der Physik VI, S. 180 u. 224, Berlin 1928: insbesondere ist 
dort auch über Verbesserungen des Zimmermannschen Ansatzes berichtet; siehe auch Prager: ZAMM 7 (197), 
S,.354, und die Arbeiten von Schleicher. Neuerdings hat M. A. Biot (J. appl. Mech. 4, 1 bis 7 [1937 |) die in der 
vorliegenden Arbeit untersuchte spezielle Frage in ähnlicher Weise behandelt (Anm. während der Korrektur). 
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zo) 
I .s 
-_ \we öwdE \we os&rdg 
mir) 7 IE): -cosex2qa& - (£)-coses.rd se 
fe 0 
(3.1), 
‚(.r) . \r S)\ecosörd: 4) AUITVRPET 
} . 7 Ds I n 6 oo. I! ne e En 8 oo. u 


h u 
wobei sehon benutzt ist, daß die Funktionen w,p,g gerade Funktionen in x sein sollen. 
Die Differentialgleichung der Balken- (Platten-) biegung: 
E' J w!V - p(a)+ g(a)=0 
(p,g hier als Zug positiv eingeführt) geht damit über in 


EJIWIE)- PIE HOY(H=O 
oder 


(S) (£ 0)’ 3 (1 m (£ 0)” E' 


= PR > Pin PR ‘).) 

P()-0(9) 1 6m-—1) P(9-0(9)12% | ne ER). 

Im elastischen Halbraum herrscht nach Voraussetzung ein ebener Formänderungs- 

zustand; es kann daher die vom Verfasser an früherer Stelle eingeführte Spannungsfunktion*) 

D zur Bestimmung seines elastischen Verhaltens herangezogen werden. Man erkennt leicht *), 
daß der Ansatz: 


I 
I (/ u ur ds 
D \(mze:+\ E es) eos &u Be ea emo Fr 
A e 4 
0 





mit den willkürlichen „Konstanten* M=M(&, N = X (£ zum Ziele führt. Rechnen wir die 
a.a.0, benutzten Lameschen Konstanten gleich auf die hier benutzten um, so erhält man’) 
[für die Spannungen und Verschiebungen in der Berührungsebene zwischen Halbraum und 
Platte (20) 





Il m A z E 
36, 0,2 (X%,0): - \(a- = „Nleos&xds, 
v0 
| I f 2 m N. R 
3, Ms 3 = „sine. dE. 
0 - 
” (3.4) 
L([ m i de 
we, (2,0) = \ [2 M A cos &ır 
T. m—2 | & 
I(/[ m m \. sd 
0) \(,, 3M+ — N)sin Ex 
J ; y y } & 


Als UÜbergangsbedingungen zwischen Platte und Halbraum kann man nun die zwei 
Gleichungssysteme ansetzen: 


(elastischer Halbraum) (nachgiebiger Halbraum) 
=», 090,=0, =0 q=0,,=Kw 
oder, für die zugeordneten Funktionen geschrieben: 
m ."s6m—-1P(&$) —-0(d 
m — 2 B G m E00) a BD). 
I2KOPIS) VS) 
M— —,N=5.'00), k" z9)" 


2M+mN=0, 


2M 


JE) 





>) K. Marguerre: Spannungsverteilung ... . Platte. Ing.-Arch. 4 (1993), S. 332. 
ı) A. a.0.,S. 339. 
5) A.a.0. Gl. (2.54). 
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Man erhält also den Druck @(£) zwischen den Medien als Funktion der gegebenen Belastung 
P(:) in der Form: 


P(&) P(£) 
VE) (4 Q(E = ER’ 5 re 
® & $13 e E A 
l- 66," )) It jox5° 


und daraus folgt, wenn wir durch 


/2EI_ 1/6 #7 E' 
) m -=dÖ 8.7 
ı | 7? | 66, ; | Ei; * | I? Ko 3.0), 
zwei „Bezugslängen“ /, und IT, einführen, und zur Abkürzung setzen 
l Z 2. Hu 
Sı el, 2 I. ’ Sr &l,, u l. 
I rs 
rg ee GA I, ine I 1 EA Ni Till (38) 
fı er I, 7 + &°® Il» “en l, 7 1+£,? . . . R 
0 0 


Tritt an Stelle der Last p (=), speziell die Schneidenlast P’ längs #—=0, so ist‘) P(£) 
einfach zu ersetzen durch P’ (vor dem Integral); die Aufgabe ist damit auf die Diskussion der 
Integrale 


’“ an SL 
ER l Sm E e ds, ts) 1 T-- tn. R. kon &. “ „2 (3.9) 
7 I+£, an) 1+% An. lI+5, 
0 v — ı 


zurückgeführt. 


4. Diskussion der Integrale. Vergleich. Die Auswertung der Integrale (3.9) erfolgt zweck- 
mäßig mit Hilfe des Residuensatzes in einer komplexen = £+in-Ebene: 


| oil, X dE, : | " oifaxa d &, 
0,(8,) = Me r- \ i+# (2 (X) Re, * \ Be + (4.0). 
0 — 60 


Das Integral 7, läßt sieh (in Übereinstimmung mit bekannten Ergebnissen) dabei in 
geschlossener Form angeben, für q, erhält man einen geschlossenen Ausdruck und ein 
numerisch leichter auswertbares Integral. Die Pole (Unendlichkeitsstellen des Nenners) der 
Integeranden liegen bei: 


| 4 6 ET | art a | 
mr - 2 e” 2 ] 3 . Aa 1, ag v2 (li), u ) ( l+i), 
> (3) | ! +) = 3) 1 z ö = (y l . ui 

Sı 2 2 ı», A v2 l ), 2 =yg i). 


q, läßt sich, wegen des Poles bei —1, nicht, wie q, längs der negativ reellen Achse 
integrieren man wählt hier zweckmäßig den Weg über die positiv imaginäre Achse. In 
beiden Fällen kann man den Integrationsweg durch einen großen Kreisbogen vom Radius R 
in der positiv imaginären Halbebene schließen. Das Integral über den Kreis verschwindet 
mit R>v;: man hat daher: 





N 
(x " a > >» (1) i e koche d ıı — yr > > - (DD . ”» (2) 
Ylr,) Re [2 es (I, )+ ey Tin 4: (,) = Ne [i Res (£,’) +i Res (£,®)]. 
Bu; 
re BE in D(Ö 
Da die Pole des Integranden V(d einfach sind, ist das Residuum unmittelbar durch N’ (ö 


vereben; man erhält also: 


6) Vgl. z. B. Anm. S. 225, a. a. O., S. 337 unten. 
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. A M ei-ı'® Xı l » P N; Xı (i u,” a ( j elsa'’ Xo | ei ee 
4,(%,) „HEIL 5 7# (ma 3 an; (2: \0,) = NE | - ma TVo /= 13 
5 1 3 “1 ) JT R I+,' S(C, ) ll, | 
0 
E % a Ks Bu 1 Sr 1 ” 7 2. u" 
_—e 3 v3cos— tsınz : V2e '*’lkcoso > +tsn zz)? (41). 
, \ _ ze 4 \ | - ] — 
S- 
l [e-"*Xın’dn, 
u l+n,° 
0 





Die Frage nach dem oszillierenden Charakter der Lösungen ist auf Grund der Formen (4.1) 
leicht zu beantworten. Bezüglich 4, ist keine Erörterung nötig; bezüglich q, erkennt man zu- 
nächst, daß der oszillierende (Residuen)anteil bei x, =0 überwiegt. Denn das Integral lautet 
für 2, —=®%: 


1 'n ’dn, Il ndn 1 
. 0) = ® . Pe er ae Fr a 2): 
= \i fm aRligm u (4.2) 
0 0 
es ist also 
=. dE, | ; 2 E ,.. a 
Aloe 53 V3-9V3=5V3 =0,385, 
»i ® ® 


was man übrigens auch unmittelbar durch elementare Integration bestätigen kann. 


Man sieht aber auch, daß der oszillierende Anteil bei seinem ersten Minimum 
absolut genommen größer ist als das monoton fallende Integral; das Minimum liegt an der 


A 1 ee >  - x 4 
Stelle, wo die Ableitung — —e ? 1 3sin 5 + cos „| zum ersten Male verschwindet, also 


3 ya y 
0 ) 
bei - - —— d. h. 
2 6 
(i) „ pr )* r 4 [7 
2, 5 1=b523 A RR Te: De © 


Der Betrag der Funktion an dieser Stelle ist 








| > 7 V3 l f ‚4 
se’ ze 3.522.103 
Für das Integral aber gilt: 
” F- 
| eo 1X,” | | u | 
7 \ l+n,° dn,< 7 \e n1%Xı u,” d 7; de e' 256 - 10 3 
0 A m 


Der oszillierende Anteil entscheidet also das Vorzeichen. 


Die Lage der Nullstelle von g, (x) läßt sich graphisch (d. h. näherungsweise) ermitteln; auch 
für diese graphische Auftragung des Verlaufes von 9, (x,) ist die Zerlegung des Integrals (3.9) 
mit Hilfe des Residuensatzes sehr vorteilhaft. Denn die direkte Auswertung ist, auch ab- 
gesehen von grundsätzlichen Schwierigkeiten, für große «-Werte wegen der heftigen Oszillation 
des Integranden mit £ sehr unbequem. Das Integral J, aber läßt sich gerade für große x 
sehr leicht mit Hilfe der Simpsonschen Regel für das Intervall zwischen O0 und einer 


ie 0) 
geeignet gewählten oberen Grenze a=a(x,) auswerten, da der Beitrag, den \ liefert, sehr klein 
a 


Ist. Man kann diesen Anteil übrigens (worauf hier nicht näher eingegangen werden soll)abschätzen, 


3 


’ 1 
t rs ; näherungsweise ersetzt durch Bart durch den Integrallogarithmus®) 
N Iı . 


in geschlossener Form ausdrücken läßt. 


da er sich, wenn man 


Ein unmittelbarer Vergleich des Gesamt-Verhaltens der beiden Unterlagen ist, schon 
weil die Bettungsziffer X [kg/em’] und die Elastizitätszahlen [kg/em?] verschiedene Dimension 
haben, nicht möglich. Wohl aber kann man die beiden g(#)-Kurven in ihrem typischen Ver- 
lauf vergleichen, indem man z. B. auf gleiche Maximalspannung 9, bezieht. 


8) Vgl. z.B. Jahnke-Emde: Funktionentafeln, S. 78, Leipzig 193. 
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Wenn wir die bei der Diskussion der Inteerale unterdrückten Konstanten wieder ein- 


führen. so eılt: 





4,1) T y,| n 4) 2] 1%.) 
pi = =. eo. @\ , p' x 5 j (4.4). 
Ur LE emand 7 Face ’ 9 ia, sgre Yılhlosoa tina) 
er -] ,. ) = ) tg 
wie es sein muß, in Übereinstimmung mit bekannten 


(Das Resultat für 4, befindet sich, 
rzebnissen.) Mit Hilfe dieser Gleichungen lassen sich die Parameter 1,. I, ausdrücken durelı 


das Verhältnis von Last P’ zu Maximalspannung 4. in der Form: 


pP' „pP 
l,- 0,385 . =0,358 ee ee en A 
Im Im 
ls wırd also: 
- ah; .)=.) er * 
My 0,589 .r, Mr (2.9) X, , : : : ( 1.6). 
Im Im 
In Abb. 2 ist daher », gm I’ 
0.385 .r, 0.353 .r, 


als Abzisse gewählt. Ordinate ıst g/m 


I; ; 
4" pP’. 0,385 12" pP. 0,353° 
T. 
Die Gleichgewiehtsbedingung | y(.r) de = P’ geht in diesem Koordinatensystem über ın die 
SI 


(Kontroll-) Bedinzung: 
De 


gm da, 1 


Für die Nullstellen ergibt die numerische Auswertung der Integrale (3.9) 
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und beinachgiebigem — Verhalten des Baugrundes. 


Verlauf der Normal-Pressnunge bei elastischem 


\hh. ?. 
In Abb. >? ist der Vollständiekeit halber noch die nieht oszillierende Funktion 


L 


Ya 23 (cos? x, 5 |leRinmsdeE [e 200 (#, 0) PEN 
0 7 \ | + & u T in. #3 le r2!xXo ( r,<0) ee 
[m nu SER Ts er), - Ho 


aufgenommen. Der Faktor 2 im Exponenten ergibt sich aus der Forderung (4.7), der die 
wenn sie vernünftigerweise vergleichbar sein soll, genügen muß. 


Funktion 4,. 
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Der oben rechnerisch nachgewiesene oszillatorische Charakter der Funktion (3.8), der 
unabhängig ist vom Verhältnis der elastischen Konstanten, läßt sich auch leicht physikalisch 
einsehen: Die Annahme verschwindender Schubspannung zwischen Balken und Unterlage, 
die auf den Maximalbetrag der Normalspannungen fast olıne Einfluß ist, beeinflußt den 
Verlauf entscheidend. Wenn keine Schubspannungen angreifen, müssen nämlich die 
Normalspannungen oszillieren, weil sie sonst den Balken als ganzen verbiegen würden: die 
Enden* würden eine von Null verschiedene Neigung aufweisen (s Abb. 3), was offensichtlich 





Abb. 3. Typischer Verlauf der Balkendurceh 
biegeung unter Einzellast und elastischem Gegen 
druek: a,e) bei niehtoszillierendem Druck ohne 
und mit Riüekbiegrung infolge der Tangential 
kräfte; b) bei oszillierendem Druck. 


mit den Randbedingungen im Unendlichen in Widerspruch steht. Die Lösungen von 
Boussinesq und Michell’), für Halbraum und Halbebene aber zeigen, daß bei Anwesenheit 
von Schub die Normalspannungen nicht zu oszillieren brauchen. 


5. Die Verschiebungen. Achsensymmetrisches Analogon. Die Normalverschiebungen des 
Balkens sind im Falle der nachgiebigen Unterlage mit Hilfe der Gl. (3.5) und (3.5) sofort 
angebbar: 


P'i1l(cosöxde . 
1°, (.,0) 7’ \ (sa ee er 
0 Ei 


Bei der Bestimmung der Einsenkungen, die der elastisch gelagerte Balken erleidet, 
stößt man auf die gleiche charäkteristische Schwierigkeit, die die Michellsche Lösung für 
die Halbebene unter Einzellast bezeichnet’): die Verschiebune wird für große .“-Werte 
logarıthmisch unendlich. In unserem Fall äußert sich das formal darin, daß man ein 
nicht-konvergierendes Integral erhält; aus den beiden letzten GlIn. (3.5) ergibt sich nämlich 


m . m IO(8) 
2M =N= =, Biso: 
m —.& m (1, 
P’ m II cos (Eur) ds r 
mw, (#,0) = \ - (9.2). 


G, moon 


Diese Schwierigkeit fällt beim achsensyinmetrischen Analogon weg; dort findet man: 


0) m—1l 1 \ J,(Er)de ze 
w,(r,0)=— — 2.9). 
Zu, GG m Zn) 1+1?& ass, 
N) 


In diesem Falle lassen sich also auch die Einsenkungen in den beiden Fällen unmittelbar 
vergleichen. Dafür ist aber die numerische Auswertung mit Hilfe des Residuensatzes (Bestimmung 
des Verlaufes mit r) schr erschwert, weil es für Real- und Imaginärteil der Besselschen 


‚Verl. z. B. Love-Timpe: Lehrbueh der Elastizität, S. 235 und 250, Leipzig 1907. 
») Vgl. z.B. Love-Timpe: A. a. O. S. 250, 
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(bzw. Hankelschen) Funktionen vom Argument r (1/2 +i/2y5) keine Tafeln gibt‘). Da 
erundsätzliche Verschiedenheiten gegenüber dem ebenen Problem aber auch nicht zu erwarten 
sind, soll von dem Versuch eines zahlenmäßigen Vergleichs Abstand genommen werden. 


Interessant ist dagegen der Vergleich der Maximal-Durchsenkungen und -Spannungen, 
die sich elementar errechnen lassen. Man erhält: 


., | L 
00 Pm—-ii |] dE, | (0.0) P1]1 \ &,dE, 
m m, (VU) - y r W,mT W\VUVU) = 7775 er 
1 m 1 GG, m I, 22)j1+2° 7 ‚ Kl22a)1l+3&, 
0 0 
0.985 we —ı P | | P 
dm) e y . | BR Be r Du | 
; nm 2G,l, | S Kl} 
| 
: 0" : | r Ö° ’’ 
weren (Gi, 61; (+ also: | wegen K: 121, E’ also: 
(5. ı,’p | sm-—1l 12 P N 
) .). . rg | m, . in 7 . . e),: 
ı ı m ‚ie m N’ (\ | em t m Ay ({ ( )» 
v. | 00 
Äh PıX 262, P1i1[2,85, 
ferner: YA 1:2 a | +83 ferner: Pe 2a \ IE 
0 0 
P P 
0,355 - 0,25 (5.5) 
ur; Fr u 
21, 2l, 
Kliminiert man die Parameter ! aus den Gl. (5.4) und (5.5), so ergibt sich: 
| m—1 P? m--1 P* 
( :0.221.: . dUms = :V.094 : — 2.6). 
Im ı We m (0° Im: We ° m  (0" a 


Bezogen auf gleiche Maximal-Durchsenkung »v „ ist also der Maximaldruck beim elastischen 
Halbraum über doppelt so groß wie beim nachgiebigen. Umgekehrt verteilt sich also die 
Pressung über den elastischen Halbraum sehr viel weniger weit. 


6. Die maximalen Biegespannungen. Von besonderem praktischem Interesse ist der Vergleich 
(ler maximalen Biegespannungen in der Fundamentplatte, wie sie sich aus den beiden Hypothesen 
über Verhalten des Baugrundes ergeben. Die größten Spannungen an Ober- und Unterseite 
der Platte ergeben sich nach 
m E' 


' » 
(G ww 0=75 


mw” ee 
m 1 Er er 


Om 


(ebener Fall) aus den Gln. (5.1) bzw. (5.2) in einfacher Weise: 





SL N 
er (( lecosäörfds pr E' 6 kr i &dE 
ee = +1? mE En) Ir 
Ü Ü 
(6.1). 
LI Rn 
Pd5G 1 Vers 0,6, d &, P’öE l jr E I. Pig d F 
u 98 1+£/? WERE 3 1 + £,' 
0 0 
Wegen der Gl. (3.7) ist daher das absolute Maximum (für «=0) gegeben durch: 
P'6l, „ 2331 Pl, P'6l, 1 2,12 P'l, ' 
0 0,385 = On e — er 
er 0 0° Ö° 0 0 2Yv2 O0” 


Führt man an Stelle der Parameter ! wieder das Verhältnis zwischen Last und größter 
Normalspannung q,„ nach Gl. (4.5) ein, so erhält man: 


»’2 pP’? 


Ö —()SUJ. = OÖ 0.7D- = . ‘ r ö . . r ö 6.3). 

m N &. m Ay E Ym ( ) 

Die gefährliche Biegespannung im Fundament liegt also bei elastischer Unterlage um 

rd. 20°/, höher als bei nachgiebiger, wenn die maximale Pressung in beiden Fällen als gleich 
anzenommen wird. 


10) Betr. das Argument r(1+i) vgl. Jahnke-Emde: Funktionentafeln, S. 296 ff., Leipzig 1939. 
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Im achsensymmetrischen Fall wachsen die Biegespannungen unter der Last über alle 
(„renzen: 


n LH 
”w, 129w Pm—-iIl 1 “in dE, 0? mw Iodw PıI | (& dE, 
or "rdr Go m |’ 2n)i+E£? or? roör K |’ 2n)1+2£,'' 
v 0 


Es ist aber interessant, daß das Verhältnis beider Ausdrücke unabhängig von allen 
physikalischen Konstanten den Wert 1 hat. 


Zum Schluß mag die Gl. (6.3) noch mit den Ergebnissen einer für technische Zwecke 
meist ausreichenden Näherungsrechnung verglichen werden. Ersetzt man den Einfluß der 
Bettung auf die Biegespannungen durch den einer verteilten Dreieckslast von gleicher Gesamt- 

! 
intensität (| g (x) dx = qm:1==P’) (vgl. Abb. 4), so erhält man für das Maximalmoment 
a 
M „x 4m: 17]6, also für o,: 
P I 
u u un ER » 
u Et „1 N (6.4). 


p' 











Abb. 4. Ersatz der Bodenpressung nach Abb. 3 durch eine Dreieckslast. 


Mit Hilfe der Gleichgewichtsbedingung für die vertikale Riehtung läßt sich I durch 
max Ym ausdrücken, und man findet: 


p'? 
Ö’ Gm 


Om = (6.9), 
ein Ergebnis also, das größenordnungsmäßig vernünftig ist und außerdem auf der sicheren‘ 
Seite liegt. 


In den Anwendungen wird allerdings q4,. ebensowenig unmittelbar bekannt sein wie I; 
man kennt nur (meist auch sehr näherungsweise) den Modul E, bzw. die Bettungsziffer Ä, 
d. h. in unserer Schreibweise die Parameter !, und !,. Wählt man nun nach Gl. (4.6) 





37 a 
I: 0-1, l N | 2 l, (6.6), 
d. h. ersetzt man die Kurven der Abb. 2 durch Dreiecke von der Basis 2.*, so kommt: 
P’'I PL ah 
Om 4,0 N? Om = 3,39 N? j . . . . . . . . . (6.4 ). 


Vergleich mit Gl. (6.2) zeigt, daß auf diese Weise die Näherungsrechnung, vor allem für 
die elastische Unterlage, zu ungünstige Ergebnisse liefert. Gl. (6.5), die sich auf gleiche Höhe 
{m des Dreiecks bezieht, ist besser, wenn auch der für wirklichen Baugrund zu wählende 
Zahlenwert wohl genauer zwischen 0,50 und 0,55 liegt. 


7. Zusammenfassung. Ein durch eine Einzelkraft belasteter langer Träger ruhe auf einem 
nicht-starren Baugrund. Es wird der Verlauf der Pressungen zwischen Träger und Unterlage, 
wie er sich aus den beiden Annahmen „elastischen“ (Hooke) und „nachgiebigem* (o. x m) Ver- 
halten des Baugrundes ergibt, (die unmittelbare Gegenüberstellung der Auswirkungen beider 
Hypothesen ist die hauptsächlichste — methodische Absicht der vorliegenden Veröffent- 
lichung) bestimmt, ferner die Maxima der Einsenkungen und der Biegespannungen (ebener 
und achsensymmetrischer Fall). Insbesondere zeigt sich, daß die Normalspannungen in beiden 
Fällen oszillierend abklingen, nur befindet sich im Falle elastischen Verhaltens des 
Baugrundes die erste Nullstelle in größerem Abstand vom Lastpunkt, und die Spannungen 
entgegengesetzten Zeichens sind viel weniger ausgeprägt. — Der physikalische Grund für das 
oszillatorische Verhalten ist die Annahme verschwindenden Schubes zwischen Träger und 
Baugrund, 655 
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Über die Knickung eines Balkens durch Längskräfte. 


Dem Andenken an Erich Trefftz gewidmet. 
Von Otto Blumenthal in Aachen. 


m ‚Jahre 1598 hat L. Vianello'!) ein Näherungsverfahren angegeben, um die Knicksicher- 
I heit eines dureh ein System von an verschiedenen Stellen angreifenden Längskräften be- 
anspruchten Balkens festzustellen. Man geltt aus von einer wahrscheinlichen Form der 
Biegzungslinie und unterwirft diese einem Iterationsverfahren. Die überzeugende Kraft von 
Vianellos Abhandlung liegt in einer Reihe durchgeführter Beispiele, bei denen sich heraus- 
stellt, daß das Verfahren auch bei recht roher Ausgangskurve in ganz wenigen Schritten zum 
/iel führt. Demgegenüber möchte ich hier zeigen, daß dem Verfahren nicht unbedingt ge- 
traut werden kann, daß es nämlich in gewissen, nicht trivialen Fällen entweder schlecht 
konvergieren oder sogar zu verkehrten Resultaten führen kann, die eine wesentlich höhere 
Knieksicherkeit liefern, als der Balken tatsächlich besitzt. Wie aus dem Folgenden lıervor- 
sehen wird, liegt dies daran, daß unter Umständen die dem kleinsten Knickwert ent- 
sprechende Biegungslinie eine unerwartete Gestalt hat. 

Diesen Beweis führe ich in $ 2 an einem Beispiel. Im $ 3 entwickele ich dagegen 
einiee hinreichende Kriterien für die Anwendbarkeit des Verfahrens, die Vianellos Bei- 
spiele ın der Hauptsache deeken. Zum Verständnis muß ich in $ I die von E. Trefftz?) 
entwickelte Theorie der Erscheinung vorausschicken, die auf einer linearen Integralgeleichung 
beruht. Ich kann dabei die Trefftzschen Entwicklungen ein wenig abkürzen. 


$1. Zurückführung der Aufgabe auf eine Integralgleichung. Der Balken habe die Endpunkte 
V und 7, die Abszisse sei mit s, die Ordinate der Biegungslinie mit % bezeichnet’). Die 
Verteilung «der Längskräfte wird nach dem Vorbild von Trefftz stetig angenommen, von 
der Diehte pts), wobei positives p(s) im Sinne wachsender s gerichtet sei, außerdem wirke 
eine stetige (uerbelastung von der Dichte q(s)'Y). Auch Elastizitätsmodul # und Träg- 
heitsinoment ./ können wir, ohne Komplikationen einzuführen, als stetige Funktionen der 
Abszisse voraussetzen, sie sind überall positiv. Das Biegemoment sei m(s). Als Rand- 
bedingung werde an jedem der beiden Enden eine der folgenden 4 „natürlichen* zugelassen: 


a) Stützung: yzm=U), 

b) Einspannung: ye=y -d, 

e) freies Ende: m = m’ = 0°), 

dd) Einspannung in eine in senkrechter Richtung verschiebliche Backe: „= m’=V. 


Sei /ıs, f) das von der am Punkt f angreifenden senkrechten Einzellast 1 am Punkte s 
bei den gegebenen Randbedingungen hervorgerufene Biegungsmoment. Es ist 


| en M, (), Ss S ! | 
Mis, ea Ehe ce are en 
| | 9,8 —(s—ih) st | 


0 


wo), und 9, Einspannungsmoment und Stützkraft am Punkte O0 bezeichnen. AM (s, t) ist also 
eine stückweise lineare Funktion von s. Kompbinieren wir also 


d’y m (Ss) 


) 
ds? E.J (1,2) 


lt 
d’ Mis, A 


ds? 


0 s=t, 


so erhalten wir nach dem bekannten Verfahren des Greenschen Satzes 


I 
m 


E.J Mis,t)ds 


d N () ds 


dy dMis,t\ dy dis, \ 
! 


(15 ae |, te Da,5 y 


N 

1, ZVDI 42 (1808). S. 1436 bis 1443. Das Verfahren wird auch anderwärts benutzt, siehe z. B. E. Pohlhausen: 
Berechnung der Eigensehwingungen statisch bestimmter Fachwerke, ZAMM 1 (1921), S. 34 bis 2. 

2) ZAMM 3 (1923), S. 272 bis 275. 

3) Zwisehenstützen, die Trefftz zuläßtz lassen wir zur Vereinfachung der Formeln weg. 

ı) In der Voraussetzung stetiger Diehten p(s) und q(s) sei eingeschlossen, daß an den beiden Enden keine 
konzentrierten Kräfte wirken. Das vereinfacht einen Schluß. Die Endformel (1,6) giltaber, wie man dureh 
Grenzübergang zu konzentrierten Kräften leicht zeigt, ohne jede Einschränkung der angreifenden Kräfte. 

5) m’ 0 kommt heraus, weil nach unserer Voraussetzung (Fußnote %) an den Enden keine konzentrierten Längs 
kräfte auftreten. 
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und daher, weil die Randglieder wegfallen und weil nach (1,1) 


dMis, h d Mis,t) 
| ds ) | ds j 
ist, 
l 
" mis) . 
vW=\yyg nd ee 
0 


Andererseits berechnen wir, unter Berücksichtigung von Fußnote *), das Moment mm (s) 
aus den Längs- und Querkräften (m, Einspannungsmoment am Ende 0): 


m (s) == m, + \iy(s) — y(o)y p (0) do \(s - o)g(o)do. 
0 0 


Daraus folgt: 


l 
m’ ()= , (u (s) Ps) 3 ee 
mit 
P(s) : \p(o)do > A KH ON ah SEE Dr Ve 4 Y5,9 
0 


Da aber bekanntlich — m’’ (s) diejenige Dichte der Querbelastung ist, die das Biegemoment nm (s) 
hervorruft, so ergibt sich aus dem Prinzip der Addition der Querbelastungen 


er I I 
l | 
m (s) \m’in Mis,Nat: \.,veor@M (Ss, At+\g(d Mis,Nat. 
0 Ö 0 


Das letzte Integral ist das von der gegebenen (Querbelastung herrührende Moment m, (s). 
Das erste Integral wird durch Produktintegration umgeformt, wobei die ausintegrierten 
Glieder wegen P(V)= P(l=0 (siehe Fußnote *?)) verschwinden. Also resultiert 

! 
oO Ms, t) 


m()=\y(d) ri T 


M) 


Bu ea re A 


Setzen wir hierin für y’(t) den Wert aus (1,3) ein, so kommt 





I I 
\Mıls. "mÄr) Od Mir, h) 
m(s) M,\S) -\ Pin Of d \ En. of dr 
M 0 
(1.6). 
l I 
m(r) 5,0 Ms, do Mr, t) 
ms) + \ En \ A ne rl. 
0 0 


Dies ist für m(s) eine lineare Integralgleichung mit symmetrisierbarem Kern. In der 
Tat, führen wir die neue Funktion 
' m(s) 
n(s)= 5 
Y EJ(s) 


ein, so genügt diese der Integralgleichung 


I 
m.,(s) i ” 
RAE teen mar Ge Gr NE 
0 


mit dem symmetrischen Kern 
I 
0 M(s, 0 Mir, t) 


\r Ben N 


((r,s): A’ A 


] E J (r) E .J (S) a 
0 
Die zugehörigen Eigenfunktionen genügen den Gleichungen 
l 
y;(s) k;\wy;(r) G (r,s) dr a 
0 


wobei die 4; nach der Größe ihrer Beträge numeriert sein sollen. 
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Das Vianellosche Verfahren besteht darın, daß man, ausgehend von einer plausibelen 
(0) 


Momentenlinie n‘” (s), die ihrerseits aus einer plausibelen Biegurgslinie y‘” (s) gewonnen 
wird, das Iterationsverfahren anwendet: 





l ! | 
n”’(s)=\n"(r)G(r,s)dr, nn? ()=\in"’(v)alr,s)dr=\n"(r) 6? (r,s)dr, 
hr () u" 
( 
| | (1,9), 
"(lin kh)alr, drin" Ir)" lr,s)dr,... 
() 0 


wo G') (r,s) den kt" iterierten Kern bezeichnet. Es habe nun die Ausgangsfunktion in den 
Kigenfunktionen y die (formale) Darstellung 


n‘” (s) UI PALIN IN FT Hr a er es (1,107), 


womit ausgedrückt wird, daß die ersten Koeffizienten e,,...,ei—ı verschwinden können. Wir 
wollen außerdem, um ganz allgemein zu bleiben, annehmen, daß y'; (s) und yi; +1 (s) zum gleichen 
Kigenwert /; gehören®). Dann wird 


C; y; (S)—+ Ci 4 ıy; +ı[(s) U; ı>»(S) 
rı j + ri yı Z 
nk) (s) Be n ug = LT a ui: a R . . . . l 10 
A;h ee ERAOR: 
Ist jetzt |; <A: +2, so Ist das Iterationsverfahren (1,9) konvergent und ergibt 
(}) r 
h nn’) (s) 
lım i: 1,10’ 
k>« nik l I) (s) ı . ( . ), 
a ’ BE a 2 \ l 
während gleichzeitig lim A" nn (s) = e;y;(s)+ ei +ıyi-+1(s), bis auf den Faktor _,_, ‚ die 
kK>& ] E .J/ (s) 


Momentenlinie für eine mögliche Knickfigur des Balkens darstellt. Die Bedingung |4; |< |4i + >| 
ist sicher erfüllt in dem wichtigsten Fall, wenn Pf) sein Vorzeichen nicht wechselt, also 
über den ganzen Balken weg entweder nur Druck oder nur Zug herrscht. Dann ist nämlich 
der Kern @ (r,s) nach Trefftz (S. 274) definit (bei Druck positiv definit) und hat also nur 
Kigenwerte eines Vorzeichens. Bestehen aber Eigenwerte verschiedenen Vorzeichens, so können 
wir die Konvergenz des Iterationsverfahrens unter allen Umständen dadurch erzwingen, daß 
wir zum iterierten Kern G@° (r, s) übergehen, d.h. nur die Funktionenfolge n“"” (s), n“” (s),..., 
n 2) (s),... betrachten. Es gilt immer 

h n'') (s) 

lım einzel. 

k>oNn“T-)(S) 

Die gesuchte Knieksicherheit ist der kleinste positive Eigenwert. Sollte also für eine 
Verteilung Pit) der Längskräfte ein positiver Eigenwert vorhanden, aber der absolut kleinste 
Kirenwert neeativ sein, so ist das Vianellosche Verfahren im allgemeinen unbrauchbar. 
Aber diesen Fall wollen wir hier beiseite lassen. Es sei also 4, positiv. Dann ist eine evi- 
dente notwendiee Bedineune für die Brauchbarkeit des Vianelloschen Verfahrens i=1, 
d.h. die Ausgangskurve »"” (s) muß Komponenten nach den zum ersten Eigenwert gehörigen 
Kigenfunktionen aufweisen. Darüber hinaus müssen, damit das Verfahren rasch konvergieren 
soll, diese Komponenten verhältnismäßig groß sein. 

Da das Fehlen dieser Komponenten in n“” (s) eine Ausnahmeerscheinung sein wird, 
kann man saeen. daß das Vianellosche Verfahren im allgemeinen den Knickwert liefert. 

Eine große Schwierigkeit und Unsicherheit kommt aber dadurch herein, daß schon bei 
sehr einfachen Beispielen die üblichen anschaulichen Merkmale dafür, daß n‘” die erste Kom- 
ponente überhaupt oder genügend stark enthält, völlig versagen können. Dies soll im 
folgenden erörtert werden. Die Methode Vianellos verliert dadurch nach meiner Auf- 
fassung erheblich an praktischem Wert. 


0 und / gestützten Stab von unveränderlichem E.J greife im Punkte a die Längskraft + P, 
im Punkte b(> a) die Längskraft P an. Die Ordinaten der Biegungslinie an den Punkten 
a und 5b seien mit « und 5 bezeichnet (Abb. 1, S. 237). Dann findet man für die Momente 


$ 2. Ein Beispiel, in dem das Vianellosche Verfahren versagt. An dem in den Punkten 


A 


Feld Vs. a) M, = Pe 


Ss, 


6, EB. Trefftz hat mieh darauf aufmerksam gemacht, daß mehr als 2 orthogonale Kigenfunktionen zu einem 


Kigenwert nieht existieren können, weil sie Integrale einer linearen homogenen Integralgleichung 4. Ordnung sind, 
Daß andererseits ? Eigenfunktionen zu einem Eigenwert wirklich auftreten, beweist ein Beispiel des $% 2. 
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I —aqa B—a 
Feld Illaz=s=b) Mı= .P’ s+P(y—a)=P l—s)+Piy-P) 
) Aa l A ! > 
PT ls 5) r(y u 
 — a 
Feld II B<s<D Mu= P' , 
Berücksichtigt man die Randbedingungen y (0) y(l)=0 und benutzt für Feld II der Symmetrie 
halber die letztangegebene Form des Momentes, so ergeben sich mit der Abkürzung EJ y? 
die Gleichungen der Biegungslinie: 
) as" 
Feld | yı(s) gt 3 +4As, 
: Bern I a-+ b\ | a—+b 
Feld 11 Yy11(8) = ud. . s 5)+ B cos ‚(s = r C sin »|s a, 2,0). 
„E-all—s) 
Feld II Y(s) ——y! z Dis) 





Die Integrationskonstanten A, B, C, D bestimmen sich aus den Bedingungen 


la) Ynla)=a, Ynlb)= Yyn(b)=P. 
Wir erhalten: 
A „p a a” Pr p Pr. all I)? en 
uns ans dh zu: BE 2. Dt) 
! I -(b-+ ,. b-a I—(b a) 
Beos» -—- —=(ß —a) 91 Em, Csinv» (Ba) 7 2.22). 


Die Gleichungen für stetigen Übergang der Tangente, (a) = y’ı(a) und Yy(b) = y'n(b), 
lauten nach Einsetzen aus (2,11), (2,12) 


: ‚PB-aa a B-a,| u ba, CM b—a 2.29) 
v“ = az Z — bvsSsınv "7 UvVCcOosvy .uo 
ı 3 a I 2 2 
RB all b)? ’ p 7 ü h da, h Ad 
ya! ] 307 er Byrsuy —. rl —a— . u a. 


Die vier Gl. (2,2) sind homogene lineare Gleichungen für die Unbekannten B, C,a, ß, aus denen 
sich durch Elimination die Gleichung für » ergibt: 


i—! 
Fi)  2+ f u nr - )eos v(b a) 











(2,3). 
mt. E u I | 2 a ,d . tag a 
- - + - sin v a 
Fi Ta ae USE’ ß 
Wir suchen die kleinste Wurzel », d.i. den Knick wert. 
Sei zunächst a=1-b. Dann gibt es genau eine Wurzel», die ae Ist. 


In der Tat ist F'(0) positiv, dagegen 


af 7 e | a jS h er 022 M 
F|; ;) 2 Er )<0 (wegen a==l—b). 


Also ist mindestens eine Wurzel in dem Intervall vorhanden. Um zu zeigen, daß nur eine 
da sein kann, genügt der Nachweis, daß an jeder Nullstelle dieses Intervalls die Ableitung 


dF a l h 1 /1 l l | a? (l—-b)?\ 
— ee . 
dv (' b r a (a) y? | a Ih r all .) r > f b r a 


/ 


ga l I v/ a? (lb)? 
| —+ , BEN 
+| y r R 73 (l ») 3 f bh a ) (b —a)cosv(b — a) 


sinv(b - a) 





negativ ist. In der Tat, setzen wir die positive Größe 


7/1, 1 I(® d-y_ 
ir tr a (l „)+ 3 f Er da )=P 
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und beachten. daß der Koefftizient von cosr(b - a) bis auf den Faktor (b a) mit dem 
Koeffizienten von sın v(b  a)ın Fir) übereinstimmt. so erhalten wir dureh Elimination dieses 
Koeffizienten vermöre Fir)=0 nach kurzer Rechnung 

d F' dA I 7 | 


(bh (t) r -2>cosr(b dt) 


‚sın vr (fl a)<od. 
dy ——. A sın vr (b (dt) I ie 


Verwickelter liegen die Verhältnisse für a=1 b(„symmetrische Beanspruchung“). 


IHliıer ist immer rv ’ eine Wurzel von Fir)=0,. Es kann höchstens eine noch kleinere 
) { 


Wurzel geben. Ob eine solehe auftritt oder nieht, hänet von dem Werte von a ab. Wir 
geben die Diskussion genau, weil sich an sie das einfachste Beispiel für Versagen des 
Vıanelloschen Verfahrens anknüpft. 

Die Gleichune für » lautet: 





ı /2 I ard\ _ 
0=2+2c0osv(l— 2a) + —I—+-;)} a Isinv(l 2a) 
\» dt dr) +) 
2.4). 
I I | | I va\ . / 
eos v| > a)\cosr( a)+| " we =} ” |sin | 2 a) 
r a ' u nn T T . 
Der zweite Faktor ist für #» 0 positiv und hat für » ‚., den Wert 
h dl l ( 
I—-2a/l | 1 Ad | 27° 

N ’ ) 2 2 ıLI[I8 3) 

7 Fr 2.0) I—-23a3 2ararıl FT "a la - | . 1“ 
(2,9). 

| 
2 zarıl I el 





| 
fia) mimmt in dem in Betracht kommenden Intervall O<a< , ständig ab und geht von 


positiven zu negativen Werten über. Die Nullstelle liegt beı 
a -V.3193 1. 


Ist beisymmetrischer Beanspruchung das gedrückte Stabstück verhältnis- 


. . . r . „T .. 
mäßıe eroß (a<a,, dann ist der Knickwert r 9a Ist aber das eedrückte 
rw 
Stabstück kleın (a>a,, dann existiert eın Knickwert vr, 1: Füras a,hat 
ut 
7 
Fir\=0 die Doppelwurzel » ] 
m 


Wır werden nun zeieen, und darın besteht unser erstes Beispiel: Die zu dem 
kniekwert », gehörige Biegungslinie ist merkwürdigerweise nicht symmetrisch 


] ar ; 
zur Achse s = ,, sondern hat bei s° —, eine Nullstelle. Die symmetrische (konkave) 
>’ ne .. JT 
Biegungslinie gehört zu v er, 
l dad 
In der Tat. sei zunächst » — , Setzen wir in (2,24) 1—b=a und addieren zu (2.23). 
A 
} D A 
so erhalten wir wegen cosv = 0) 
[er a | 4 | | N 
I et) = „(pP ayfla): N 
31 m I is —2Z2r?! / 


d.h. #5 a, außer wenn a-=-a,, also v Doppelwurzel ist. Ist aber A=ua, d.h. 


e. er 
ereifen «die beiden Druckkräfte ın gleicher Höhe an, dann besteht (Gl. (2,0)) die Biezuneslinie 


aus 2 zur Achse s , symmetrischen Geradenstücken und einem zu dieser Achse symme- 


triıschen konkaven Stück einer eos-Linie, ist also konkav. 


JT } I, I i 
Ist dagegen r =», TEL dann Ist cos vr , =0 und daher bei !- b= a nach (2,21) 
“Ga y 
» 0. Daher ergibt Subtraktion der Gl. (2,23) und (2,24) «+ P=0. Demnach schneidet die 
Biegungslinie die Achse, und die Gl. (2,0) zeigen genauer, daß sie den Punkts= ,„, y=-0 


zum Symmetriezentrum hat. 
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Für a = a, bestehen Biegungslinien der beiden Typen gleichzeitig. Diese zeigt die Abb. 1°) 
Die Momentenkurven haben die gleichen Symmetrie-Eigenschaften wıe die Biegungs- 
linien. Insbesondere ist die zu », gehörige Momentenkurve nicht überall positiv, sondern hat eine 
I ) 7 
Nullstelle bei. Dagegen ıst die zu du) PER 7 
Daraus folgt: Für a>a, wird das Vianellosche Verfahren, wie man es „ver- 
nünftieer* Weise anwendet, unbrauchbar. Denn man benutzt als Ausgangskurve für die 
Bieruneslinie immer eine konkave Kurve, in der Erwartung, daß der kleinste »-Wert einer 
konkaven Biegungslinie entsprechen werde. Diese Erwartung ist aber füra>a, falsch. Geht 


zugehörige Momentenkurve nirgends negativ. 


Pi 7 

man insbesondere von der genauen konkaven Biegungslinie aus, die » 1-94 entspricht, so 
. a . . » . nd 2 .. 2 4 r T a 

führt das Verfahren immer wieder auf diese Linie zurück und liefert den »v-W ert; Ya 


der nieht der Kniekwert, sondern zu groß ist. Geht man von einer anderen kon- 
kaven Biegungslinie aus, dann kann die Zerlegung ihrer Momentenkurve nach den Eizen- 
funktionen auch die zu », gehörige Eigenfunktion als Komponente enthalten, aber jedenfalls 
als kleine Komponente. In diesem Falle würde zwar wie am Ende des $ I gesagt das 
Vianellosche Verfahren schließlich den riehtigen Knickwert liefern, aber erst nach einem 
langen und praktisch undurchführbaren Prozeß. 

Die Tatsache, daß die zum kleinsten Eigenwert gehörige Biegungslinie nieht konkav ist, 
erscheint zunächst verwunderlich. Man kann sie sich aber mechanisch anschaulich machen 
dureh die Erklärung, daß das gedrückte Stück (a,l--a), wenn es genügend klein ist, sich 
nieht ausbiegt, sondern als Ganzes querstellt. 





PP in 
z 
/ 
/ 
/ı 
’ >ae 
/ Qq 5 { 
F-P 
/ 
/ it 
/ de rn 
/ L - 
/ = 
; a, | 
L / 
/ AP 
u“ 
[A628.2 3] 
P- Pr 
Abh. 1. Abb. 2. Abh. >. 


Die auffallenden Formen der Biegungslinien, die wir bei symmetrischer Beanspruchung 
lestgestellt haben, treten auch bei unsymmetrischer Beanspruchung auf und sind hier sogar 
noch anschaulicher zu verfolgen. Vorauszuschieken ist, daß konkave Biegungslinien überhaupt 
nicht auftreten können, weil die Momente in den Feldern I und III entgegengesetzte Vor- 
zeichen haben. Schon dies erschwert die Anwendung des Vianelloschen Verfahrens. Es 
zeigt sich aber darüber hinaus, daß die zum 1. Eigenwert gehörige Bieruneslinie eine Null- 
stelle haben kann. Wir erhalten so weitere Beispiele für Nichtanwendbarkeit des Vianello- 
schen Verfahrens. Sei nämlich das Verhältnis z : :q==1 fest gegeben. Bei kleinem a er- 
halten wir dann zunächst zum kleinsten Eigenwert »,< , 5 E eine positive Biegungslinie mit 
endlichen Steigungswinkeln der Tangenten an den beiden Enden. Mit wachsendem a aber 
nimmt die Neigung der einen dieser beiden Tangenten zu Null ab. Für einen Wert a = 1, 
entsteht eine Biegungslinie mit horizontaler Tangente an einem Ende, als ob der Stab dort 
eingespannt wäre. Mit weiter wachsendem a wechselt diese Tangente ihr Vorzeichen und 
entsteht eine Biegungslinie, die auch negative Werte annimmt. Das Umgekehrte findet mit 

/ ge 

der zum zweiten Eigenwert », 1 g: un“ er zugehörigen Biegungslinie statt. Diese hat 
für kleine a eine Nullstelle und das Aussehen einer sin-Linie im Periodenintervall. Mit 
wachsendem a aber nimmt die Tangente an einem Ende — und zwar an dem entzegengesetzten 
wie vorher ab, geht bei einem Werte a, durch Null und wechselt ihr Vorzeichen. Für 
Werte a>a, hat diese Biegungslinie keine Nullstelle mehr. Die zu », und v, gehörigen 
Biegungslinien haben also ihren Typus für genügend großes a gerade vertauscht. Die Abb. > 
v=»,) und 3 (r=r,) zeigen die Übergänge. 


‘), Wir haben hier das in I angekündigte Beispiel, daß der erste Eigenwert mehrfach sein kann. 
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Ich führe den Beweis nur für den Eigenwert »v,. Zur Aufstellung der Gleichung F=0 
lassen sich hier, weil v»(b — a)<z ist, B und € aus (2,21) und (2,22) ausrechnen und in (2,23) 


und (2,24) einsetzen. Man erhält zwei homogene Gleichungen in a und 5 von folgender Gestalt: 


| 
alu +3) ee ee 


| 
ar Pl, „r R). Sara Tr % A 
) ’ 
aus denen sich F=-0 ın der Form ergibt 


LE B+ | v0 97 
al l ' m Zu l bh’ . . . . . . . . . . .4). 


Dabei ist 


N | E me . f Fig B.. a 1 v f = h "eo „o a nr 


3 5 I 3 75 > 
I / (l h)? »? y hb-+a b—-a y ba f R 
R=- 1 ; 5 i te» +7 tl )eot» 5” ...872). 


Wir verlangen horizontale Tangente am linken Ende, d.i. A=0. Setzen wir 
in (2,11) gemäß (2,61) 


| | 
N, =—-+ıN. — 
” 1 do P e a 
ein, SO folgt 
v” a” 
61 (2,51) 
und daraus nach (2,7) 
| v? a? 
— — n 2.82). 
R alter) een 88) 


Dies in (2,71), (2,72) eingesetzt, ergibt unter Einführung der dimensionslosen Größen und 
Bezeichnungen 


OR l h—a z2—(q+1 
(= e 2 Hr vd, z a v > LK : q } r (2.9) 


nach kurzer Rechnung: 


l x? 
(—-VDaxtgy: = 2 +BgHN oe rn er 


7 \” 6 
| a 
(4 +-Dretp= - | 2 + 2 +3q ı)\ EG E -  ° 
und durch Multiplikation dieser transzendenten Gleichungen die algebraische: 
„A gr 
3 I 2?’—-NV)+z(-!—-g+2:.@?—1))—-2e +2 )=0 .. . . RB). 
.)ı) +) 
Gl. (2,91) zeigt, weil 9< , ist, daß notwendig 


EC WHEN GE EREE RE 


sein muß. Wir wählen außerdem 4 so nahe bei 1, daß —2q’+3qg-+1 und daher der 


ne al I+y> 
Koeffizient von .*’ in (2,95) positiv ist, also 4 < ui 


Dann hat (2,93) für jedes positive z 


eine positive Wurzel @. Für z kommen nur Werte>g +1 in Betracht, für diese ist x eine 
ständig wachsende Funktion «(z), die von der Größenordnung Y z unendlich wird®). Man 
denke sich diese Funktion in (2,91) eingetragen. Für 2=g4+1 verschwindet die linke Seite, 


>, Der Beweis des zweiten Teils dieser Behauptung folgt einfach aus (2,93), der erste Teil ist für das Folgende 
unwesentlich, weshalb ich den Beweis weglasse. Die Bestimmung des Intervalls sämtlicher Werte qg, für die eine 
horizontale Tangente auftreten kann, und die Bereehnung'der zugehörigen Werte a und v (d.h. z und x) ist eine inter- 
essante algebraische Aufgabe, deren Lösung ich durehgeführt habe, aber hier unterdrücke. 
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die rechte ist positiv. Man lasse aber jetzt z wachsen, bis 9 den Wert 5 annimmt. Dann 
ist die linke Seite größer als die rechte. Somit gibt es ein z, d. h. ein a, derart, daß die 


zugehörige, dem ersten Eigenwert v,< z entsprechende Biegungslinie am linken Ende eine 


>h 
horizontale Tangente hat. Eine zusätzliche Betrachtung zeigt, daß die Tangente mit wachsen- 
dia A) 


dem a ihr Vorzeichen wechselt, denn 


| 3 
1 ax u 


/ 


re ergibt sich als positiv. Somit ist für Werte g 


unsere Behauptung bewiesen. Für die Abbildungen wurde q = 1,25 gewählt. 


$3. Hinreichende Kriterien für Anwendbarkeit des Vianelloschen Verfahrens. Da, wie die 
Gegenbeispiele des $ 2 zeigen, das Vianellosche Verfahren zu falschen Ergebnissen führen 
kann, wäre es wichtig, genügend allgemeine hinreichende Kriterien für seine Gültigkeit zu 
kennen. Es ist mir trotz vieler Rechnungen nicht gelungen, in dieser Richtung so weit zu 
kommen, wie ich wollte; ich kann aber doch einige Ergebnisse angeben. Sie folgen aus einem 
zuerst von Kellogg*’) bewiesenen Satze, dessen einfachster Fall so lautet: 


(K,).. Wenn der stetige, symmetrische Kern @G(r,s) füralle r,s des 
Intervalles nicht-negativ und für r=s mit Ausnahme der Grenzen 
positiv ist, dann ist der kleinste Eigenwert positiv und die einzige 
zu ihm gehörige Eigenfunktion ist im Innern des Intervalles positiv. 


Kellogg hat diesen Satz verallgemeinert durch Einführung der „assoziierten Kerne“ '®) 


Gi): Fe 


5 FO. EA Rn 
s u. (i (Yu, s,) ... (1 (ru, Sn) 
Er findet: (X). Wenn für 
Er SS Fir u, <S:..<a, 
die sämtlichen Kerne 
ZW N Bi lI<sisn 
S So o Si 


nicht-negativ und im Falle der Gleichheit aller Paare r,s sogar mit 
Ausnahme der Grenzen positiv sind, dann sind die ersten n Eigen- 
werte von G(r,s) positiv, zu Jedem gehört nur eine Eigenfunktion, und 
diese haben der Reihe nach 


0, 41, n—1 

Vorzeichenwechsel in dem Intervall, Ihre linearen Kombinationen 
haben nicht mehr als n—1 Vorzeichenwechsel. Alle anderen Eigen- 
funktionen haben mindestens n Vorzeichenwechsel, 


A. Den Satz (A,) kann man bei dem beiderseits gestützten Balken anwenden. 
Vianellos Verfahren geht hier von einer überall positiven und konkaven Kurve aus, so- 
daß m“ (s) nicht-negativ ist. Ist @ (r, s) nirgends negativ, dann sind auch m” (s), m® (s),... 
nicht-negativ, das gleiche gilt daher auch für die Eigenfunktion, gegen die sie streben ''), Da 
aber laut (A,) bei nieht-negativem und für r-=s positivem G(r,s) die zum 1. Eigenwert 
gehörige Eigenfunktion positiv ist und deshalb wegen der Orthogonalität keine andere nicht- 
negativ sein kann, so führt das Verfahren in der Tat auf den ersten Eigenwert, w. z. b. w. 


»), Am. J. of Math. 38 (1916), S. 1 bis 5, und 40 (1918), S. 145 bis 154. 
10) Ene. Math. Wiss. II 3, S. 1385. 
11) (Genauer gesagt: jede konvergente Teilfolge aus ihnen. 
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Um das Kriterium auszuwerten, berechnen wir @G(r,s), wobei wir zur Vereinfachung 
der Schreibweise E.J = const = 1 annehmen. Es ist 


M(x,)= ae 








und daher, für rs, 


r s I 
NE le y\roaı li N PiDdt+ .\ Pddt.... 082. 
6 r s 


ezeiechnen wir mit P,,P,,P’, Mittelwerte von P in den Intervallen (0, r), (r, s), (s, D), so läßt 
sich @(r,s) in eine der beiden folgenden Formen setzen: 


G(r,s) 1 \rili ip, P)+[(1 NP+ SP 











pas Barth 
| FL Dear ie 2 aualr & 3 u 


und daraus ergibt sich: Wenn der beiderseits gestützte Stab überall gedrückt 
ist (P>0) und die Diehten der Längskräfte p(f) überall gleiches Vorzeichen 
haben (7’(f) nicht zunehmend oder nicht abnehmend), dann liefert Vianellos Verfahren 
die Knieklast. Das Kriterium gilt z. B. im Eulerschen Fälle und in Vianellos Beispiel 1. 
Es ist nieht erfüllt (in unserem Gegenbeispiel und) in Vianellos Beispiel 2, doch läßt sich 
bei diesem zeigen, daß das allgemeine Kriterium @ (r, s) 20 hier besteht. 


B. Falls eines der beiden Enden eingespannt ist, kommt man mit so einfachen 
Betrachtungen nicht aus, weil dann m (s) sicher sein Zeichen wech$elt. Man muß dann Sätze 
(K,) oder (X,) heranziehen. 


Wir beschränken uns auf den Fall EJ=1 (der allgemeine Fall liege sich nach der 
eleichen Methode behandeln). 


Falls das linke Ende eingespannt, das rechte gestützt ist, findet sich 

















f? ” / ge ) 
an n(1- r\=Mı(@) «>21 
Mir, ft) er a (3.5) 
Med (-%) e<t 
er, 93=l1-)lı 25 \rneli jr 
l I’). 21 
0 
| | (34), 
s\» RN r\» ! 
\ | ‚\ (nelı EnKL \ )F\4 nilı ,)a 
+ /'(r,s) 
l 
\Pidt rs 
ir,s)—i, ee 
Eddie ’ZSs 
Weil m (r) die negative 2. Ableitung der an den beiden Enden verschwindenden Aus- 


biegung y(r)ist und die erste Ableitung für r =0 verschwindet, ergibt sich 


I r I 
0 vl HOW) =\dr\m()do-\(l - r)m(r)dr, 


0 0 0 
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und da diese Orthogonalitätsbeziehung für jede Eigenfunktion des Kerns @ (r,s) gilt, so gilt 
sie nach dem bekannten Satze von E. Schmidt auch für den Kern selbst: 

l 

Be u re a (3.5). 


0 
was sich auch leieht durch direkte Rechnung bestätigen läßt. Daher hat 


+) 
+) 


0, = ATi 1 T)+60, 8) 3,6) 
vie, Bi) au } > irrt, 8 i ; Ä j ; ; : ; i (db 
| ı/\ | 

die gleichen Eigenwerte und Eigenfunktionen wie G (r,s) und dazu noch die ihr Vorzeichen 


nieht wechselnde Eigenfunktion I j zum Kigenwert 


Wir setzen jetzt voraus, daß ein positives A sich so groß wählen läßt, daß @G,(r,s) und 


Gi - niecht-negativ sind in der Weise, daß den Bedingungen des Kelloggschen Satzes 
=. 4. 


l p. 
(K,) genügt ıst. 


Die Ausgangskurve des Vianelloschen Verfahrens genügt den Randbedingungen, wir 
haben also, wie oben, 
! 


Erna . zer ra 
ı 





Die weiteren Funktionen m’ (r), m” (r),... erfüllen diese Bedingung gleichfalls infolge 

der Eigenschaft (3,5) des Kernes. Außerdem wird die Ausganeskurve immer so gewählt, 

daß die Krümmung und also m” (r) nur einmal das Vorzeichen wechselt. Wir wollen der 
m” (r) 


Einfachheit halber darüber hinaus voraussetzen, daß mr eine monoton wachsende Funktion 
ist, eine Bedingung, die sich leicht erfüllen läßt und bei den üblichen Zeiehnungen erfüllt 
sein wird. | 


Wir benutzen die folgende, leicht beweisbare Identität'’). Seien 9,(r), 9, (r) beliebige 
stetige Funktionen und A(r,s) ein beliebiger stetiger Kern, dann ist 


I l 

\Air,s)g,(n)dr \Kir,s)go,(n\dr Mi 

0 1) if, fa 

I I .) \ \.er)a A erplr,)) K|,' ja r,dr,. 
r " u u 

\Kir,s)y,s(r)dr \K(r,s,)y,(r)dr iz 


Wir wenden sie an auf 


e ö 
o,(r)=1 Y,lr) = m” (vr), K(r,s)=G,(r, s). 
Ks ıst dann 
y S 
\Kinsoun)dr Ali r): 
0 
I ! ! 
' 3. s\(/ r\ 
\Kinsg.in)ar: 7 1 ı\ N uetndrt\eo, s)m"(r)dr 
0 0 0 
! 
\Gfr,s) m" (r) dr = m“ (s) nach (3,50). 
0 
' mn‘ (r) 
Da nun wegen des monotonen Wachstums von en der Ausdruck 
(o,(r) 9; (r,) — p. (r.) 9, (r,))=( — r,) m” (r,) —- (1 — r,) m‘” (r,) (r,<r,) 
sicher positiv ist, und der positive Charakter von ls) vorausgesetzt wird, so ist auch 
u 
1 S, (1) /a N S;| da u 
m nm (S,) | jjm (S,) (Ss, <S,), 


12) Siehe Kellogg: Am. J. Math. 40 (1915), S. 151. 
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die linke Seite der obigen Integralidentität, sicher positiv. Ebenso weiter die entsprechenden 


Ausdrücke, die mit m“, m®’,... gebildet sind. Hieraus ist aber leicht ersichtlich, daß 
alle m’ (vr), m” (r),... nur einen Vorzeichenwechsel haben. Das gleiche gilt also von m (r) 


als Grenzfunktion '?). 


Nun haben nach dem Kelloggschen Satz (K,) die beiden ersten Eigenfunktionen des 
Kernes @,(r,s) keinen bzw. einen Vorzeichenwechsel, während alle weiteren Eizenfunktionen 
mindestens zwei Vorzeichenweechsel haben. Deshalb gehört das nach dem Vıianelloschen 
Verfahren gebildete m (r) sicher zu dem 1. Eigenwert von @G(r,s). Also: Falls bei dem 
linksseitig eingespannten, reehtsseitiggestützten Balken sich ein positives 

r.r, ' 
Aso groß wählen läßt, daß @, (r,s) und ef) nicht-negativ entsprechend 
ua Be, 
den Bedinzeuneen des Kelloreschen Satzessind, dannliefert das Vianellosche 
Verfahren in der üblichen Form den Knick wert. 

Die Auswertung dieses Kriteriums wird überraschend einfach. Zunächst ist sofort zu 
übersehen, daß @,(r,s) durch genügend große Wahl von A positiv gemacht werden kann, 
außer wenn rl oder s= 1! ist, in welchem Falle es verschwindet. Betrachten wir weiter 


;> £ ra) 3 )( =) r,s 
Y’.W, N ‚\i ru ‚\ mir = rG(r,,s,) 
else # R a 


A r 1 ı )+teo,s) A ; hi Al 7) HG (r,, s,) 


und ordnen nach Potenzen von A, so fällt das Glied mit A? heraus, der Koeffizient von A aber 


. 
A > . + ee y . . \ 
erhält bis auf den Faktor ] nach einigen Vereinfachungen die Form 


Cr. 18.8.) hi hi Tr s)+[! al res) 


\ 


Di? rs [i al rir.s.). 


/ 


Betreffs der gegenseitigen Lage der Punkte r,,r,;, s,, s, können bei Berücksichtigung der 
Bedingungen r,< r,.s, <s, sechs Fälle eintreten, von denen aber wegen der Symmetrie des 
Kernes G (r.s) nur 3 betrachtet werden müssen: 


ll. ,<n. Su, <a; Fhr,e) "r.s)=iPlhdt: Pr,s)=ftw,;#) rat, | 
0 
Che, 6,6 | (r, zu 7 ‚\rnat E v\ronan 
| £ 8, | 
inne, sl alz P,), 


wo P’, den Mittelwert im Intervall (s,, s,), 7, denjenigen im Intervall (s,, !) bezeichnet. 


! I l 
>» 1.5.54 AO RHE IE | 7 75 Fer.s)=\P(hdt; "(r,s)=I{r,s)=|\Pfltdt, | 
I I ra 
" | | s,\\ . S;\( , | r,\ Ss ) 
en) le li ‚\rmaı 1 ‚\rwnan+(i ‚I ‚\\romaı 
No St Sr 
| NS, 1 S | 
nn) 5) 1-7) Py)+|i Jh Ten —-sB., | 


wo P’, und P, die gleiche Bedeutung wie oben haben und P, der Mittelwert in (s,,r,) Ist. 


> 
- 


13) Siehe Fußnote #1), 
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3. , au <a. a0 "er,s.)=\PYdi; Fe, )=) Pd; Ir, )=Flr„e)=|Plidi, 
s1 S2 ro 
I I I 
Sn S, S, 1 > N, > N Sy N, N 
RC RE u 6 ‚)\rwaı 1 ‚)\rwar+ \rnaı 
r2 So ra 
P, / S, Fr 1 
{1 =)(i \\rmaı 


(s, Sl = =(P"—P))+ 7 Pi" +l1-7)P, 





wo P, und P,” die Mittelwerte in den Intervallen (s,, r,) und (r,,) bezeichnen, während ?, die 
eleiche Bedeutung wie oben hat. 


Daraus schließen wir: Wenn die Pnicht-negativsindundvomeingespannten 
zum gestützten Ende nicht abnehmen, dann ist das Vianellosche Verfahren 
sicher anwendbar. 


Die Koeffizienten C,, C,,C, sind nämlich dann im Kelloggschen Sinne nicht-negativ. 

| r,T, ee 2 

Berechnen wir außerdem das konstante Glied von 6 “ eine Rechnung, die ich hier 

FrX "u 

unterdrücke —, so zeigt sich, daß dieses, abgesehen von dem nicht-negativen Summanden 
I\r,s)I'(r,, s,) - Ir, s,)/'(r,,s,), in jedem der drei Fälle eine lineare Kombination der 
S, r, 
5 | EEE ist wie das entsprechende €. 
Wenn daher in- einem © die Koeffizienten dieser Potenzprodukte von Null verschieden sind, 


gleichen Produkte von Differenzen vr, r,, 8; —8,1- 


. .. » Y ® Y r N", .. » \ 
dann läßt sich durch genügend große Wahl von A sicher G, ” i überall außer an den Grenzen 
u ar ., 

positiv machen. Es bleibt also nur noch zu beweisen, daß, dieses Ergebnis auch bestehen 
bleibt, wenn der Koeffizient eines Potenzprodukts Null ıst. Hier zeigt nun die Betrachtung 
des konstanten Gliedes, daß nur das Verschwinden der Differenzen P,— P, oder P,’ —- P,' 
Schwierigkeiten macht. Diese Differenzen verschwinden dann und nur dann, wenn Pit) in 
dem Intervall (s,!) bzw. (s,,l) konstant ist. Unter diesen Umständen ist aber, wie die 
Rechnung zeigt, das konstante Glied von Gl „ selbst positiv, so daß auch dann die Be- 

« ı ?2 


dingungen des Kelloggschen Satzes für @,(r,s) sichergestellt sind '”). 


C. Zum Schluß sollen noch drei leichter und vollständiger behandelbare Randbedingungen 
besprochen werden, wobei immer EJ=1 gewählt wird. 


Falls das linke Ende eingespannt, das rechte frei ist, wird 


(." t u A 


M (x, f) I r ; 1 E 
Daraus folgt 
! 
Gir,s)\ Pat (rs) 


14) Der gewonnene Satz umschließt insbesondere den Eulerschen Fall. Leider ist er nieht anwendbar auf 
Vianellos einschlägiges Beispiel 5, wo nämlich die P vom eingespannten zum gestützten Ende abnehmen. Es ist mir 
auch nieht gelungen, durch Übergang zu iterierten Kernen (die schärfere Kriterien liefern müssen) die Richtigkeit des 
Verfahrens bei dem Beispiel 5 zu erweisen. 

Im Falle des beiderseits eingespannten Balkens scheinen sieh aus dem Kelloggschen Satz keine hin- 
reichenden Kriterien ableiten zu lassen. Die Momentenlinien m(r) haben hier mindestens 2 Nullstellen, so daß ein 
Kriterium (Ka) anzuwenden wäre. Formal bietet sieh hierzu auch ein Weg, indem man den Kern 


Gy) = 14 B[I- )6- tens 


heranzieht. Dieser kann aber bei den gewöhnlichen Formen der Biegungslinie — insbesondere im Eulersehen Falle 

das Kriterium (Ka) nieht erfüllen. Denn im Eulerschen Falle ist die zu der einfachsten Biegungslinie gehörige Momenten- 
linie eine Kosinuslinie im Periodenintervall und wird von gewissen Geraden in diesem Intervall in 3 Punkten geschnitten, 
entgegen dem Kellogzgschen Satz. 
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und ist für Pit)-20 nieht-negativ. Das Vianellosche Verfahren, ausgehend von einer 
überall konkaven Biegungslinie, führt also zum Ziel’). (Vianellos Beispiele 3 und 4.) 

Ist das linke Ende gestützt, das rechte in einer senkrecht beweglichen Backe eingespannt 
(Falld, $ 1), dann ıst 


ak Ze] 
und also 
G(r,s)=\P(f dt (rss). 


\ueh in diesem Falle ıst also beı überall gedrückten Stäben (Pt) =0) das Vianellosche 
Verfahren, auszehend von einer überall konkaven Biegungslinie, gerechtfertigt. Bekanntlich 
bestehen diese Randbedingungen für die linke Hälfte eines beiderseits gestützten, durch Druck - 
kräfte symmetrisch zur Mitte beanspruchten Stabes von der Länge 21, falls der Stab symmetrisch 
auskniekt. Die Knieklast eines gedrückten, beiderseits gestützten, symmetrisch beanspruchten 
Stabes entspricht also, entweder der rein konkaven symmetrischen Ausbiegung oder einer 
antisvmmetrischen Ausbiegung. Bei der antisymmetrischen Ausbiegung erfüllt die Biegungs- 
linie des halben Stabes an beiden Enden die Randbedingung der einfachen Stützunge. Nehmen 
also die Druckkräfte vom linken Ende zur Mitte dauernd zu oder ab, so hat nach «den Ergeb- 
nissen von Abschnitt A. diejenige antisymmetrische Biegungslinie den kleinsten Eigenwert, 
die links von der Mitte rein konkav verläuft. Im ganzen also zeigt sich: Wenn ein beider- 
seits gestützter Stab von Druckkräften beansprucht ist, die zur Mitte 
symmetrisch sind und von den Enden zur Mitte monoton zu- oder abnehmen, 
dann ist die der Knieklast entsprechende Ausbiegung entweder die ein- 
[achste symmetrische oder die einfachste antisymmetrische. Dieser Satz gibt 
für den genannten, praktisch wohl wichtigsten Belastungsfall eine willkommene Aufklärung 
und Ereänzunge zu den Betrachtungen des $ 2 und zeigt, daß man in diesem Fall durch An- 
wendung des Vianelloschen Verfahrens auf zwei Auszangslinien, nämlich auf die linken 
Hälften der symmetrischen und der antisymmetrischen Biegungslinie, die Knicklast sicher 
bestimmen kann. 


In ähnlicher Weise läßt sich auch der beiderseits eingespannte, symmetrisch 
beanspruchte Balken behandeln, dem wir wieder die Länge 21 geben wollen. Ist die 
Biegungslinie symmetrisch, so besteht für die linke Balkenhälfte bei O0 die Randbedingung b, 
bei ! die Randbedingung d. Damit ergibt sich 








! 
Mia, tı 
z | 
21 _ 
und 
! 
I; I; Ic 
ns Pinrat 7 \Pintat 7 \Pitat+\Pibat (re). 
Krzänzen wir den Kern zu G, (r,s)= 1+ G(r, s), so erfüllt dieser das Kelloggsche Kriterium (A, 


für alle nieht-negativen P. Der Kniekwert des gedrückten, beiderseits eiın- 
gespannten, symmetrisch beanspruchten und symmetrisch ausknickenden 
Balkens läßt sieh daher nach Vianello ermitteln, wenn wir voneiner ständig 
wachsenden Momentenkurve auszehen'*). 


Ist aber die Biegungslinie antisyvmmetrisch, so erfüllt das rechte Ende der linken Balken- 
hälfte die Bedineuneen eines eestützten Endes. und wir sind also auf den oben behandelten 


KallB. zurückgeführt, der unter Umständen wenn nämlich an dem symmetrischen Balken 
die Längskräfte von den Enden nach der Mitte zunehmen -— eine Entscheidung durch Vianellos 
Verfahren zuläßt. 628 

15) In diesem und den folgenden Fällen ist zwar die Kelloggsche Ungleichung A (r, r) > 0 nieht immer erfüllt, 
dann nämlich nieht, wenn Pif) auf gewissen Streeken identisch verschwindet. Diese Ausnahmen sind aber einfacher 


\rt und beeinträchtigen die Gültigkeit der Resultate nieht, wie man dureh leichte Zusatzbetrachtungen erweisen kann. 


x wie unter B. zu führen sind, lasse ich weg. 


1 Die Beweise, die ganz analı 
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KLEINE MITTEILUNGEN 


Näherungsformel für das Auskühlen 
einer Platte mit konstanter Anfangs- 
temperatur. 


Verwendete Zeichen: 
2 X [m] Dicke der Platte in Metern. 


ı.h !.Grad !! Wärmeleitzeschwindigkeit 


(les Plattenstoffes. 


/)\keal-m 


kz-m#| dessen spez. Gewicht. 


ce |keal-kg"'-Grad !| dessen spez. Wärme. 


”. Grad) Anfangstemperatur der 


Platte. 
#,|Grad| Umgebungstemperatur. 
t\h]| Zeit in Stunden. 


(). \keall Anfänglicher Wärmein- 
halt der Platte, über der 
Umgebungstemperatur ge- 
INIeSSen. 

Platte 


Wärmeinhalt der 


zur Zeit f. 


(),;,  keal 


#,, \Grad| Wandtemperatur. 


alkeal-m ?-h '-Grad '| Wärmeüberzangszahl von 
der Plattenoberfläche auf 
die Umzebung. 
7. a 
Abkürzungen: a ch 
C+Y v. 


Hier soll eine Näherunesformel für den Fall des 
Auskühlens bzw. ÄAnwärmens einer unendlich aus 
eedlehnten. ebenen Platte endlicher Dieke aus dem 
Stationärzustand mit überall gleicher Temperatur %. 
dureh Absenken bzw. Erhöhen der Umzebungs- 
temperatur auf den konstanten Wert %, angegeben 
und ihre Herleitung kurz erläutert werden. Zweck 
dieser Formel ist es. diese häufige Aufgabe der 
Bereehnune mittels Reehenschieber und e X-Tafel. 
wie sie in den meisten Handbüchern zu finden ist. 
zueänglich zu machen. Die zu diesem Zweck von 
Prof. Dr. Gröber!)errechneten und vonG.Pösch]?) 
erweiterten Kurvenblätter haben: den Nachteil. 
dal man den Verlauf der zu bereehnenden Größen 
(Wärmeinhalt und  Oberflächentemperatur) nur 
punktweise erhält. wobei noch gewisse Bereiche der 


Variablen A X stiefmütterlich behandelt sind. da 
es eben unmöselieh ist. den zesamten Bereich von 
Ih X 0 bis AX x auf einem Kurvenblatt dar 


zustellen. 


Nach der Fourierschen Theorie ergibt sich der 
Wärmeinhalt der Platte zur Zeit ft: 


/ 4 n 
a rs 7 
sin?” O7, Ok 
(J, Je \ 2 - . € = 
O1,” d,- sin O7; 608 O7, 
7 l 


wobei die Of; 
Gleiehun®e 


die Wurzeln der transzendenten 


cotg 0: 


!,Gröber-Erk: Die Grundgesetze der 
tragung, Verlag Springer, Berlin. 


?) ZAMM Ba. 12, 19532, Heft 5, S. 


Wärmeüber 


2s0, VDI-Verlag, Berlin. 


die sich aus den Anfangsbedingungen ergibt, be- 
deuten. ‚Jede der wellenförmigen Temperaturver 


teilungen, aus denen die anfängliche Temperatur 
einen 


verteilune aufgebaut ist, liefert Anteil zum 


q [An N a, 
700 L. | 700 3 
700 





504 





\ 
s% 
N Dyiy y 
ST Penfempergtyp 
_— Un 
R N re w 


Abb. 1. Verlauf der Oberflächentemperatur und des 

Wärmeinhaltes bei AX=10. (Die gestrichelten Linien 

eeben den tatsächlichen Verlauf an, während die voll 

ausgezogenen Linien den dureh die Formel errechneten 
' Verlauf angeben.) 
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Abb. 2. Verlauf des Wärmeinhaltes bei IX n, 


Wärmeinhalt der Platte. der um so rascher ver- 
schwindet, je größer ö, ist, das heißt. je kompli 
zierter die entsprechende Einzelverteilung ist. 


Diese unendliche Reihe zeigt für eroße Werte 
von AX und für kleine Werte von f eine sehr 
schlechte Konvergenz, jedoch liefert der erste 
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Summand auch bei unzünstiester Konvergenz, das 
m 
- Ä 
ist bei AX vo und Z=0, noch immer —,— 81°/o 
des Wärmeinhaltes der Platte. Der Anteil des 
ersten Summanden an Q. ist eine Funktion von Öj,, 
also von AX. wird mit G (Grundverteilung) be- 
zeichnet und aus einer unten anzerebenen Formel 
mit hinlänzlicher Genauirkeit berechnet. 


Der Exponent der e-Potenz, nach welcher die 
Grundverteilunzabklingt. ist ebenfallseine Funktion 
von A X und außerdem proportional der Fourier- 


dl» a , = 
schen Kennzahl v: Sämtliche anderen Einzelver- 


teiluneen wurden zusammengefaßt und die zeitliche 
inderun® dieses Restes der Fourierschen Kenn- 
zahl proportional und als Funktion von AX derart 
festeeleet,. daß der Fehler möglichst klein bleibt 
(Abb. 2). Die Formel für den Wärmeinhalt nimmt 
dann die Gestalt an: 


( 
too 
( 


zu 


((.e hit 4- (100 G).e—hket, 


Im Interesse der Einfachheit der Formeln für 
G, k, und A, erwies sich eine Teilung des Inter- 
valles als notwendig. 


Der praktisch am meisten in Frage kommende 
Bereich liert zwischen AX=0 und AX 10. da 
bei dieken Wänden aus schlechtleitendem Material 
ein guter Wärmeübergang selten ist. In diesem 
Bereiche liefert das anzerebene Näherungsverfahren 
unbedingt genauere Werte als die Ablesung aus 
Kurvenblättern, da der maximale Fehler nur auf 
den Anfang des Vorganzes beschränkt bleibt und 
stets kleiner ist als 1°/, des ursprünglichen Wärme- 
inhaltes der Platte. (1°/, erst bei A X = 1.) 


Bei AX xo- tritt aus Konvergenzeründen der 
überhaupt größte mögliche Fehler auf. Dieser 
Fehler beträgt 3°/, von Qe. ist also für praktische 
Anwendungen belanglos, zumal er auf den Beginn 
des Vorganges beschränkt bleibt und rasch kleiner 
wird (Abb. 2). Denn bei praktischen Anwendungen 
liegt die Hauptschwierigkeit in der richtigen An- 
eabe der Stoffwerte und vor allem der Wärme- 
übergangszahl a, wobei bedeutend größere Un- 
renauiekeiten auftreten. 


Da die Näherungsformel für Q; beim Differen- 
zieren einen einfachen Ausdruck ergibt, der dieselben 
Größen wie die Formel selbst enthält, läßt sich für 
die Oberflächentemperatur ein einfacher Ausdruck 
herleiten, der allerdings keine so genauen Werte 
liefert wie die Formel für Q,, aber den tatsäch- 
lichen Verlauf der Oberflächentemperatur sehr 
charakteristisch und mit praktisch genügender 
Genauirkeit wiedergibt, indem er das rasche 
Abfallen der Oberflächentemperatur bei Beginn 
des Vorganges dureh eine Unstetiekeit ersetzt 
(Abh. 1). Ein Teil der Wand von der Dicke 2X, 
der Breite Y und der Höhe Z mit dem Wärmeinhalt 
().=0:y-® +2 X Y-Z gibt seine Eirenwärme durch 
eine Oberfläche von der Größe F=?2-.Y-Z an die 
Umgebung ab. Zur Zeit ? fließt durch diese Ober- 
fläche ein Wärmestrom 


dd 
- r_[k,:G: kıtı+ k.-(100-G)- hut], 
TE Tee Ri er Ä 


Hieraus errechnet sich die Oberflächentemperatur: 


era ec F hit t RR, (100 (i)-e lvo ur 


F,n X 
oder: 100 “ 


U ah 


[k-Gre—kıt+ k,-(100- G)-e— kat), 


Wie aus Abb. I ersichtlich ist, liefert diese Formel 
erst brauchbare Ergebnisse, wenn der Wärmeinhalt 
der Platte auf 96°/, seines Anfangswertes gesunken 
ist. 


Die Näherungsformel und ihre Anwen- 
dung. Bei der Stellunz der Aufgabe müssen be- 
kannt sein: Die Stoffwerte /. y und ce der Platte, 
die Dieke der Platte 2X, die Anfangstemperatur %.., 
die Umgebungstemperatur %, und die Wärmeüber- 
eaneszahl «@. 


i 2 j a 
Hieraus berechnet man «@ IX 
C+Yy 


/ 


-X und 


- 


dlann je nach der Größe von AX: 
für AX kleiner als 1'2: 


(= 100 13 (h X), 


| 


h X (hi X)? a 


MBH ET “.; 


hi, — 14 


PETE TER. 
1 ( 1 hX | 4° 
157 1” 
k 11% 
RN 
a | 500 
e, 50 
er 3 a 3 = r| 


Mit den so erhaltenen Werten geht man in die 
Näherungsformeln: 


( 
10, — (re —kıt 4- (100 


ce 


G)-e—hzt 


(in °/, von Q.), 


Pa hi 


A. G-e—kıt 4 k,-(100 
Öe ah” ab Zn ne 


G)e—hkzt] 


(in °/o von u, gilt erst ab Q, =. 96°/, von Q.). 


Diese Formeln lassen sich einfach zeichnerisch 
darstellen. Sie gelten natürlich auch für den Fall 
der einseitig vollkommen isolierten Platte, wobei 
dann die Plattendieke mit X zu bezeichnen ist. 


Prag. F. Nistler. 690 
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BUCHBESPRECHUNGEN 


Professor JULIUS RATZERSDORFER, 


Dr.-Ing. 


Die Kniekfestirkeit von Stäben und 
Stabwerken. IX-+2321 S, m. 151 Abb. Wien 
1036, Verlag Julius Springer. Preis geb. 28,80 M. 


Seit der .„Knickfestigkeit“ von R. Mayer (er- 
schienen Berlin 1921) erschien in deutscher Sprache 
keine zusammenfassende Behandlung der für den 
Ineenieur immer wichtirer werdenden Stabilitäts- 
fraeen. Da seither in unseren Kenntnissen über 
Stabilitätsgrenzen und Labilitätserscheinungen 
elastischer Tragwerke zahlreiche und erhebliche 
Fortschritte erzielt wurden, wurde das Fehlen einer 
ınodernen Darstellun®e dieser Fragen immer mehr 
als störend empfunden, Das Buch von Ratzersdorfer 
füllt einen großen Teil dieser Lücke aus, wobei es 
trifft, daß es 


sich  elücklieh sich mit zwei 
andern gleichzeitig erschienenen Büchern über 
Kniekune auf vorteilhafte Weise ergänzt. Das 


buch von Mayer behält übrigens wegen der ab- 
weiehenden Zielsetzung besonders für viele prak- 
tische Fragestellungen seinen Wert unvermindert bei. 


Die Grundlage der Ausführungen des Verf. bildet 
die ausführliche Besprechune der Knickfestigkeit 
des zeraden Stabes mit konstantem Querschnitt für 
verschiedene Fälle der Randlagerung und der B»- 
lastune. Es wird sowohl die elastische wie die 
unelastische Kniekunge betrachtet. Hinsichtlich der 
Abhängiekeit der Kniekspannungen bzw. des Knick- 
moduls von der Querschnittsform und den Baustoff- 


eieenschaften im unelastischen Bereich wäre eine 
eingehendere Darstellung erwünscht gewesen. Das 


neuere Schrifttum ist teilweise nur ungleichmäßig 
berücksichtiet worden. 

Der nächste Abschnitt behandelt einige Fälle des 
eedrückten geraden Stabes mit gleichzeitiger Bie- 
eungsbeanspruchung. Hier sind jeweils auch die 
in diesen Fällen besonders wichtigen zweiten Nähe- 
runeen für die Ausbierungen angegeben. Sodann 
folet ein Abschnitt über den geraden Stab mit stetig 


veränderlichen Drucekspannungen und stetig ver- 
änderlichem Stabquersehnitt. 
Der eirentliche Wert des Buches liegt in den 


Abschnitten über die Stabilität des geraden Stabes 
bei kontinuierlicher Lagerung, die Kniekfestigkeit 
von räumlichen und ebenen Stabwerken sowie das 
Problem «der Seitensteifirkeit. Insbesondere das 
Gebiet der Fachwerkknickung, zu dem vom Ver- 
fasser wertvolle eirene Beiträge gegeben wurden, 
hab hier eine schöne abgerundete Darstellung 
[unden. 


OP- 
gt 


Die Behandlung der Knickunz von Rahmen und 
l"achwerken beschränkt sieh in der Hauptsache auf 
len elastischen Bereich. bei dem großen prak- 
tischen Interesse, das die unelastische Kniekung 
hat, wäre bei einer Neuauflage eine Ergänzung nach 
(dieser Riehtung erwünscht. In diesem Zusammen- 
hanre erschiene auch eine einzehendere Behand- 
lune der z. B. für schwierigere Fälle zur Verfügung 
stehenden Annäherungsmethoden zweckmäßig. 

Den Schluß Buches bildet ein kurzer Ab- 
schnitt über einige Sonderfälle der Stabilität von 
borenträrern in ihrer Ebene, die einer einfachen 


des 


mathematischen Behandlung zueänelich sind. Die 
 B. bei den weitgespannten Bogenbrücken auf- 


tretenden Fragen sind jedoch nicht berücksichtigt. 

Die oben ausgesprochenen Wünsche vermindern 
nicht den großen Wert des Buches, das sieh bald 
seinen Platz besonders für das wichtige Gebiet der 


Kniekune von Fachwerken, Rahmen und Stab- 
werken sichern wird. 
Hannover. F. Schleicher. 688 


CARL STORMER, Resultats des ealeuls 
numeriques des trajectoires de cor- 
puscules @elestriques dans le champ 


d’un aimant &l@ementaire, IV. Fais- 
ceaux de trajectoires de ]’infini, 
avec asymptotes paralleles A l’axe 


des X 95 S. V. Faisceaux 
toires avec asymptotes 
l’axe magnetiques et 
axe 9M. 
Dybwad. 


de trajee- 
paralleles ä 
normales ä cet 
Oslo 1936, 1 Kommisjon Hos Jacob 


In den Jahren 1913 und 1914 hat Stormer 
eine Reihe von durch numerische Integration 
eines Systems von Differentialgleiehungen zweiter 
Ordnung gewonnenen Resultaten über den Verlauf 


von Elektronenbahnen im Felde eines Elementar- 
magneten veröffentlicht. Diese Rechnungen hat 


er im Jahre 1930 zusammen mit einer Reihe jüngerer 
Fachgenossen wieder aufgenommen. Dabei hat eı 
die Gleichungen aber nicht in der früheren Form 


benutzt, sondern hat verschiedene für die nume- 
rische Integration praktische Umformungen vor- 


genommen. Der Zweck dieser neuen Rechnungen 
ist einmal, das Studium der periodischen Elektronen: 
bahnen zu vertiefen, ferner soll eine neue Reihe von 
Bahnen «durch den Anfangspunkt-gewonnen werden 
und endlich soll eine Schar von Bahnen berechnet 
werden, die aus dem Unendlichen kommen und die 
für die Untersuchung der kosmischen Strahlung von 
Wichtigkeit sind. Die beiden vorliegenden Mittei- 
lungen bringen das durch die außerordentlich mühe- 


vollen numerischen Integrationen bis jetzt ge- 
wonnene Zahlenmaterial Willers. 69 


Dr. GERHARD HERZBERG, V’rofessor für Physik 
an der T. H. Darmstadt, z. Zt. Carnegie Gastpro. 
fessor für Physik an der Universität Saskatsche- 
wan, Kanada. Atomspektren und Atom- 
struktur. Eine Einführung für Chemiker, Phy- 
siker und Physiko-Chemiker. — Wissenschaftliche 
Forschungsberichte, Naturwissenschaftliche Reihe, 
herausgegeben von Dr. R. E. Liesegang, Frankfurt 


a. M. Band 37. XVI + 188 S. mit 79 Abb. und 21 
Tabellen. Dresden und Leipzig 1936, Verlag Th. 
Steinkopff. Preis geb. 14 M. 


Die vorliegende Einführung ist für den Lernen- 
den, vor allem den Studierenden, geschrieben. Die 
Theorie der Spektren wird in den allgemeinen 
lsehrbüchern der Physik meist kurz behandelt. da 


ihre Grundlagen relativ einfach sind und im 
wesentlichen in der Anwendung der Quanten- 
theorie auf freie Atome bestehen. Andererseits 


ist eine eingehendere Kenntnis der Tatsachen und 
Gesetzmäßigkeiten für viele Nachbargebiete in 
Physik und Chemie unerläßlich und das Erscheinen 
eines Buches, das auch für den Chemiker ver- 
ständlich ist, daher sehr zu begrüßen. — Es werden 
zunächst die Bohrsche Theorie und die Grund- 
tatsachen der Quantenmechanik in Verbindung mit 
den Einelektronenspektren behandelt. Im 2. Kap. 
folgen Mehrelektronenspektren, Elektronenspin, Zee- 
man- und Starkeffekt, usw. Dann bringt Kap. 3 
mit dem Pauliprinzip den Zusammenhang von Spek- 
tren und periodischen System, Kap. 4 Auswahl- und 
Intensitätsregeln, Seriengrenzen und Intervallregeln 
für Multipletts, Kap. 5 Hyperfeinstruktur. Das 
letzte Kap. enthält Tabellen, Para- und Diamagnetis- 
mus und Anwendungen in der Chemie (Valenz- 
beeriff usw.). 
Potsdam. 


Theodor Sehmidt. 696 
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948 Nachrichten Math. und Mech. 
Ferner sind bei der Schriftleitunz folgende Dr. ROBERT SAUER, a. 0. Prof. a. d. Techn. 

Bücher eingeranzen (ausführliche Besprechung Hochseh. Aachen. Projektive Linienrzeo- 

hleibt vorbehalten): metrie. (Göschens Lehrbücherei, I. Gruppe, reine 


KURT MEUSEMANN, Kapitänleutnant. und Dr. 
KURT STANGE, Studienrat a, d. Marinescehule in 
Flensbure-Mürwick. Mathematik und Ma- 
rine. (Mathem.-Physik. Bibl. Reihe I. Nr. 89.) 48 S. 
m. 37 Abb. i. Text. Berlin u. Leipzig 19936 
B. G. Teubner, Preis kart. 1.20 M. 


Dr. H. WEINREICH, Oberstudiendirektor am 
Schillerrealeymnasium in Stettin. Physikali 
he Beobachtunes- und Denkaufg: 
\lltags. (Mathem.-Physik. Bibliothek 
Nr. 90.) 52 S. m. 5 Fig. i. Text. 
1936. Verlae B. G. Teubner, 


5% 
ben (des 
Reihe 1. 
u, Berlin 
1.20 M. 


H. NEUBER, 
‚rundlaeen 
recehnune. VI 


Preis kart. 


Kerbspannungslehre. 
für genaue Spannungs 


I60 SS. m. 106 Abb. i. Text 


und auf einer Tafel. Berlin 1937. Verlage Julius 
Springer. Preis brosch. 15 M. 
EDMUND LANDAU, Über einiee neuere 


Fortschritte der additiven Zahlen 


theorie. (Cambridlee Traets in Mathematics anıl 
Math« matieal Ph St 5. NT. 35. 34 N, London 037 
Cambridge University Press. Preis 6 sh. 


a.d. Techn. 
(eometrie. 


Dr.-Ine, habil. ULRICH GRAF, Dozent 
Hochsech. Berlin. Darstellende 


Verla 


Leipzig 


und angewandte Mathematik, Bd. 23.) 194 S. m. 
36 Abb. Berlin u. Leipzie 1937. Verlae Walter de 
Gruyter. Preis geb. 9 M. 


KURT MICHEL, Vom Flohelas zum Elek 
tronenmikroskop. (Deutsches Museum, Ab 
handlungen und Berichte. 9. Jahre.. Heft 1.) I + 
37 >. Berlin 1937. VDI-Verlae G. m. b. H. Preis 
M. 


0.90 

JOSEPH MILLER THOMAS, Professor of Mathe- 
maties, Duke University, Differential Sy 
ste Ameriean Mathematical Society  Collo 
Publieations 


Vol. XXL) IX 119 S. New 
1937. publ. N 


bv American Mathematical So- 
Preis 2 $. 


ms. 
quium 
York 


eletv. 


Helsinki vom 21. bis 
25. Juli 1936 abgehaltenen neunten Tagung der 
Baltischen Geodätischen Kommission. Reiligiert 
von Generalsekretär ILMARI BONSDORFF. 294. 
Helsinki 1957, Osakeyhtiö Weilin & Göös Aktiebolae. 


Verhandlungen der in 


Jahrbuch 1936 der Lilienthal-Gesellschaft für Luft- 
fahrt. 647 S. m. 441 Abb. München und Berlin 1937. 


Kommissionsverlae G. Oldenbourz. Preis zeb. 2OM. 


AIR D IND GE: 


TIBOR RADO, Subharmonie Funetions. 
yreebnisse der Mathematik und ihrer Grenzgebiete, 
hrse. v. d. Schriftleitune (des „Zentralblatt für 
Mathematik" 5. Bd. 1.) V+56 8. Berlin 1997. 
\erla@ Julius Springer. Preis brosch. 6,60 M. 


NACHRICHTEN 


174 S. m. 281 Abb. Leipzig 1937, Verlag Quelle X& 
Mever. Preis geb. 4 M. 
Dr. A. Kupper f. 
Am 12. Juni starb im Alter von 32 Jahren in 


\usübung® seines Dienstes bei der Erprobung eines 


Sereltlugzeuges den Fliegertod Dr. phil. August 
Kupper. Abteilungsleiter im Institut für Flug 
mecehanik der Deutschen Versuchsanstalt für Luft 


fahrt. Die DVL. und mit ihr die technische Wissen 
sehaft. verliert in ihm einen der seltenen Männer, 
lie hohes technisches Können mit strenzem wissen 
sehaftliehem Denken in sieh vereinizen — befähigt 


auf beiden Gebieten selbständige Arbeit zu leisten, 
und vor allem AÄnrerungen aus (lem einen Gebiete 
lem anderen fruchtbar zu machen. 

Geboren wurde X. Kupper am 14. März 1905 ın 


P’asine bei München. Herbst 1923 berann er (das 
Studium des Maschinenbaues und der technischen 
’hvsik, 1028 promovierte er an der Münchener 
Universität mit einer Arbeit aus der mathe 
matischen Physik. Schon seit 1924 stand er in 
naher Beziehung zu der ‚Junkers-Flugzeugwerke- 
\. G., der er dann von 1930 bis 1932 angehörte. 
Nach zwei Jahren selbständiger konstruktiver unıd 


versuchsteehnischer Tätirkeit trat er Anfang 1934 


in «lie DVL ein. wo er (seibst praktischer Flieger 
als Leiter an der gleichzeitir versuchstechnisch- 
fliegerisch und theoretisch tätiezen Abteilung für 


"luemeechanik ein weiter, seinen vielseitiren Fähie 
keiten aneemessenes Betätieungesfell hatte. Be 
sondere Vertdienste erwarb er sieh in dieser Stel- 
lung um die Förderune der Zusammenarbeit zwischen 


ler DVL und hervorraeenden Hochschulwissen 
schaftlern war es (loch seiner lebhaften Initiative 
mit zu danken. daß z. B. Prof. Trefftz die Füh- 


lunz mit der DVL aufnahm, ja sich sogar entschloS, 


seine eanze Kraft dem in der DVL bearbeiteten 
P’roblemkreise zu widmen ein Plan. den der 


Tod sehmerzlieh vereitelte. 
Für den Textteilverantworilich: Professor Dr. Fr. A.W illers, Dresden-A.20, Dorotheenstr. 12. 
Copyrieht 1937 by VDI-Verlag G.m.b.H., Berlin NW 





Seine Freunde haben in dem Freunde einen Men- 
sehen verloren, «der. bei aller Leistung als In- 
eenleur umd Wissenschaftler. die Kraft hatte, auch 
in unserer Zeit des groben Hastens ein ganzer 


\iensch zu sein. K. Marruerre 732 


Geseilschait für angewandte Mathematik und 
Mechanik. 


Die «diesjährige wissenschaftliche Tagung der 
Gesellschaft findet vor dem auf die Zeit vom 19. bis 
24. September in Bad Kreuznach anberaumten 


13. Deutschen Mathematiker- und Physikertag eben- 


falls in Bad Kreuznach am Freitag, den 17. und 
Sonnabend, den 18. September 1957 statt. Die 
u PER . a Ten inf an IC Qoantamhar ahon 
ilatapdi VebDSeribiiiibtaih,. \vııtl Alil iO, ar pt 11 U 1 Aruı—N 

halten werden. Es wird webeten. Vortrags 


anmeldungen bis spätestens zum 1. August d, J. an 


den Geschäftsführer der Gesellschaft für ange- 
wandte Mathematik und Mechanik, Prof. Dr. ©. 
Weber, Drescden-A. 16, Hindenburzufer 15. gelangen 


zu lassen. 

Persönliches. 
Das Preußische Staatsministerium hat die von der 
Preußischen Akademie der Wissenschaften in Berlin 
vollzorene Wahl des ordentlichen Professors an «(ler 
Universität Göttinzen Dr. Ludwie Prandt] zum 


auswärtigen ordentlichen Mitglied ihrer physi- 
kalisch-mathematisehen Klasse bestätigt. 

Der Ernst Abbe-Gedenkpreis der Carl Zeiß- 
Stiftune wurde Prof. Dr. Wilhelm Blaschke- 


Hambure verliehen. 

Dr. phil. habil. Günther Sehulz- Berlin wurde 
die Dozentur für reine und angewandte Mathe- 
matik in der mathematisch-naturwissenschaftlichen 
Fakultät der Universität Berlin verliehen. 
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